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Tne .chemical  and  vibrational  kinetics  of  the  H2-F2  chain  reaction  were  studied 
to  imorove  the  description  of  the  evolution  of”the  various  vibrational  levels 
in  order  to  better  predict  chain-reaction-laser  performance.  Two  types  of 
experiments  were  performed:  (1)  flew  tube,  and  (2)  oulsed  reaction'.  In  the 

now  tube  experiments,  fluorine  atoms  were  crocuced  by  passing  through  a 
mass^heater.  ^The  F2  and  F emanating  from  the  heater  were  mixed  with  additional 
F£\  Ar,  and  H2  to  achieve  the  desired  stoichicmetry .xThe  concentrations  (OVER) 
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20.  ABSTRACT 

of  the  HF  vibrational  levels  from  0 through  6 were  monitored  as  the  mixture 
flowed  through  a 5 cm  diameter  tube  through  the  use  of  a gain  probe,  fundamen- 
tal band  IR  emission  spectroscopy  and  Lys 3 overtone  emission  spectroscopy.  In 
the  pulsed-reaction  experiments  a predetermined  mixture  of  F2,  H2  and  Ar  was 
sealed  in  a small  constant-volume  test  cell.  The  reaction  was  initiated  by 
the  dissociation  of  a small  fraction  of  the  F2  with  a frequency  doubled  pulsed 
ruby  laser.  Pressure  changes  were  monitored  with  a Kistler  gauge  and  ground 
state  HF  with  a probe  laser.  Vibrational  levels  1 through  6 were  observed  by 
fluorescence  measurements.  In  addition,  gain  was  measured  on  the  1-0  and  2-1 
transitions  and  laser  tests  performed  which  achieved  lasing  as  high  as  the 
5-4  transition. 

The  results  of  the  study  indicate  no  significant  changes  need  be  made  to  the 
rates  normally  used  for  the  hot  and  cold  reactions.  Vibrational  deactivation 
rates  must  be  scaled  up  rapidly  with  quantum  number,  however,  to  adequately 
describe  the  measured  distributions.  A viable  deactivation  model  Is  presented 
and  possible  improvements  discussed.  It  appears  that  the  HF  chain  reaction  is 
of  limited  usefulness  as  a cw  laser  device. 


UNCLASSIFIED 


SUMMARY 


In  order  to  better  understand  the  physical  mechanisms  important  in  chain 
reaction  lasers  and  their  effect  on  performance,  a basic  study  was  initiated, 
comprised  of  identification  of  the  problem  areas,  experimental  testing,  and 
subsequent  improvement  of  theoretical  models. 

The  major  problem  areas  were  identified  as  uncertainty  of  the  values  for 
and  effects  of  (l)  the  chemical  pumping  reaction  rate  coefficients  (H  + F2  -» 

HF  + F,  F + Hp  -»  HF  + H) , (2)  HF-KF  vibrational  deactivation,  and  (3)  H atom 
deactivation  of  HF.  Potentially  significant  problem  areas  also  included  wall 
effects,  back  reactions,  hot-atom  effects,  and  chain  branching. 

The  laboratory  experiments  to  study  these  problem  areas  fell  -under  two 
categories:  cw  flow-tube  experiments  and  stationary-gas  pulsed-react ion 

experiments . 

The  cw  flow-tube  chain-reaction  experiment  was  initiated  by  the  thermal 
dissociation  of  small  quantities  of  F.p,  subsequently  mixed  with  F0,  Eg,  and 
Ar  in  a flowing  situation.  As  the  reaction  progressed  through  the  flow  tube 
population  distributions  of  the  various  v and  J levels  were  measured. 

Initial  F- atom  concentrations  were  measured  with  metered  HC1  titration 
and  HF  probe  laser  absorption.  This  technique  was  verified  by  separate  electron- 
spin-  resonance  tests  performed  under  the  subject  contract.  Initial  temperature 
profiles  were  measured  with  a thermocouple  probe,  and  velocity  profiles  with  a 
pitot  probe.  These  results  were  incorporated  in  the  analysis.  Ground  state  HF 
was  measured  with  probe  laser  absorption,  whereas  for  vibrational  level  popu- 
lations in  v = 1,  2,  3,  ^ absolute  IF.,  Av  = 1,  emission  spectroscopy  was  used. 
Vibrational  levels  v = 3>  k,  5,  6 were  determined  with  overtone  emission 
spectroscopy  using  hv  = 3* 

Tests  were  performed  over  a range  of  &,/ F„  ratios  with  a fisted  initial 
F-atom  concentration  in  order  to  obtain  reaction  rate  and  vibrational 
deactivation  information.  Some  HF  was  added  initially  to  enhance  self- 
deactivation  effects.  Subsequent  tests  involved  variation  of  the  initial 
F-atom  concentration  while  trying  to  maintain  the  same  net  HF  production  rate 
in  order  to  isolate  the  influence  of  H-HF  collisions. 

The  pulsed  reaction  experiments  were  used  to  study  H-.-F'  chain  reaction 
kin  tiCfl  in  a :losed  volume  containing  premised  Hp  and  F- . The  reaction  was 
initiated  by  producing  fluorine  atoms  in  the  mixture  by  dissociating  ?-  with 
pulsed  'J1/  radiation  from  a doubled  frequency  ruby  laser.  Three  types  of 


experiments  were  performed:  First,  a study  was  made  of  the  temporal  development 

of  populations  in  vibrational  states  of  HF  produced  from  the  reaction.  Fluo- 
rescence measurements  were  used  to  track  excited  state  populations  and  absorption 
of  cw  radiation  to  follow  ground  state  formation.  Second,  small  signal  gain 
was  determined  for  v(l-C)  and  v(2-l)  transitions  by  use  of  a gain  probe.  Third, 
pulsed  laser  radiation  from  the  reacting  mixture  was  studied  using  a laser 
configuration. 

The  rationale  for  these  experiments  was  that  they  would  serve  as  an 
experimental  check  on  the  laser  computer  program  as  modified  by  results  from 
flow-tube  experiments  and  other  available  data.  Pulsed  experiments  were  chosen 
as  being  relatively  simple  in  that  the  problems  attendant  with  mixing  and  flow 
are  not  present,  the  reaction  is  initiated  over  the  total  volume  in  a short 
dme  period  and  changes  in  temperature  and  pressure  could  be  monitored  in  a 
constant  volume  cell. 


■Analysis  of  the  data  from  the  two  sets  of  experiments  has  indicated  that 
no  drastic  changes  need  be  made  to  the  chemical  reaction  rates  now  in  use.  The 
cold  reaction  rate  was  found  to  lie  between  1.0  and  2.0  times  the  nominal  value, 
and  the  hot  reaction  rate  between  0.5  and  1.3  times  its  nominal  value.  Other 
effects  (back  reactions,  chain  branching,  hot-atom  enhancement)  could  not  be 
easily  isolated,  but  appear  to  be  of  minor  importance.  Under  some  circumstances 
(such  as  these  experiments)  wall  effects  can  be  nornegligible. 

The  original  vibrational-deactivation  model  was  found  to  be  inadequate 
in  describing  upper  level  vibrational  populations.  Much  better  agreement  has 
been  achieved  by  adopting  a model  allowing  for  rapid  deactivation  of  upper 
vibrational  levels  via  multiquanta  V-T  transitions.  Single  quantum  V-V 
exchange  is  then  found  to  be  somewhat  slower  than  expected,  but  in  any  case  to 
weakly  influence  the  population  distributions.  Remaining  deficiencies  in  the 
fit  of  the  model  to  the  data  indicate,  however,  that  multiquanta  V-V  exchange 
may  be  the  key  to  even  better  agreement.  Such  a possibility  should  be  explored 
f -.rthe  r . 

The  direct  influence  of  H atoms  on  vibrational  deactivation  proved 
difficult  to  isolate  due  to  the  strong  influence  of  HF  self-deactivation, 
which  is  still  not  completely  understood,  ir.  the  test  conditions  feasible  in 
these  experiments.  Hence,  no  new  conclusions  were  reached  with  regard  to  such 
rates,  now  believed  to  be  very  rapid  for  v £ 3- 


The  laser  model  as  modified  was  employed  to  predict  the  performance 
potential  of  chain-reaction  lasers.  It  was  fo'ind  that  operation  in  a fluorine- 
rich  condition  tended  to  give  the  best  gains  and  powers,  in  agreement  with 
experimental  observation.  Nonetheless,  the  enhanced  vibrational  deactivation 
was  suffici'  nt  to  reduce  predicted  output  power  to  the  extent  that  chain-reaction 
las  not  competitive  with  cold-reaction  lasers  for  high  altitude  applicati  ns. 
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SECTION  I.  INTRODUCTION 


The  work  described  in  this  report  covers  the  research  carried  out  under 
Phase  I of  the  subject  contract.  Phase  I was  a technology  development  pro 
designed  to  improve  the  ability  to  predict  analytically  the  performance  po 
of  chain  reaction  lasers,  and  to  assess  their  feasibility,  particularly  with 
respect  to  high  altitude  missions.  The  program  comprised  three  tasks:  (l) 

identification  of  kinetic  problem  areas  in  model  predictions,  (l)  experimental 
testing  to  study  those  problem  areas,  (3)  analysis  of  the  data  and  resultant 
modification  of  existing  computer  models  to  assess  performance  potential. 

The  details  of  the  various  experiments,  the  subsequent  analysis  and  cede 
modification,  and  the  resulting  conclusions  will  bo  presented  below. 


c+  (It) 


SECTION  II 


CW  PLOW  TUBE  REACTOR  EXPERIMENT 


Introduction 

The  CW  flow  tube  reactor  was  selected  as  one  of  the  experiments  to  be 
conducted  because  populations  in  the  various  V and  J levels  and  the  ground 
state  are  readily  obtained  under  steady  state  conditions  (Refs.  1-3)*  Results 
of  Task  I,  Problem  Area  Identification  (Ref.  4)  indicated  that  the  chemical 
and  relaxation  processes  probably  account  for  a major  share  of  the  observed 
discrepancy  between  predictions  and  experiment.  Results  of  Task  I also  indi- 
cated a practical  range  of  conditions  over  which  flow  tube  results  likely  could 
provide  rate  information  for  the  overall  reaction  and  molecular  energy  trans- 
fer processes  that  are  occurring.  A discussion  of  the  flow  tube  research 
conducted  to  obtain  this  information  is  presented  in  the  following  sections. 

In  general,  significant  data  results  on  vibrational  level  populations  and 
total  HE  formation  are  included  in  the  discussion  under  Section  IV  of  this 
report . 

Description  of  Experimental  Configuration 

The  apparatus  components  employed  in  the  tests  were  essentially  those 
described  in  References  1 through  3 with  suitable  modifications  of  a miner 
nature  to  achieve  additional  required  data  and  improved  accuracy  of  results. 

A schematic  of  the  major  components  as  employed  here  is  shown  in  Figure  1. 

Fluorine  mixed  with  a constant  purge  flow  of  argon  makes  three  passes  in 
a nickel  tube  assembly  through  the  heated  zone  of  a well  regulated  commercial 
tube  furnace.  At  low  pressure  significant  dissociation  of  ?2  (ujx  = 1 — ' 
percent)  can  occur  thermally  (Ref.  3)  at  temperatures  under  lGCC°n,  a safe 
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3.  Bcedeker,  L.  R. , J.  F.  Verdieck,  and  R.  J. 
Final  Technical  Report,  Contract  N000l4-7u 
R75-951S83-4 , East  Hartford,  CT,  May  1975 • 


Hall,  "HF  Chain  Reaction  Laser," 
C-0379  (HRL  MS).  UTRC  Report 
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FIGURE  1 HF  CHAIN  REACTION  FLOW  TUBE  EXPERIMENT 


working  level  for  nickel.  The  hot,  partially  dissociated  F2  passes  rapidly 
out  of  the  heated  zone  through  an  array  of  seven  holes  and  the  translational 
temperature  is  reduced  rapidly  to  near  room  temperature  by  mixing  with  argon 
introduced  through  a fine  scale  injector  ring  placed  close  to  the  end  of  the 
heated  length.  This  injector  is  supported  by  a length  of  copper  pipe  which 
is  water  cooled  (not  shown  in  Fig.  l)  and  does  not  touch  the  hot  nickel  tube 
assembly  except  at  the  flange  joint  3C  cm  downstream  of  the  heated  zone  and 
outside  the  tube  heater  where  all  components  are  water  cooled.  The  trans- 
lational temperature  of  the  F0  is  lowered  rapidly  enough  here  to  preserve  ? 
atoms  in  flow  from  the  mass  heater  to  ohe  H2  injector  (Refs.  1-3 ).  In  order 
to  adjust  the  overall  level  of  dissociation,  additional  Fg  is  added  to  the 
argon  diluent  flow  which  thus  bypasses  the  mass  heater. 

An  auxiliary  injector  ring  was  located  in  the  F2  dilution  section 
through  which  metered  flows  of  EF  or  HC1  were  added  upstream  of  the  Ep 
injector. 

Before  the  Ho  injector,  a transition  section  was  placed  to  adjust  the 
diameter  of  the  flow  smoothly  to  the  diameter  of  the  reaction  zone  d = 5 cm. 

The  Kj  mixer  is  an  array  of  small  tubes.  Figure  2,  with  small  injection  holes 
angled  30°  downstream  to  avoid  impingement  of  the  injectant  or.  adjacent  tubes 
and  baekstreaming  effects.  The  K2  is  mixed  with  a large  amount  of  argon  and 
large  penetration  occurs.  The  characteristic  mixing  tends  to  be  between 
layers  separated  by  about  1.5  mm,  the  hole  spacing  along  each  tube,  thus 
assuring  mixing  in  1 to  2 cm  (Refs.  l-a). 

Viewing  ports  were  placed  along  the  reaction  tube  a z x = 4,  3,  12,  24,  31, 
bl,  51  cm  from  the  H2  injector.  One  viewing  port  was  placed  1 cm  upstream  of 
the  E;  injector  to  check  for  baekstreaming.  Each  set  of  viewing  ports  was 
purged  with  a metered  flow  of  argon.  The  viewing  ports  contained  narrow  slots 
cut  through  to  the  reaction  zone.  Thus,  maximum  spontaneous  emission  could  be 
collected  by  aligning  the  spectrometer  slit  with  the  long  side  of  the  slot 
while  minimizing  distortion  of  the  flow  velocity'  profile. 

The  inside  of  the  reaction  tube  and  transition  were  coated  with  Teflon 
to  inhibit  H atom  recombination  on  the  wall  (Ref.  5)-  The  particular  coating 
that  gave  best  application  results  contained  about  13  percent  TiSo  as  pigment. 
Substantial  reduction  in  E anom  recombination  would  be  expected,  however, 
compared  to  pure  copper  based  on  exposed  surface  area  of  Teflon.  Present  tests 
are  also  at  a factor  of  4 higher  pressure  than  previously  (Ref.  3)  reported. 
Experimental  results  must  b4  examined  closely,  however,  to  see  if  an  H atom 


Wood,  R.  J.  and  H.  Wise,  ''Kinetics 
Surfaces."  J.  Them.  Phys . , 65.  197b 


of  Hydrogen  Atom  Recombination  or 
. 19c 1. 


ratios . 


or  other  wall  effect  exists  at  large  dis 
calculations  and  examination  of  results 
effects  to  be  a factor  unless  the  H atom 
would  be  limited  to  25  percent. 


tances  or  Fp/Ho  ratios.  At  present, 
in  Section  17  do  not  indicate  wall 
loss  was  unity  and  then  the  effect 


Test  Program  and  Instrumentation  for  CW  Plow  Tube  Reactor  Experiment 

Flow  tube  experiments  were  directed  towards  (l)  acquisition  of  F atom 
titration  data  with  HC1  addition  and  HF  probe  laser  absorption  in  order  to 
calibrate  the  mass  heater  and  verify  establishment  of  desired  low  levels  of 
initial  F^  dissociation  and  (2)  acquisition  at  several  stations  along  the 
flow  tube  reaction  zone  of  ground  state  and  vibrational  level  population 
information  from  probe  laser  absorption;  dv  = 1 and  hv  = 3 emi  ■ ion  spectros- 
copy data  for  conditions  determined  from  Task  I studies  in  order  to  provide 
relevant  fits  to  the  model  that  would  yield  improved  information  on  rates  for 
the  hot  reaction;  HF  V-7,  V-T  self-relaxation;  and  K atom  deactivation.  In 
initial  work,  the  Proposed  Research  Plan  (Ref.  h)  for  flow  tube  experiments 
was  prepared,  and  based  on  the  results  of  Task  I ani  discussion  with  AFWL  and 
Aerospace  personnel,  increased  emphasis  was  placed  on  verification  of  velocity 
and  temperature  uniformity  with  most  testing  conducted  at  higher  press-ure,  up 
to  p = 20  torr,  at  velocity  v ~ lC"1  cm/sec  with  a useful  reaction  zone  length 
L ~ 50  cm.  A summary  of  instrumentation  is  given  in  Table  1.  The  overall 
arrangement  of  the  instrumentation  and  the  apparatus  is  shown  in  Figure  3. 
Detailed  discussion  of  each  diagnostic  and  results  is  given  after  the  overall 
conditions  that  were  attained  are  summarized  in  the  next  section. 


Summary  of  Conditions 


The  objectives  of  the  experiment  were  sought 


two  series  of  tests: 


(l)  In  the  first  series  operation  over  a wide  range  of  Fp/EA  about  unity 
was  sought  maintaining  fixed  F in  order  to  obtain  information  about  the  hot 
and  cold  reaction  rates  near  room  temperature.  Conditions  attained  are 
summarized  in  Table  2.  This  series  also  provided  HF:KF  interaction  informa- 
tion since  HF  production  varied  over  a wide  range  and  generally  KF:HF  appeared 
to  dominate  H:HF  interaction  where  effects  were  measurable.  Such  EF:EF 
information  was  supplemented  by  direct  addition  of  HF  to  the  flow  tube  for 
these  cases. 


(2)  In  the  second  series  operation  was  sought  at  a higher  level  of  F and 
?2  iissociatior.  maintaining  constant  Ho  and  nearly  constant  HF  production. 
Such  a condition  was  sought  by  decreasing  the  level  of  proportional  to  the 
Lncr  tuse  in  F.  The  relative  effect  of  K:KF  collisions  is  enhanced  as  a resull 


i 
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TABLE  1.  INSTRUMENTATION 

1.  Probe  Laser  Absorption:  v = 0,  initial  F 

2.  Overtone  Spectroscopy  Av  = 3>  SMI  S-l  Fhotodetector  (T  = -50°C 

controlled):  v - 2, , k , 5»  6 

3.  IR  Spectroscopy  Av  = 1,  FbS  Detector  (Ul0  cooled):  v = 1,  2,  3>  *+ 

Pitot  and  wall  pressure  distributions:  gas  velocity  and  pressure 

5.  Gas  Temperature  Radial  Distribution 

6.  Gas  feed  rates,  inert  argon,  Eg  and  corrosive  Fp,  HF,  EC1:  species 
concentrations . 
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TAELS  2.  SUMMARY  OF  COHDITIOIIS* 


I)  Chemical  Rates  and  HF  additive 

p =20.0  ± O.U  torr  T.  = 308  ± L°K 

v = ( graph  given)  + 0 

- 6 i 


Initial  Conditions 


Case 

- 

f2~ 

scc/min  (torr) 

h2 

scc/min  (torr) 

Ar 

sl/min 

F/2F2 

7 

0.0021  ± 

.002 

1130  ± Lo 

(.0950) 

135  ± 3 
(.0113) 

238  ± 12 

0.011 

2 

0.0021  ± 

.002 

6U0  ± 30 
( .0538) 

200  ± 5 
(.0168) 

238  ± 12 

0.020 

3 

0.0021  ± 

.002 

250  ± 15 
(.0210) 

UC0  ± 10 
( .0336) 

238  ± 12 

0.050 

L 

0.0021  ± 

.002 

60  ± 5 
( .0050) 

800  ± 20 
(.0672) 

238  ± 12 

0.210 

♦These  conditions  are  the  average  values  and  maximum  uncertainty  of  quan- 
tities that  were  recorded  as  each  case  was  repeated  for  the  separate  probe 
laser,  over*one  emission,  IR  emission  and  F atom  titration  measurements. 

♦■♦Total  ^2  feed  to  ex;  eriment 
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FIGURE  3 FLOW  TUBE  REACTOR  EXPERIMENTAL  ARRANGEMENT 
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Actual  initial  conditions  attained  for  the  two  resulting  cases  are  summarize i 
in  Table  5-  Results  for  HF  production  and  vibrational  level  changes  are  given 


in  j^ction  iv. 


Vacuum  Purnpii:.-  arrangement 


Expe  r ime  r.t  and  Results 


The  experiment  was  connected  through  a 1 in.  vacuum  line  (with  in-line  11b 
trap.,  isolation  valve  and  throttl0  valve)  to  a 3C0C  Cr?4  mechanical  blower  backed 
by  two  ISO  CFM  trochiodal  p’jmps . 

Retting  of  Van  Flows 

Diluent  argon  for  3 £m_d_iameter  f low_tube  at_p_=_2f  torn 

Argon  was  selected  as  the  diluent  for  the  experiment  instead  of  helium 
because  its  higher  molecular  weight  at  constant  velocity  results  in  higher 
Reynolds  number  flow,  thinner  boundary  layers  and  a more  uniform  velocity 
profile  in  the  reaction  zone.  'This  minimizes  corrections  that  have  to  be 
applied  in  the  analysis  to  account  for  diffusion  of  HF(v)  to  the  walls  and 
for  streamlines  in  the  flow  having  different  velocities  and  time  histories. 
Uniform  velocity  is  enhanced  for  the  present  flow  tube  since  the  pressure  and 
velocity  are  moderately  high,  p ~-20  torr  and  v ~ 1C4  cm/sec.  The  combination 
of  these  conditions  permits  operation  here  with  a tube  diameter  small  relative 
to  the  reaction  zone  length  D/L  = 0.1. 


A listing  of  the  argon  feed  distribution  along  the  flow  tube  is  given  in 
Table  1:  (l)  Argon  flow  to  the  mass  heater  is  needed  to  overcome  a slight 

ejector  effect  at  the  Fp  dilution  injector  and  maintain  the  pressure  in  the 
mass  heater  higher  than  in  the  F2  dilution  stage.  This  assures  uniform  dis- 
tribution in  the  mass  heater  ?2  effluent  in  the  seven  houe  array  (Fig.  l)  in 
he  presence  of  small  radial  pressure  gradients  at  the  ?-  dilution  injector, 
uch  purge  also  serves  to  clear  the  system  of  F?  rapidly  when  establishing 
reference  levels  for  absorption  measurements,  (2)  The  argon  flow  to  the  F- 
dilution  injector  is  limited  at  present  by  a pressure  limit  in  the  Fp  system 
which  supplies  additional  cold  F-  for  mixing  with  this  diluent.  (;)  The 
argon  flow  in  the  He  injector  is  higher  than  in  the  F?  dilution  injector  “:o 
max:-  cp  for  this  limit  in  ?2  dilution  injector  flow.  Also,  velocity  profile 
tests  gave  slightly  more  favorable  results  with  more  argon  in  the  Hr  inj-'-ctor. 
(4)  Purge  in  the  annular  dead  zone  between  the  Fo  dilution  injector  ar.d  nickel 
mass  heater  was  thought  to  be  necessary,  but  had  to  be  k-  pt  to  a low  value  to 
maintain  temperature  uniformity  as  discussed  later  in  this  section.  (c) 

Steady  purge  of  the  HF(HCl)  injector  assured  attainment  of  a fast  probe  laser 
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TABLE  3.  SUMMARY  OF  CONDITIONS* 


II)  H Atom  Deactivation 


p = 2 ±0.4  torr  Ti  = 308  ± 4°K 


V 


(graph  given) 


+ 0 
- 6* 


Case 


F 

torr 


Initial  Conditions 


f2 

see 


/min 


h2 

scc/min 


Ar 

sl/ain 


5~* 


0.0019  ± .0001 


4 10  ± 20 
0.034 
±10*0 


210  ± 10 
0.013 
±10fo 


238  ± 12  .028 


0.0052  ± .0003 


140  ± 6 
0.012 
±10$ 


210  ± 10  238  ± 12  .217 
0.018 
±10$ 


♦These  conditions  are  the  average  values  and  maximum  uncertainty  of  quan- 
tities that  were  recorded  as  each  case  was  repeated  for  the  separate  probe 
laser,  overtone  emission,  IR  emission  and  F atom  titration  measurements. 

♦*F2  feed  to  mass  heater  ♦ F2  feed  bypassing  mass  heater 

♦♦♦Reduce  HF^v=l)  from  Av  - 1 spectra.  Since  F-  hot  and  cold  feeds  were  each 
about  5 percent  too  high.  Hence  F^Ho  product  was  10  percent  too  high  for 
that  data. 
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TA£LE  4.  AEGON  FEED  DISTRIBUTION  TO  FLOW  TUBE 


Location 

1 Mass  Heater  Inlet 

2 ?2  Dilation  Injector 

3 H2  Injector 

4 Purge  between  mass  heater  and 
Fg  dilution  stage 

5 Steady  purge  of  HF(HCl)  Injector 

6 Total  viewpoint  purge 


?a~. e calibrated  Rotometers) 
std  liters/min 

3 

63  ± 3 
146  ± 5 
3 


233  ± 12 


zero  reference.  Such  purge  flow  could  not  be  much  higher  than  the  HF  (and 
probably  HCl)  flow  rates  in  order  to  achieve  steady  and  accurate  HF  absorption 
signals  due  to  the  additive;  this  condition  is  probably  indicative  of  an  unstable 
chemical  wall  effect  in  the  copper  feed  lines  at  high  dilution.  (6)  Purge  to 
each  set  of  viewports  was  controlled  separately  with  needle  valves.  The  total 
flow  to  all  eight  sets  is  given  in  Table  4.  The  purge  flow  through  each  port 
into  the  flow  tube  represents  a 0.0036  fraction  of  the  total  flow  passing  by  the 
slot  inside  the  flow  tube;  a small  disturbance  to  the  velocity  profile.  Such 
small  purge  proved  adequate,  however,  to  remove  ground  state  HF  from  the  window 
volume  and  establish  stable  levels  of  absorption  insensitive  to  changes  in 
purge  flow  about  this  value. 

Reactive_Gases_and_Mditive  3 

Gases  Fg,  HF  and  HCl  were  metered  with  monel  mass  flow  transducers 
obtained  from  Matheson.  These  flow  transducers  were  calibrated  and  checked 
on  a regular  basis  using  a standard  volume  in  the  feed  system,  isolation 
^®-T.es,  a preoSure  gauge  and  a stopwatch.  For  H?  (UHP  grade)  a rotcmeter 
«.ave  reliable  calibration  results  at  the  feed  pressure  employed  and  was 
preferred  over  the  slow  responding  mass  flow  transducer.  Commercial  ?2 
was  used  directly  from  LA.  cylinders  with  a UTRC  safety  enclosure  system:  an 
in-line  sodium  bifluoride  trap  was  used  to  remove  HF  from  F£ . Distilled 
HF  was  obtained  in  a monel  cylinder  from  Ozark  Mahoning  and  was  cooled  with 
LH2  and  pumped  on  prior  to  use;  HCl  was  used  directly  from  a #3  cylinder 
of  technical  purity. 

Pressure  Distribution  in  Flow  Tube 

The  pressure  data  shown  in  Figure  4 include  repeatability  obtained  from 
tests  on  separate  days.  Data  were  obtained  with  a wafer  fluid  scanner  and  a 
pressure  transducer.  Reactants  have  only  a small  effect  on  the  flow  pressure 
for  these  very  dilute  cases.  As  discussed  above,  the  pressure  in  the  mass 
heater  is  maintained  at  a level  above  the  value  in  the  F2  dilution  section. 

A small  pressure  rise  &p  ~ 1 torr  occurs  across  the  H2  injector  which  indicates 
that  a very  small  ejector  effect  exists  with  this  angled  injector  (Fig.  2). 

The  pressure  in  the  reaction  zone  is  relatively  flat  with  an  average  value  of 
about  20.0  torr.  (Calibration  of  pressure  was  with  a Wallace  and  Tiernen 
gauge . ) 

Velocity  Profiles  Along  Flow  Tube,  Check  on  Velocity  Results  With  Metered 
Argon  Flow 

Velocity  profiles  at  several  stations  along  the  flow  tube  were  measured 
with  a l/l6  in.  diameter  pitot  probe  and  a MK5  Baratron  differential  transducer. 
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FIGURE  4 PRESSURE  DISTRIBUTION  ALONG  FLOW  TUBE 


A radial  traverse  assembly  was  used  as  shown  in  Figure  5 that  could  be  installed 
at  each  set  of  viewports.  Data  were  obtained  in  nonreactive  flow.  Pitot 
pressure  results  Ap  ~ 1.5  torr  maximum  corresponded  to  a subsonic  core  Mach 
number  in  argon,  M ~ 0.3*  Probe  results  were  converted  to  velocity  using 
standard  compressible  flow  tables.  The  probe  was  long  enough  in  the  flow 
direction  and  alignment  was  accurate  enough  such  that  little  probe  correction 
is  needed  and  none  has  been  applied.  A check  was  made  with  argon  feed  rate 
for  one  case  as  discussed  below. 

Results  at  x = 11,  20,  40  cm  are  shown  in  Figure  6 and  at  x = 76  cm, 

Figure  7-  Results  were  very  steady  and  repeatable.  In  general,  the  profiles 
were  as  uniform  as  desired  and  steepest  velocity  gradients  were  within  one 
probe  diameter  of  the  wall.  Close  to  the  injector  the  data  show  a local 
peaking  in  velocity,  probably  a complex  characteristic  of  the  details  of  the 
mixer  and  flow  tube  geometries.  Since  such  a perturbation  could  be  accounted 
for  by  incorporating  this  profile  in  the  analysis,  the  effort  required  to 
correct  it  by  altering  the  geometry  has  not  been  made.  The  probe  was  subse- 
quently located  at  x = 76  cm  downstream  of  the  reaction  zone  of  interest, 
but  where  the  tube  diameter  was  still  D = 5 cm.  Velocity  results  still  show 
a large  uniform  core  flow.  The  average  velocity  obtained  by  a ten  step  inte- 
gration over  this  radial  profile  was  0.88  x 104  cm/sec  which  for  p = 1.96  torr, 
the  pressure  at  x = j6  cm  corresponds  to  a flow  rate  of  253  std  liters /min, 
very  close  to  the  metered  flow  rate  240  std  liters /min. 

Gas  Temperature  Uniformity 

Gas  temperature  was  measured  with  a simple  thermocouple  probe  inserted 
through  a port  in  the  transition  section  about  4 cm  upstream  of  the  Hq  injector 
(Fig.  l).  With  the  aid  of  relations  in  Reference  6,  the  error  introduced  by 
probe  conduction  was  estimated  to  be  only  2°C,  a small  value  here.  Operational 
experience  has  shown  that  temperature  uniformity  is  influenced  greatly  by  the 
flow  purging  the  annular  region  between  the  mass  heater  and  ?2  dilution  stage. 
This  purge  enters  near  the  downstream  end  of  the  Fo  dilution  stage  where  the 
mass  heater  wall  is  cold,  flows  in  the  upstream  direction  and  enters  the  main 
flow  through  the  small  gaps  between  the  seven  hot  exhaust  ports  in  the  mass 
heater  and  the  cold  Fo  dilution  injector  (Fig.  l) . Temperature  studies, 

Figure  8,  show  that  this  purge  flow  must  be  low  in  order  to  preserve  the  tem- 
perature uniformity  provided  by  high  F_  dilution  and  :old  dilution  stage  liner 
vails . 


6.  Rizika,  J.  W.  and  W.  M.  Rohsenow,  "Thermocouple  Thermal  Error," 
Industrial  and  Engineering  Chemistry,  kU,  1168,  1952 • 
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FIGURE  8 EFFECT  OF  PURGE  GAS  FLOW  RATE  ON  FLOW  TUBE  TEMPERATURE  DISTRIBUTION 


With  the  dilution  purge  flow  set  at  a low  value  the  effect  of  the  mass 
heater  flow  on  the  gas  temperature  distribution  was  investigated  by  turning 
the  mass  heater  flow  on  using  argon  instead  of  Results  (Fig.  9)  show  a 

relatively  uniform  increase  of  1°C  across  the  flow  tube  due  to  the  mass  heater 
which  indicates  reasonably  uniform  midng  has  occurred.  For  actual  conditions 
Tables  2,  3>  and  4 > this  gas  temperature  measurement  rose  to  about  62°C  (335°K) 
and  uniformity  while  not  checked  was  probably ± 4°c.  When  mixed  with  room 
temperature  argon  in  the  Ho  injector,  the  flow  tube  initial  temperature  is 
estimated  to  be  T = 3C8  - 3°K. 

HR  Probe  laser  Absorption 

Measurement  Technique 


An  H?  probe  laser  absorption  diagnostic  (Fig.  10)  was  assembled  under  the 
present  program  employing  a CW  HF  single  mode  probe  laser  similar  to  the 
one  described  in  Reference  7.  The  probe  laser  was  mounted  on  an  optical  table 
isolated  with  an  air  suspension  system  with  well  regulated  table  height.  The 
sensitive,  stable,  double  beam,  single  detector  arrangement  of  Figure  10  was 
selected  based  on  the  expected  requirements  for  the  F atom  titration  measuremen 
A large  area  room  tempera-cure  FbS  detector  coupled  with  a 5 cm  diameter 
aluminum  sphere  provided  adequate  signal  and  was  relatively  insensitive  to 
small  angular  perturbations  in  beam  position.  The  overall  arrangement  of 
components  on  the  4 ft.  x 12  ft.  optical  table  is  shown  in  Figure  3.  When 
driving  the  PZT  stack  on  the  laser  cavity  with  a ramp,  a stable  power  profile 
over  the  doppler  broadened  gain  region  was  achieved  as  viewed  on  che  FbS  detec- 
tor (Fig.  ll).  A Lamb  dip  is  clearly  evident,  further  indication  of  a stable 
single  longitudinal  node.  In  general  for  absorption  measurements  here,  the 
cavity  PZT  voltage  was  set  manually  near  line  center;  the  output  frequency 
would  remain  within  the  Lamb  dip  region  for  periods  of  minutes  without  active 
feedback.  Oene rally,  frequency  was  stable  to  ±30  MHz,  sufficient  for  about 
3 percent  accuracy  in  absorption  (Ref.  7)- 

The  operational  characteristics  of  this  dual  beam  absorption  technique 
with  a single  detector  are  shown  in  Figures  12  and  13 . Note  in  Figure  13  that 
stability  of  better  than  1 percent  of  full  scale  was  achieved  for  1-minute 
periods  (typically  stability  was  0.-+  percent).  Such  stability  allowed  measure- 
ment  of  F atom  titration  signals  to  about  1C  percent  accuracy  for  the  lowest 
?2  feedrates  of  interest  into  the  mass  heater,  ~ 3C  sec,  min  (Fig.  13).  Not? 
also  that  the  window  purge  flow  was  on  and  set  as  discussed  in  the  section  on 
ras  flow  rates . 


7.  Hinehen,  J.  <7.,  "Determination  of  Vibration-Rotation  Line  Strengths  for 
HF  and  DF  by  Use  of  an  HF.  DF  CW  Las-  r, " 7.  .dw  . Soc«  Am. , • - . ■ . I'm, 
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FIGURE  EFFECT  OF  MASS  HEATER  ON  FLOW  TUBE  TEMPERATURE  DISTRIBUTION 
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FIGURE  13  PROBE  LASER  ABSORPTION  DATA  ACQUISITION 


Calibration  Check  With  Direct  Hr  Addition 


Results  ox'  tests  in  which  HF  was  added  to  the  flow  tube  through  a separate 
injector  ring  (Fig.  l)  are  shown  in  Figure  lh.  These  measurements  are  pretty 
much  in  the  linear  range  ox’  the  absorption  and  agree  very  well  with  calculations 
using  Hinchen's  data  for  doppler  broadened  line  center  absorption  coefficient 
(Ref.  ?)  with  the  10  cm  total  absorption  length  in  this  double  pass  arrangement 
(Fig.  10).  Pressure  broadening  due  to  Ar  does  not  se^m  to  be  a factor  her®  at 
p = 20  torn  and  based  on  results  of  pulsed  experiments,  Section  III,  would  be 
at  most  a 5 percent  effect.  Such  data  (?ig.  1+)  suggest  also  that  there  is 
little  loss  of  ground  state  HF  to  the  walls  between  the  HF  injector  and  obser- 
vation point  x = ll  cm.  HF  diffusion  would  be  slower  than  H and  even  for  total 
disappearance  at  the  wall  would  probably  be  limited  to  a 5 to  10  percent  effect. 

Titration  for  F Atoms 


The  HC1  titration  technique  x’or  F atoms  was  checked  out  with  F3R  tests 
as  part  of  this  program  (See  Section  V).  This  technique  was  applied  in  the 
flow  tube  studies  by  detecting  the  appearance  of  HF  with  the  probe  laser  as 
HC1  was  added.  At  first  excess  HC1  (>  1Q0C  scc/mir.)  was  added  to  the  argon 
flowing  to  the  Eq  injector  with  R;  = 0 (Fig.  13).  A saturation  level  of  HF 
formed  was  noted  well  downstream  and  checked  for  repeatability  for  each  run. 
This  absorption  saturation  level  varied  with  F ^ feed  to  the  mass  heater  as 
expected  (Fig.  13)  and  at  the  expected  levels.  Very  little  Interaction  of 
this  saturation  level  with  cold  Fp  flow  was  noted  with  excess  HC1,  consistent 
with  the  slow  Cl  + ?2  reaction  rate  at  room  temperature  (see  Section  V). 


Since  some  slight  changes  in  the  excess  HC1  titration  result  (£l0  percent) 
were  noted  periodically  with  cold  F- , an  attempt  was  made  to  check  the  stoi- 
chiometry of  the  titration  reaction  with  metered  HC1  addition.  In  order  to 
accomplish  this,  the  HC1  was  added  through  the  injector  ring  that  had  been 
used  x’or  HF  (Fig.  l)  since  the  concentration  of  HC1  in  the  x’eed  lines  to  this 
ring  could  be  kept  high  to  avoid  having  difficulties  with  spurious  loss  of 

to  be  a concern  in  KF  additive  studies, 
first  a 1:1  HF:KC1  correspondence  does 
;'e  occurs  x’rom  such  a linear  response 
ith  HC1  flow  rate  thereafter.  This 
*■  F?  reaction.  This  system  has  not  been 
r now  it  has  seemed  reasonable  to  extra- 
its intercept  with  the  linear  respons 3 
curve  to  define  the  level  of  F atoms.  This  point  for  cases  1 through  - of 

town  in  Figure  15  while  ? levels,  cases  5 and  o of  Table  TIT,  are 
hown  in  Figure  it . An  overall  calibration  of  the  mass  heal  sr  from  such  int  r- 
cepts  is  given  in  Figure  17. 
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FIGURE  14  PROBE  LASER  ABSORPTION  CALIBRATION  CHECK  BY  HF  ADDITION 
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FIGURE  15  HCI  TITRATION  FOR  F 
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FIGURE 
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FIGURE  17  INITIAL  F ATOM  CONCENTRATION  FROM  HCI  TITRATION 


?1 


amission  Spectroscopy  Measurements 


The  emission  spectroscopy  instrumentation  is  discussed  extensively  in 
References  1,  3 and.  the  addendum  to  Reference  1.  The  same  principles  apply  to 
the  fundamental  spectra  obtained  here  and  further  discussion  is  contained  in 
Reference  8.  Note  that  in  general,  it  was  planned  to  compare  the  experimentally 
determined  Einstein  A coefficients  for  overtone  transitions  determined  experi- 
mentally by  Siieo  in  his  Ph.D.  thesis  at  Cornell,  January  1975»  and  published 
in  Reference  9 with  those  calculated  by  Herbelin  and  Emanuel  (Ref.  10 ) by 
comparison  of  Av  = 1 and  Av  = 3 results  for  v = 3>  4 in  the  present  experiment. 
In  general,  Herbelin  and  Emanuel  values  gave  excellent  agreement,  while 
values  quoted  by  Siieo  and  Cool  (Ref.  5)  are  closer  to  Herbelin  and  Emanuel's 
values  than  they  had  been  at  the  time  of  the  thesis  work. 

Vibrational- rotational  level  populations  obtained  from  the  data  are 
presented  in  Section  IV. 

Contaminants 


Possible  contamination  of  the  gas  mixture  has  been  considered  from  the 
standpoint  of  total  leak  rate  of  the  experiment  and  purity  of  gases.  The 
system  could  be  pumped  to  less  than  0.001  torr  blank-off  pressure  through  an 
in-line  Hi?  trap  to  reduce  condensable  pressure  even  further. 

With  respect  to  hydrogen  bearing  contaminants,  e.g.,  HgO  and  CxHy  types 
it  is  possible  to  use  the  probe  laser  to  estimate  total  H contamination  by 
noting  the  HF  absorption  level  that  occurs  when  is  added  to  the  mass  heater 
but  H2  or  HC1  are  not  present.  Any  HE  formed  indicates  some  of  the  F atoms 
have  reacted  with  a hydrogen  bearing  contaminant.  Generally  such  spurious 
HF  formation  is  below  0.0003  torr  and  was  checked  for  all  tests.  Hence,  such 
contamination  generally  is  less  than  15  percent  of  the  lowest  values  of  initial 
F used.  Total  leak  rate  into  the  experiment  was  at  or  below  0.0004  torr  for 
operating  flow  and  pressure,  and  since  this  leak  is  primarily  N2  the  spurious 


8.  Lynds,  L. , "Investigation  of  Master  Oscillator  Power  Amplifier  Systems 
for  Chemical  Lasers,"  Final  Report,  Contract  F29601-73-C-C072 , UTRC 
Report  IT9 11653-24,  East  Hartford,  CT,  August  1974. 

9.  Siieo,  R.  N.  and  T.  A.  Cool,  "Overtone  Emission  Spectroscopy  of  HF  and 

DF:  Vibrational  Matrix  Elements  and  Dipole  Moment  Function,"  J.  Chem. 

P4ys.,  62,  117,  1976. 

10.  Herbelin,  J.  M.  and  G.  Emanuel,  "Einstein  Coefficients  for  Diatomic 
Molecules,"  J.  Chem.  Phys . , 60 , 689 > 1974. 
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level  of  EF  is  likely  related  to  contaminants  in  the  argon  for  these  highly 
dilute  experiments.  The  argon  used  was  obtained  from  standard  UTRC  sources 
and  the  purity  is  about  equivalent  to  a high  purity  or  prepurified  grade  so 
the  total  contaminant  level  of  the  argon  at  20  torr  would  be  between  0 . OOGh 
and  0.CC10  torr,  with  some  fraction  EoO  or  CH^ . The  hot  ?2  probe  laser  contam- 
inant test  can  identify  cylinders  of  argon  which  obviously  have  hydrogen 
bearing  contamination  levels  much  higher  than  typical  values.  On  one  such 

occasion  spurious  EF  signals  were  of  the  order  of  the  Fp  pressure  and  were 

related  to  the  argon  by  systematically  turning  off  sections  of  the  multi- 
cylinder argon  tank  farm  that  supplies  the  experiment;  in  this  case  two  of 
the  ten  LA  cylinders  in  the  farm  were  found  which  when  shut  off  caused  the 
high  spurious  EF  signal  to  go  away  and  these  were  discarded. 

In  summary,  contaminants  0E,  CE  and  spurious  EF  are  0.0003  torr  for 
this  experiment  while  is  probably  s O.OOOU  torr.  (Smaller  amounts  of  0P 
are  likely  present.)  Due  to  the  fast  nature  of  E?:E?  and  E?:E  interactions, 
the  fast  formation  of  initial  EF'  levels  of  the  order  of  F .002  torr)  it  is 

not  likely  that  these  contaminants  levels  will  have  a significant  effect  or. 

interpretation  of  present  experimental  results. 

Discussion  and  Recommendations 

The  mode  of  data  acquisition  here  has  required  that  the  conditions  for 
each  case  be  repeated  at  least  three  times  in  order  to  acquire  the  three  types 
of  experimental  information  that  are  basic  to  determination  of  the  production 
of  EF  and  the  individual  vibrational  level  populations  up  to  v =6;  (l)  probe 
laser  absorption,  (2)  absolute  IE  emission  spectroscopy,  (3)  absolute  overtone 
emission  spectroscopy  for  Lv  = 3*  In  the  present  program  the  application  of 
EC1  titration  to  determination  of  ? atom  level  was  evaluated  and  improvements 
made  during  the  course  of  the  program  such  that  the  final  determination  of  F 
atom  partial  pressure  with  metered  EC1  was  made  by  repeating  each  case  a fourth 
time.  Such  a mode  of  data  acquisition  requires  that  great  care  be  taken  in  the 
calibration  and  application  of  all  instrumentation  used  to  define  conditions 
and  that  cross  checks  be  made  periodically  among  the  types  of  measurements. 

Such  care  and  cross  checks  were  made  to  the  extent  possible  with  a manual  method 
of  data  acquisition.  For  example,  all  flowmeters  were  calibrated  with  known 
volumes,  accurate  pressure  gauges  and  a stopwatch  and  were  checked  periodically 
during  the  course  of  data  acquisition.  Calibration  of  the  emission  spectrometer 
with  a standard  source  was  possible  at  any  location  of  the  spectrometer  system 
along  the  flow  tube  and  such  calibration  checks  were  performed  before  and  after 
each  set  of  data.  Also,  when  one  type  of  measurement  was  being  made,  say  IE 
emission,  then  the  probe  laser  was  placed  at  a downstream  location  and  the  1 - 1 
of  absorption  noted.  The  velocity  probe  was  placed  at  a downstream  location  and 
th<  ;ent  r line  pitot  pressure  noted  for  each  run.  In  addition  where  possible 


within  the  scope  of  this  program  sets  of  data  were  repeated,  principally  here 
the  probe  laser  absorption  for  cases  1 through  4 of  Table  2.  Some  emission 
spectra  for  cases  1 through  4 were  repeated  during  the  course  of  HF  additive 
tests  since  the  zero  HF  additive  condition  was  referenced  again.  Results  gave 
reasonable  agreement  with  prior  data.  Repeat  data  for  the  H atom  cases,  5 and 
6,  were  not  acquired. 

Nevertheless,  the  above  procedures  are  difficult  to  follow  in  order  to 
arrive  at  a set  of  experimental  results  with  the  degree  of  internal  consistency 
required  to  evaluate  both  production  and  vibrational  level  population  at  various 
stages,  of  the  development  of  the  chain  along  the  flow  tube  in  analytical  terms 
as  related  to  the  conditions  of  the  experiment.  Such  a procedure  was  followed 
here  because  it  was  justified  relative  to  the  level  of  analytical  understanding 
of  K2/F2  chain  reaction  kinetics  at  the  time  and  was  consistent  with  the  experi- 
mental stage  of  the  art  for  the  thermally  initiated  flow  tube. 

In  terms  of  the  methods  employed  to  acquire  the  data,  the  experimental 
state-of-the-art  and  the  scope  of  the  program,  the  results  obtained  here  and 
from  data  analysis,  Section  IV,  are  very  encouraging.  Of  particular  signifi- 
cance for  absolute  measurement  interpretations  are  the  good  calibration  check 
of  the  probe  laser  with  known  HF  addition  (Fig.  14),  the  near  1:1  correspondence 
between  HF  formed  and  HC1  added  during  the  metered  titration  tests  (Figs.  15  and 
16),  and  the  agreement  between  ultimate  HF  levels  predicted  from  measured  con- 
ditions with  measured  ultimate  levels  of  HF  using  the  probe  laser  (Figs.  43  and 
4-7  for  case  l).  Absolute  emission  measurements  have  been  checked  several  times 
and  the  agreement  between  Av  = 1 and  Av  = 3 for  v = 3,  ^ is  encouraging  since 
the  present  understanding  of  transition  probabilities  (Refs.  9 and  10)  seems  to 
have  narrowed  the  discrepancies  a great  deal.  For  the  present  cases  a large 
fraction  of  the  total  production  of  HF  early  in  the  flow  tube  close  to  the  H2 
injector  is  in  vibrationally  excited  levels  v = 1,  2;  hence,  the  accuracy  of 
emission  data  is  important  in  determining  total  production  as  well  as  the 
accuracy  of  probe  laser  data.  At  the  end  of  the  reaction  zone  ground  state  HF 
is  dominant  and  the  probe  laser  is  of  greatest  significance.  In  between,  both 
measurements  are  relevant.  Further  work  could  be  done  to  scale  cases  such  as 
3 and  4 with  pressure  and  velocity  to  obtain  chemical  rate  information  based  more 
extensively  on  probe  laser  data  only.  Probe  laser  studies  with  the  HF  additive 
at  one  downstream  station  showed  a simple  sum  relation  and  such  measurements  at 
earlier  stations  should  show  enhanced  reliance  on  probe  laser  data  for  HF  pro- 
duction because  of  enhanced  deactivation  to  the  ground  state.  Data  at  various 
levels  of  initial  F should  be  extended  as  a simple  check  on  HF  production  results. 
In  this  regard,  case  6 results  need  to  be  repeated  and  an  experimental  check  made 
on  the  ultimate  HF  level. 
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Present  results  give  a good  indication  of  the  detailed  type  of  physical 
behavior  of  HP: HP  interaction  that  is  occurring  and  do  indicate  a chemical  rate 
of  HF  formation.  In  light  of  the  present  state  of  experimental  results  for 
cases  5 and  6 (Table  3 ) , however,  and  the  need  for  some  revisions  to  HP: HP 


modeling,  the  identification  of  an  overall  role  for  H atoms  in  the  chain  is 
incomplete  at  this  time. 


For  the  future,  further  analytical  work  is  needed  to  compare  with  present 
experimental  results.  When  defining  future  experiments  of  this  type  certain 
changes  are  more  or  less  obvious,  particularly  the  HP  level  very  close  to  the 
Eq  injector  must  be  examined  to  document  the  level  of  HF  that  may  be  present 
early  in  the  chain  due  to  small  recirculation  patterns  in  the  wake  of  the 
injector  bars.  In  this  regard  slight  modification  of  the  Eo  injector  should 
be  possible  to  purge  such  recirculation  zones.  It  would  also  be  desirable  to 
speed  up  the  data  acquisition,  take  all  types  of  measurements  during  the  same 
run  and  to  automate  the  data  acquisition.  In  addition,  scaling  of  results 
with  diluent  pressure  would  be  desirable.  However,  beyond  this  it  is  likely 
that  further  progress  will  be  linked  to  the  ability  to  provide  better  knowledge 
and  control  of  overall  conditions  on  close  to  an  instantaneous  basis;  progress 
of  this  type  as  well  as  automated  data  acquisition  could  best  be  done  by 
linking  the  experiment  to  a small  computer  in  the  laboratory. 

Summary  and  Conclusions 

A CW  flow  tube  H2/F2  chain  reaction  experiment  was  conducted  with  thermal 
initiation  to  measure  population  distribution  in  the  various  Y and  J levels  and 
the  ground  state  of  KF  over  a range  of  conditions  for  which  results  likely  coal 
provide  rate  information  for  the  overall  reaction  and  molecular  energy  transfer 
processes  that  were  uncertain. 


Hot,  partially  dissociated  F0  from  a mass  heater  is  adjusted  to  room 
temperature  by  nixing  with  argon  rapidly  enough  to  preserve  F.  Additional 
cold  F-  is  added  to  adjust  overall  ?2  dissociation.  Hydrogen  and  additional 
diluent  are  injected  through  a fine  scale  mixer  and  the  chain  reaction  proceeds 
under  dilute,  near  isothermal  conditions  typically  09.7  percent  argon  with 
IT  < 2C°C.  The  flow  tube  diameter  is  D = 5 cm,  ? = 1C  torr  and  velocity 
v ~ Vj*  ;n,  •.  Relatively  uniform  :onditlons  pr  vail  over  the  reaction  length 

bserved  L = 50  cm:  diffusion  of  HF  and  K to  the  walls  are  estimated  to  b°  abou 
percent  and  25  percent  effects  by  x = pC  cm  and  HF(v)  and  H loss  to  the  walls 
or  inhibited  further  by  coating  the  walls  with  Teflon.  The  reaction  was  obs°r 
at  ;even  viewports  along  the  reaction  length:  each  viewport  was  a narrow  slot  i 


tube  coupled  t; 
was  added  to  -s 


: port  s al  : with  a sapphire  window.  Met<  r 1 pui 


t j . 


Initial  F pressure  was  measured  with  metered  HC1  titration  and  HF  probe 
laser  absorption.  Ground  state  HF  was  measured  with  probe  laser  absorption. 
Vibrational  levels  in  H?  v = 1,  2,  3,  ^ were  measured  with  absolute  I?., 

Av  = 1,  emission  spectroscopy.  Vibrational  levels  v = 3>  5,  6 were  measured 

with  overtone  emission  spectroscopy  using  Av  = 3*  Temperature  profiles  were 
measured  with  a thermocouple  probe.  Velocity  profiles  were  measured  with  a 
pitot  probe  and  results  incorporated  in  the  analysis.  Partial  pressures  of 
reactants  ?2  ^2  were  calculated  from  measurements  of  total  pressure  and 

?2  a-nd  H2,  and  argon  flowrates  using  calibrated  mass  flow  transducers  and  roto- 
meters.  Pressure  distribution  along  the  flow  tube  was  measured  with  a trans- 
ducer and  wafer  scanner. 

The  objectives  of  the  flow  tube  experiment  were  sought  in  two  series  of 
tests:  (l)  operation  over  a wide  range  of  F2/H2  with  fixed  F and  HF  addition 

to  obtain  reaction  rate  and  HF:HF  interaction  information  and  (2)  setting  of 
higher  F maintaining  constant  and  nearly  constant  HF  production  to  enhance 
the  effect  of  H:HF  collisions. 

The  desired  conditions  and  data  accuracy  were  obtained  to  a large  extent 
except  within  the  remaining  uncertainty  in  modeling  of  KF:HF  interaction  the 
HF  production  was  not  close  enough  for  the  two  cases  selected  in  (2)  to  draw 
any  definite  conclusions  about  H:HF  interaction  at  this  time. 

Results  as  analyzed  give  a good  indication  of  the  detailed  type  of 
physical  behavior  of  HF:KF  interaction  that  is  occurring  on  both  a short  and 
long  time  scale  in  the  development  of  the  chain  and  allow  a conclusion  to  be 
reached  about  the  type  of  modeling  that  must  be  employed.  For  chemical  rates 
results  give  a strong  indication  of  a lower  limit  to  the  hot  reaction  rate  of 
about  1/2  the  present  value  at  room  temperature . Further  investigation  to 
clear  up  a small  uncertainty  about  initial  HF  levels  very  close  to  the  injector 
and  some  extension  of  the  conditions  investigated  could  narrow  the  chemical  rate 
uncertainty  even  further. 
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SECTION  IIx.  FJL5ZD  REACTION  EXPERIMENTS 


Fu.sed  reaction  experiments  were  used  to  study  the  chain  reaction  in 
premixed  H and  F . The  reaction  is  initiated  by  producing  fluorine  atoms 
in  the  mixture  by~dissociating  F,„  with  pulsed  IT/  radiation  from  a doubled 
frequency  ruby  laser.  Three  types  of  experiments  were  performed:  First, 

a study  was  made  of  the  temporal  development  of  populations  in  vibrational 
states  of  HF  produced  from  the  reaction.  Fluorescence  measurements  were 
used  to  track  excited  state  populations  and  absorption  of  cw  radiation  to 
follow  ground  state  formation.  Second,  small  signal  gain  was  determined  for 
v(i-O)  and  v(2-l)  transitions  by  use  of  a gain  probe.  Third,  pulsed  laser 
radiation  from  the  reacting  mixture  was  studied  using  a laser  configuration. 

The  rationale  for  these  experiments  is  that  they  would  serve  as  an 
experimental  check  on  the  laser  computer  program  as  modified  by  results  from 
flow-tube  experiments  and  other  available  data.  Pulsed  experiments  were 
chosen  as  being  relatively  simple  in  that  the  problems  attendant  with  mixing 
and  flow  are  not  present,  the  reaction  is  initiated  over  the  total  vol  one 
in  a short  time  period  and  changes  in  temperature  and  pressure  could  be 
monitored  in  a constant  volume  cell.  The  experimental  results  have  been 
compared  with  computer  calculations  of  population  development  as  the 
kinetics  were  tested  with  variations . in  initial  F atom  concentration  and  in 
the  H /F  ratio.  These  comparisons  are  presented  in  Section  IV  of  this 
report. 


Pulsed  Fluorescence  Experiment 


The  chain  reaction  between  fluorine  and  hydrogen  was  initiated  by 
dissociating  a small  fraction  of  the  fluorine  with  radiation  from  a frequency 
doubled  ruby  laser.  Diagrams  of  the  experiment  are  shown  in  Figures  lo 
and  19.  Fluorescence  from  the  chemically  formed  and  vibrationaiiy  excited  HF 
was  detected  using  appropriate  band  filters  with  an  AuSe  detector  for 
v(l-0.,  an  Ir.Sb  detector  for  v(2-0)  and  with  a Varian  L5E  photomultiplier 
for  v ' 3-0 ) , v(-.-_  , v(5-2)  and  v(6-3).  Intensity  changes  were  measured  to 
determine  population  changes  in  the  various  levels  with  time.  Ground  state 
HF  v=0 ) was  determined  by  absorption  measurements  using  a cw  HF  probe  _aser. 
Data  were  recorded  on  scope  pictures.  Rise  times  for  the  detection  systems 
were  better  than  0.5  usee.  Pressure  and  temperature  increases  during  the 
reaction  were  measured  isii  ' a Kistler  pressure  transduc  r in  a : nstan* 


detector 


FIGURE  18  FLUORESCENCE  EXPERIMENT 


volume  cell  momentarily  closed  ot’f  with  solenoid  valves.  The  cell,  made  of 
Kel-F,  was  15-2  cm  long,  1.2  cm  in  diameter,  and  was  fitted  with  sapphire 
windows.  The  Kistler  gauge  was  calibrated  against  a capacitance  nanometer 
using  pressure  pulses  over  the  pressure  range  of  interest.  A sample  trace 
of  the  gauge  outputs  is  shown  in  Figure  20. 

Fluorine  issociation  was  accomplished  using  a ruby  laser  that  delivered 
up  to  0.6  Joules  of  pulsed  radiation  at  69b2$  within  50  to  70  nsec.  The 
frequency  of  this  radiation  was  doubled  using  a temperature  tuned  RDA  crystal 
to  provide  up  to  0.l3  Joules  of  energy  at  3b71n.  Pulse  energy  was  measured 
with  a Scientech  energy  meter  calibrated  with  a biackboay  using  a shutter. 

The  pulsed  beam  had  a diameter  of  10  to  11  mm  which  closely  matches  the  cell 
IS  (12  mm)  for  excitation  of  the  whole  gas  volume.  The  concentration  of  F 
atoms  produced  was  determined  by  carefully  measuring  the  energy  (with  correc- 
tions for  window  losses)  and  calculating  dissociation  using  the  absorption 
coefficient  for  F„  determined  by  Steunenberg  and  Vogel  (Ref.  11 ).  The  value 

O L-  " *4  O . 

for  a at  3-"lA  is  1.21  x 10  at  273  K;  corrections  to  this  value  were  made 
for  the  actual  experimental  temperature.  Attempts  to  determine  F atom  content 
by  reaction  with  RC1  and  measurement  of  HF  produced  from  the  reaction  were 
unsuccessful.  In  a Kel-F  cell  the  degree  of  reaction  with  KC1  was  found  to 
be  variable.  In  an  aluminum  cell  the  reaction  did  not  go.  Although  this 
method  was  successfully  demonstrated  as  described  in  Section  V,  it  appears 
to  be  not  suited  to  small  volume  - static  flow  experiments  and  wall  inter- 
actions apparently  have  a significant  effect. 

Gas  mixtures  were  made  by  flowing  H , F_ , and  He  or  Ar  first  through  a 
mixing  section  consisting  of  three  meters  of  aluminum  tubing  in  the  form  of 
a coil,  next  through  the  fluorescence  cell  fitted  with  Kel-F  solenoid  valves 
with  a by-pass  valve  and  finally  through  a calibrated  orifice  to  the  vacuum 
pumps.  The  fluorine  was  passed  through  NaF  to  remove  HF.  Analysis  of  the 
F,  by  gas  chromatography  showed  0^  content  of  less  than  0.1  percent.  Gas 
pressure  was  measured  with  a capacitance  manometer  which  had  been  calibrated 
with  a mercury  manometer.  The  best  defined  fluorescence  data  were  obtained 
using  premixed  F^  in  argon  and  H„  in  argon.  These  mixtures  were  made  by 
adding  argon  to  lb. 7 psi  F_  for  a total  pressure  of  300  psi  and  adding  argon 
to  lb. 7 psi  H _ for  a tctal"pressure  of  500  psi. 

In  a flowing  system  where  cr.e  gas  is  added  after  another,  pressure 
measurements  do  not  give  a true  measure  of  '•oncentration  because  of  a change 
in  throat  area  of  the  exit  aperture  as  th»  boundary  layer  changes  with 
increasing  pressure.  The  flow  system  was  'alibrated  for  this  effect  by 

11.  Steunenberg,  R.  K. , and  R.  C.  Vogel,  "The  Absorption  Spectrum  of  Fluorine," 

J . Am . Chem . Soc.,  73,  901,  1956. 
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FIGURE  20  CALIBRATION-KISTLER  GAUGE  VS  CAPACITANCE  MANOMETER 
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introducing  argon  through  two  valves  in  sequence,  turning  the  valves  off  in 
the  opposite  sequence  while  reading  pressures,  ihe  dependence  of  flow  rate 
on  pressure  was  obtained  in  the  following  way: 

Let  p = test  cell  static  pressure 

A = exit  valve  effective  throat  area 
f = molar  flow  rate 

Now  for  choked  flow  through  the  valve 

f/  = Cp 
'A 

where  C is  a proportionality  factor  depending  on  the  properties  of  the  gas. 
If  the  gas  properties  remain  essentially  the  same,  this  factor  will  remain 
constant.  Similarly,  the  effective  throat  area  (dependent  on  boundary  layer 
thickness)  will  then  depend  only  on  the  pressure.  To  find  A = A(p)  consider 
the  following: 

Add  a gas  mixture  (say  Ar  + Fg)  to  a known  pressure  p.  Then  we  have 

f = CpA 

Suppose  another  valve  is  opened  to  allow  another  mixture  (Ar  + H , 
raising  the  pressure  to  p' . Then  ^ 

f* » _ r'-r' t r t 

■*»  - ^ y a 

Now  f and  f are  not  known  absolutely  but  the  ratio  f/f'  may  be  found. 
Suppose  f^  is  the  flow  rate  of  the  first  gas  added  to  the  cell  and  f that 
of  the  second . Then  ^ 

f = f 
X1 


Also  let  p,  and  p^  be  the  partial  pressures  of  the  two  gases  in  tne  final 
mixture . Now 
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Thus  if  the  partial  pressure  of  the  first  gas  in  the  final  mixture  is 
known,  the  change  in  effective  area  is  known.  A series  of  tests  over 
various  pressures  were  run  tc  obtain  A = A(p)  and  the  values  were  used  tc 
construct  the  correction  curve  shown  in  Figure  21. 
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pressure  of 


:an  be  found  by  closing  off  the  test 


and  initiating  the  reaction  of  F_  + H..  The  resultant  KF  formed  can  then 
be  related  to  p-  or  p..  , whichever  is  smaller.  Knowledge  of  the  relative 
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Absorption  by  HF  of  radiation  from  a small  cv  probe  laser  was  used  to  check 
the  actual  Fp  and  Hp  concentrations  by  determining  the  amount  of  HF  produced  in 
Fp  rich  and  in  Hr  rich  mixtures . This  radiation  was  also  used  in  the  pulsed 
xpe rime nos  for  following  the  kinetics  of  ground  state  HF  production.  A small 
cw  HF  r.rtfce  laser  was  developed  for  these  purposes  which  has  excellent  ampli- 
tude and  frequency  stability  (Ref.  12).  The  laser  radiation  is  single  mode 
and  can  be  tuned  over  the  ioppler  broadened  power  curve  of  about  300  MHz.  The 
capability  of  tuning  to  line  center  is  necessary  because  absorption  measure- 
ments made  at  frequencies  removed  from  line  center  will  significantly  err 
,?ef.  This  was  especially  so  in  mixtures  containing  100  to rr  of  argon 

vhich  caused  pressure  broadening  and  line  shift  in  the  HF  absorption  curve, 
in  order  to  evaluate  this  effect,  absorption  measurements  were  made  for  pure 
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FIGURE  21  VARIATION  OF  EXIT  THROAT  AREA  WITH  CHAMBER  PRESSURE 


HF  and  HF  diluted  with  100  torr  Ar  as  the  laser  frequency  varied  over  its 

linewidth.  A sample  measurement  of  I and  for  0.33  torr  HF  is  shown  in 

Figure  22.  Similar  measurements  for  HF  at  0.13  torr  in  Figure  23  snow 

that  for  Dure  KF  a plot  of  in  I/I  vs.  frequency  Droduces  an  abscroticn 
* " p 

curve  centered  at  the  laser  v , wnereas  for  the  sample  containing  100  torr 

o ^ 

of  argon  there  is  a freauency  shift  65  MHz  from  v and  considerab-e  ores- 

o . 

sure  oroadenmg.  These  results  are  m agreement  with  previous  measurements 
of  line  shifts  and  broadening  (Ref.  13).  In  the  latter  case,  measurement  of 
absorption  at  the  laser  line  center  results  in  a value  over  50  percent  lower 
than  the  value  for  pure  HF.  By  making  many  such  measurements  for  laser  lines 
1PU  and  1?6  the  plot  shown  in  Figure  2h  was  constructed  and  used  to  interpret 
absorption  measurements.  Using  these  corrections  to  the  absorption  measure- 
ments, the  KF  produced  by  a cw  arc  lamp  from  K?  rich  and  F^  rich  mixtures 
with  100  torr  Ar  was  determined.  In  an  aluminum  ceil  only  50  to  60  percent 
of  the  expected  HF  was  found,  indicating  that  the  aluminum  walls  inhibited 
the  reaction.  Up  to  95  percent  of  the  expected  HF  was  found  in  a Hel-F 
cell,  and  this  cell  was  chosen  for  the  pulsed  fluorescence  studies. 


Pulsed  Fluorescence  Data 


The  data  for  the  fluorescence  experiment  were  obtained  as  scope  traces 
of  outputs  from  the  detectors  for  pressure  (Kistler  gauge),  fluorescence 
from  KF  levels  v = 1,  v = 2,  and  v z 3,  and  for  intensity  of  transmitted  cw 
laser  radiation  (for  v = 0).  An  example  of  these  traces,  which  were  recoried 
simultaneously,  is  shown  in  Figure  25  for  a gas  mix  containing  0.2  tcrr  H„ , 
1.0  torr  F_  and  100  torr  Ar.  The  formation  of  KF  (v  = 0)  is  evidenced  as'a 
decrease  in  intensity  of  transmitted  probe  laser  radiation  which  occurs  over 
a relatively  long  time  of  milliseconds.  Immediately  after  the  initiating 
pulse  a small  gain  can  be  seen  in  this  trace.  Also  the  formation  and 
deactivation  of  HF  (v  = l)  and  the  change  in  pressure  due  to  temperature 
increase  from  the  reaction  are  on  this  long  time  scale.  However,  the  forma- 
tion and  decay  of  KF  (v  = 2)  and  HF  (v  > 3)  occurs  in  a much  shorter  time. 

In  Figure  2b,  fluorescence  traces  for  the  individual  vibrational  transitions 
v(6-3),  v(5-2),  v(l-l)  are  shown.  These  pictures  which  were  taken  in  three 
individual  pulse  experiments  using  the  same  pulse  energy  demonstrate  the 
extremely  fast  formation  of  HF  in  levels  above  , v = 2)  ar.d  the  very  fast 
relaxation  of  these  levels. 


Attempts  were  made  to  measure  the 
order  to  determine  absolute  populations 
needed  to  be  placed  right  against  the  ce* 
was  not  possible  to  iefine  the  source  dis 

13.  Jaffe,  , et  ai . , "Pressure  Induced 
Due  to  Noble  daces,"  J . P.  — . . Phys . , 


absolute  intensity  of  fluorescence  in 
in  each  level.  Because  the  ie* -.-etors 
1 windows  for  adequate  signal,  it 
tance  and  aperture  needed  for 


FIGURE  22  ABSORPTION  OF  1P4  RADIATION  BY  HF  AND  HF  IN  ARGON 
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FIGUHt  26  FLOURESCENCE TRACES 
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calibration.  Also  attempts  to  measure  fluorescence  through  a 1/4  meter 
spectrometer  were  not  successful  because  of  inadequate  fluorescence  strength. 

Fluorescence  experiments  were  run  with  four  kinds  of  mixtures:  ( l) 

mixtures  of  E?  and  ?2 , (2)  mixtures  of  H2  and  F2  vith  He  added,  (3)  mixtures 
of  Hq  and  Fg  with  argon  added,  and  (4)  mixtures  of  premixed  Ho  and  A r with 
premixed  ?2  Ar.  For  the  experimental  cases  with  the  premixed  gases,  the 
gas  mixture  composition  was  accurately  known  and  the  dilution  was  sufficient 
to  keep  the  temperature  increase  within  moderate  limits  of  4 0 to  90°C. 

These  two  criteria  were  not  met  for  the  first  three  mixtures.  For  this 
reason  the  data  for  the  premixed  gases  will  be  presented  here  and  are  the 
ones  compared  with  the  computer  model.  The  other  mixtures  will  be 
discussed  below. 

Variations  were  made  in  Kp  and  F content  using  mixtures  containing  (A) 

C.2  torr  Ho,  1.0  torr  F2 , 98.8  torr  Ar;  (E)  2.0  torr  Hg , 1.0  torr  Fg,  97 

torr  Ar;  (c)  0.5  torr  H2 , 1.0  torr  ?2,  98.5  torr  Ar;  (D)  0.5  torr  H p,  0.5 

torr  F2,  99  torr  Ar;  (S)  1.0  torr  Ho,  0.5  torr  F- , 98.5  torr  Ar.  For  each 

of  these  mixtures  variations  were  made  in  the  initial  F atom  concentrations 

by  modifying  the  exciting  pulse  energy. 

An  example  for  each  of  th-  se  mixtures  of  the  KF  (v  = 0)  number  density 
and  gas  temperature  time  dependence  is  shown  in  Figures  27,  28,  29,  30,  31 
and  32.  These  graphs  derive  from  scope  traces  such  as  Figure  25  where 
temperature  is  calculated  from  measured  pressure  rise,  and  absorption  of 
the  cw  beam  is  used  to  calculate  number  density  with  corrections  for  line 
broadening  and  temperature  dependence  of  the  absorption  coefficient.  Generally, 
it  is  found  that  the  temperature  peaks  at  or  before  the  maximum  KF  production. 
For  mixtures  A,  3 and  C,  the  HF  produced  reasonably  approximates  that  of  the 
maximum  amount  possible  through  complete  reaction.  For  cases  D,  2 and  F, 
the  HF  amount  is  somewhat  less.  It  was  found  that  for  each  mixture  the 
amount  of  KF  produced  was  dependent  on  initial  F atom  content  (pulse  energy) 
and  tended  to  produce  the  theoretical  maximum  at  high  F atom  content.  This 
effect  is  discussed  in  the  Data  Analysis  section  and  is  apparent  in  Figures 
57  through  62.  A complete  listing  is  contained  in  Table  5 of  all  the  experi- 
ments for  which  data  were  obtained  using  the  mixtures  A,  B,  C,  D,  £ and  F. 

These  data  are  further  discussed  in  the  Data  Analysis  section. 

Table  6 contains  a summary  of  the  experiments  that  wen  run  for  the 
other  kinds  of  mixtures  excluding  premixed  gas  mixtures.  The  first  experi- 
ments (Table  6-A)  were  with  undiluted  R~  and  Fc . Examples  of  fluorescence- 
traces  obtained  with  th°se  mixtures  are  shown  in  Figures  33  and  31*.  Many 
of  these  runs  were  analyzed  to  get  rate  coefficients:  however,  beca use  of  the 
extreme  temperature  rise  of  up  to  10C  ' , we  have  littl  confidence  in  ohes- 

analyses  and  they  are  not  presented  he  . Cue  curious  aspect  of  thes-  xperi- 
ments  is  associated  witn  prereaction  which  was  •oinmonly  observ-i.  In  'as-  s 
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0.105J  H2  1.02  TORR  F 2 0.94  TORR  Ar  98.04  TORR 


FIGURE  27  HF  NUMBER  DENSITY  AND  TEMPERATURE  HISTORIES 
PULSED-FLUORESCENCE  CASE  A10 


0.1  34 J H2  2.22  TORR  F2  1.02  TORR  A.  96.70TOHH 


(Xo)l'3yniVd3dlAl31 


FIGURE  28  HF  NUMBER  DENSITY  AND  TEMPERATURE  HISTORIES 
PULSED-FLUORESCENCE  CASE  B4 


OHM 
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FIGURE  30  HF  NUMBER  DENSITY  AND  TEMPERATURE  HISTORIES 
PULSED-FLUORESCENCE  CASE  DH4 


• ' - 


Il2  0.24  TORR  1.01  TORN  Ai  98.75  TOHH 


FIGURE  32  HF  NUMBER  DENSITY  AND  TEMPERATURE  HISTORIES 
PULSED-FLUORESCENCE  CASE  FF2 
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TABLE  6.  SUMMARY  OF  PRELIMINARY  PULSE  EXPERIMENTS 


I 


A.  Mixtures  of  Ho  and  F?  (no  diluent ) 


Mix 

I 

1.0  Torr  Fg 

0.2 

Torr  H2 

Pulse 

Energy  (Joules): 

0.025 

0.04 

0.08 

0.035 

0.045 

0.08 

0.035 

0.05 

0.08 

0.04 

0.05 

0.09 

0.04 

0.C7 

0.097 

Mix 

II 

1.0  Torr  F2 

• 5 

Torr  H2 

Puls  e 

Energy  ( Joules ) : 

0.007 

0.024 

0.043 

0.065 

0.008 

0.024 

0.054 

0.067 

0.01C 

0.025 

0.057 

0.068 

0.011 

0.035 

0.053 

0.068 

0.012 

0.037 

0.060 

0.073 

0.012 

0.038 

0.06l 

0.035 

0.015 

0.040 

0.06l 

0.047 

0.019 

0.040 

0.063 

0.048 

0.02C 

0.041 

0.065 

0.059 

Mix 

III 

1.0  Torr  ?2 

0.6 

Torr  H2 

Pulse 

Energy  (Joules) 

1 J 

0.035 

0.06 

0.092 

0.035 

0.06 

Mix  r; 


1.0  Torr  ?2  1.0  Torr  H2 


Pulse  Energy  (Joules): 


C.035 

0.035 


fc>3 


.... 


0.056 

0.060 


0.074 

0.067 


TABLE  6.  SUMMARY  OF  PRELIMINARY  PULSE  EXPERIMENTS  (CONT'D) 

Mix  V 1.0  Torr  ?0  1.5  Torr 

Pulse  Energy  (Joules): 

0.028  0.$+3  0.057 

0.031  0.052  0.093 

Mix  VI  1.0  Torr  F2  2.0  Torr  H2 

Pulse  Energy  (Joules): 

0.031+ 

0.053 

0.091 

Mix  VII  2.0  Torr  F2  0.2  Torr  H2 

Pulse  Energy  (Joules): 

0.09 

Mix  VIII  8.0  Torr  F2  8.0  Torr  H2  80  Torr  N2 

Pulse  Energy  (Joules): 

0.063  0.C82 

0.063  0.085 

Mix  IX  0.5  Torr  ?2  5.0  Torr  Hq 

Pulse  Energy  (Joules): 

0.09 

. ~ : vj-.h  He linn 

1.0  Torr  F2  0.2  Torr  Hq  50  Torr  He 

Energy  (Joules): 

. . c 0.C7C 

. 0.081 

ol+ 


► 

TABLE  6.  SUMMARY  OF  PRELIMINARY  PUISE  EXPERIMENTS  (CONT'D) 

Mix  II  1.0  Torr  Fg  0.2  Torr  Hg  100  Torr  He 

Pulse  Energy  (Joules): 

0.C40  O.C45  0.C8  0.09 

0.044  0.05  0.08 

Mix  III  1.0  Torr  Fg  0.5  Torr  Hg  1.0  Torr  Ng 

Pulse  Energy  (Joules): 

0.027  0.033  0.04 

0.03  0.037 

Mix  IV  1.0  Torr  Fg  0.6  Torr  Hg  50  Torr  He 

Pulse  Energy  (Joules): 

0.05  0.065 

0.06  0.092 

Mix  V 1.0  Torr  Fg  1.0  Torr  Hg  50  Torr  He 

Pulse  Energy  (joules): 

0.028  0.056  O.G87  0.097 

0.035  0.060  0. 0914- 

Mix  VI  1.0  Torr  Fg  1.5  Torr  Hg  50  Torr  He 

Pulse  Energy  (joules): 

0.031  0.043  0.057  0.093 

0.034  0.051  0.075  0.097 

Mix  VII  1.0  Torr  Fg  2.0  Torr  Hg  50  Torr  He 

Pulse  Energy  (Joules): 

0.053  0.091 

0.088 
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TABLE  6.  SUMMARY  OF  PRELIMINARY  PUISE  EXPERIMENTS  (CONT'D) 


Mix  VIII  0.5  Torr  Fg  5»0  Torr  Hg  50  Torr  He 

Pulse  Energy  ( Joules ) : 

0.091 

0.093 

Mix  IX  8.0  Torr  Fg  8.0  Torr  Ho  80  Torr  Ng 

Pulse  Energy  (joules): 

0.025  0.063  0.082 

0.028  O.O65  0.085 


Determination  of  HF  produced  from  Fg-Hg-He  mixtures  to  determine  Fg 
Each  mixture  pulsed  3 times  for  complete  reaction. 


Mix  I 


1.0  Torr  Fo 


1.0  Torr  Ho 


50  Torr  He 


Test  1 - pulse  energy  (joules):  0.11 
Test  2 0.11 
Test  3 0.13 


0.085 

0.11 

0.10 


Mix  II 


>5  Torr  Fg 


.5  Torr  Hg 


50  Torr  He 


Test  1 - pulse  energy  (joules):  0.12 

Test  2 0.107 

Test  3 0.139 


0.i+5 

0.088 

0.142 


Mix  III  1.0  Torr  Fg 


1.5  Torr  Hg 


50  Torr  He 


Test  1 - pulse  energy  (joules):  0.105 
Test  2 0.123 
Test  3 0.095 


0.05 

0.1 

0.C8 


Mix  IV 


0.5  Torr  Fg 


1.0  Torr  H, 


50  Torr  He 


Test  1 - pulse  energy  (joules):  0.055 
Test  2 0.075 
Test  3 0.09 


0.085 

0.085 

0.09 


concentrations . 


0.11 

0.13 

0.11 


0.11 

0.11 

0.133 


0.09 

0.052 

0.098 


O.C95 

0.06 

0.1 


6c 


TABLE  6.  SUMMARY  OF  PRELIMINARY  FUISE  EXPERIMENTS  (CONT'D) 


v 


C.  Mixtures  with  Aracn 


Mix  I 100  Torr  Ar  1.0  Torr  Fg  0.2  Torr  S? 

Pulse  Energy  (Joules): 


0.C8 

0.0^7 

0.067 

0.098 

0.017 

0.05 

0.08 

0.10 

0.035 

0.05 

0.086 

0.13 

0.045 

0.055 

0.09 

0.17 

Check  of  Eg  content  by  determination  of  HE  using  multiple  pulses  with  above  mixture 


Test 

1 - Pulse  Energy 

0.06 

0.08 

0.06 

HE 

= 

C.29 

Torr 

Test 

2 

0.105 

0.11 

0.12 

HE 

= 

0.65 

Torr 

Test 

3 

0.075 

0.067 

0.06 

HE 

= 

0.38 

Torr 

Test 

4 

0.035 

0.113 

0.085 

HE 

= 

0.55 

Torr 

Test 

5 

0.12 

0.090 

O.O96 

HE 

= 

1.10 

Torr 

Test 

6 

0.05 

0.113 

0.095 

KF 

= 

0.59 

Torr 

Test 

7 

0.80 

0.65 

0.60 

HE 

= 

0.6l 

Terr 

Test 

8 

0.04 

0.08 

0.06 

HE 

= 

0.53 

Torr 

Test 

9 

0.098 

0.085 

0.12 

HF 

= 

0.95 

Torr 

Test 

10 

0.08 

o.c8 

0.065 

KF 

= 

0.55 

Torr 

Test 

11 

0.13 

0.13 

0.12 

HE 

= 

0.70 

Torr 

0.C8 

0.09 

0.11 

HE 

= 

0.64 

Torr 

Test 

12 

0.11 

0.10 

0.11 

HE 

C.41 

Torr 

0.115 

0.02 

0.08 

HF 

= 

C.53 

Torr 

Test 

13 

0.086 

0.11 

C.080 

HE 

= 

0.32 

Torr 

Test 

14 

0.085 

0.11 

0.11 

KF 

= 

0.22 

Torr 

Test 

15 

0.060 

0.080 

0.083 

HE 

= 

0.32 

Torr 

Test 

16 

0.17 

0.10 

HF 

= 

0.18 

Torr 

FIGURE  33  REACTION  OF  AN  H2-F2  MIXTURE 
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FIGURE  34  EFFECT  OF  INITIAL  F-ATOM  ON  H2-F2  REACTIONS 


where  prereaction  produced  in  excess  of  0.02  torr  of  HF  there  was  a slow  build 
up  of  HF  in  v = 1 and  above  immediately  after  the  UV  dissociating  pulse. 

After  some  delay  there  was  a sharp  increase  in  the  excited  population  which 

was  reflected  in  v = 0.  The  delay  time  was  shorter  for  increased  initial 

F atom  concentration.  These  effects  are  illustrated  in  the  scope  traces  of 
Figure  35.  It  would  be  interesting  to  look  for  gain  or  lasing  effects 

during  the  spiking.  We  surmise  the  spiking  behavior  is  due  to  inhibition 

of  the  chain  by  HF(v  = 0)  with  slow  buildup  of  temperature  until  an  explosive 
reaction  occurs. 

In  another  series  of  experiments  (Table  6-3)  He  was  added  to  the  H -F 
mixture  to  keep  the  temperature  down . Unfortunately  the  addition  of  helium 
disturbs  the  F.  concentration  because  as  3ott  has  pointed  out  (Ref.  I1*),  the 
viscosity,  average  molecular  weight  and  specific  heat  of  the  mixture  changes 
so  that  the  initial  concentrations  were  unknown.  An  attempt  was  made  to 
back  correct  the  data  for  these  effects.  However,  testing  the  mixture  by 
determining  HF  produced  through  absorption  measurement  gave  low  and  inconsis- 
tent results.  Similar  conclusions  were  reached  for  the  data  with  Ar  as 
diluent  (Table  6-C).  For  these  mixtures  there  was  uncertainty  about  the 
H concentration.  3ecause  of  the  high  temperature  or  uncertainty  about 
starting  concentrations,  data  from  the  experimental  runs  listed  in  Table  6 
were  not  analyzed  in  detail  and  were  not  used  to  evaluate  kinetics. 


Pulsed  Gain  Measurements 


Gain  measurements  were  made  for  pulsed  chain  reactions  in  mixtures  maae 
with  premixed  F./Ar  and  H^/Ar.  The  fluorescence  apparatus  illustrated  in 
Figures  13  and  V)  was  usee  except  the  cw  probe  beam  traversed  the  length  of 
the  cell  (15.2  cm)  at  a slight  angle  off  the  axis  of  the  doubled  ruby  beam. 
Intensity  of  the  cw  beam  was  monitored  with  a AuGe  detector  and  scope  pic- 
tures were  used  to  record  data.  Two  such  gain  traces  are  presented  in  Figure 
36  for  a IP 3 and  2P3  laser  line  with  peak  gains  found  at  1.3  and  l.d  percent/ 
cm  respectively.  A very  important  aspect  of  these  traces  is  that  the  gain 
typically  lasts  for  10  to  20  usee  after  the  icuble  ruby  pulse  while  the 
reaction  producing  HF  goes  on  for  milliseconds;  that  is,  most  of  the  reaction 
energy  does  not  result  in  laser  gain. 


For  a mixture  containing  0.2  torr  H^,  1.0  torr  F_,  and  98.8  torr  Ar, 
gain  measurements  were  made  for  the  lines  IP2,  1P3,  l?14,  I?6,  1?7,  2F3, 

2Pd,  2P5  and  2P6  with  variations  in  the  initial  F atom  concentration  (pulse 


11.  Bott,  J.  F. , "Gas  Dynamic  Corrections  Applied  to  Laser  Induced 
Fluorescence  Measurements  of  H?(v=l)  and  DF(V=l)  Deactivation," 
J.  Chem.  Phys.,  cl,  3^  Id , 19? d . 
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FIGURE  35  EFFECT  OF  INITIAL  F-ATOM  ON  REACTION  - 0.02  TORR  HF  PRESENT 


0.2TORR  H2,  I.OTORR  F2,  98.8  TORR  Ar 


1P3  GAIN  1.4%/CM 
PULSE  ENERGY  0.17  JOULES 


?p3  GAIN  1.3%. CM 
PULSE  ENERGY  0.19  JOULES 


FIGURE  36  GAIN  TRACES 


energy).  These  data  are  listed  in  Table  7-  In  Figure  37  the  gain  on 
various  lines  is  shown  for  experiments  using  0.12  joule  and  0.17  joule 
dissociation  pulses.  Over  the  range  of  lines  studied  the  gain  is  seen 
to  peak  at  1P4  and  at  2P3.  The  gain  on  any  particular  line  was  found  to 
monotonically  increase  with  increasing  pulse  energy  as  shown  in  Figure  38. 

For  the  mixtures  studied,  gain  at  constant  pulse  energy  of  the  1P4 
line  was  found  highest  for  0.5  torr  Ho,  1.0  torr  F0  and  98*5  torr  Ar.  The 
data  with  various  mixtures  are  listed  in  Table  7 and  plotted  in  Figure  39* 
An  analysis  of  these  data  and  comparison  with  the  computer  model  is  pre- 
sented in  Section  IV. 


Pulsed  Laser  Tests 

Laser  tests  were  made  using  the  mixtures  of  premixed  E^/Ar  and  p2/'Ar 
employed  for  fluorescence  and  gain  measurements.  Frequency  doubled  radia- 
tion from  the  ruby  laser  was  used  to  dissociate  the  fluorine  and  initiate 
the  chain  reaction.  The  laser  configuration  closely  followed  the  design  of 
Nichols , et  al  (Ref.  15). 


Previous  studies  of  pulsed  lasing  on  mixtures  include  the  work  of 

Basov  and  coworkers  (Ref.  l6)  who  demonstrated  pulsed  lasing  using  flash- 
lamp  initiation  of  the  chain  reaction  from  levels  up  to  v,  6-5,  in  1969-  In 
these  experiments  lasing  on  transitions  from  v = 4 showed  less  power  than  from 
other  levels  leading  to  the  conclusion  that  chain  branching  by  HF(v=4)  + 

Fg  -*  HP  + 2F  was  an  important  contributor.  These  experiments  also  suggested 
that  HF(v=l)  excitation  was  transferred  to  FU  because  the  1-0  transition 
was  not  observed.  Subsequent  workers  generally  did  not  observe  lasing  from 
transitions  higher  than  v = 4 -»  3-  Careful  studies  of  pulsed  lasing  were 
reported  by  Nichols  and  coworkers  (Ref.  15).  In  this  work  a frequency  doubled 
pulsed  ruby  laser  was  used  to  start  the  F?- Ho  reaction  and  laser-emission  was 
found  on  v = 3 "'2,  2 -*  1,  1 -»  0.  The  time  history  of  ten  lines  was  measured 
simultaneously  and  ground  state  KF  was  determined  with  a pulsed  probe  laser 
and  F atom  concentration  was  calculated  from  known  cross  sections  for  photo- 
dissociation. The  temporal  behavior  of  lasing  from  individual  lines  from  a 
pulsed  F, -Hg  laser  initiated  by  doubled  pulsed  ruby  has  been  reported  by 


15.  Nichols,  D.  B. , K.  H.  Wrolstad,  and  J.  Doyle  McClure,  "Time-Resolved 
Spectroscopy  of  a Pulsed  H--F-  Laser  with  Well-Defined  Initial  Con- 
ditions,” -J.  of  Appi.  Phys . , 5 5 , 12,  December  1974. 


Basov,  N.  0.,  L.  V.  Kulakov, 
and  ?.  N.  Lebedev,  "Fmission 
Hg  * Fg  Mixture,"  JETP  Lett. 


S.  ?.  Markin,  A.  I.  Mikitin, 
Spectrum  of  a Chemical  Laser 

£,  375,  1969. 


A.  N.  Uraevskii, 
Using  an 
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TABLE  7.  DATA  FROM  GAIN  MEASUREMENTS 


I.  Line  Variation  (Mixture  0.2H^  I.OF2  98.8Ar) 


Gain 

Pulse 

Gain 

Pulse 

Code 

Line  $/cm 

Enersy 

Code  Line 

$/  cm 

Energy 

144 

1P2  1.22 

• 153 

01  1P7 

• 39 

.16 

M3 

• 98 

• 137 

03 

•39 

• 15 

M2 

•90 

• 115 

02 

.11 

.12 

Ml 

.86 

.10 

N5  2P3 

1.3 

• 19 

LI 

IP3  1.35 

• 17 

N4 

1.2 

• 13 

12 

• 79 

• 15 

N3 

1-3 

• 19 

L3 

• 78 

• 13 

N2 

1.1 

.10 

N1 

1.0 

• 09 

H4 

1P4  I.49 

.165 

K2 

1.01 

.12 

K1  2P4 

1.09 

.16 

HI 

.76 

.10 

K2 

.66 

.12 

H3 

.82 

.10 

01  2P5 

.80 

• 19 

11 

1P6  .80 

.16 

02 

.48 

.11 

IS 

• 54 

.14 

• 

13 

• 55 

.12 

P3  2P6 

.80 

• 19 

PI 

.49 

.14 

J2 

1P7  -39 

.16 

P2 

.49 

.12 

J1 

J3 

.26 

.11 

II.  Mixture  Variation  (Line  LP4) 

Gain 

Pulse 

Gain 

Pulse 

Code 

Mixture 

$/cm 

Energy 

Code 

Mixture 

$/cm 

Energy 

H4 

O.2H2  1?2  98.8Ar 

1.49 

.165 

S2 

•5h2  if2  98. 

5Ar  2.32 

• 15 

H2 

1.01 

.12 

SI 

1.4o 

.11 

HI 

• 76 

.10 

H-3 

.82 

.10 

T3 

IH2  0.5F2  98 

•5Ar  .15 

• 133 

T1 

.14 

.127 

R1 

• 5H2  -5F2  99 At 

• 94 

• 15 

12 

.09 

.082 

R2 

• 54 

.11 

01 

2H2  1F2  97Ar 

no  gain,  pre- 

reac 

tion  complete 

absc 

.rption  of 

1?4 

lin*3 
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FIGURE  39  GAIN  1P4-VARI0US  MIXTURES 


Dolgov  (Ref.  IT).  The  laser  was  found  to  emit  on  transitions  v = 5 ■*  4, 

4 -»  3,  3 -*■  2,  2 -*■  1,  and  detailed  graphs  of  intensity  of  each  line  versus 
time  were  determined.  Pulsed  lasing  from  the  F^-H^  chain  reaction  was  also 
reported  by  Suc’nard  (Ref.  18)  using  flashlamp  initiation  and  with  time 
resolved  spectroscopy.  Lasing  from  transition  up  to  v = 6-5  was  observed 
using  a high  "Q"  optical  cavity  with  10  percent  output  coupling.  In  an 
earlier  study,  Suchard  (Ref.  19)  concluded  that  only  a.  7 percent  of  the 
total  laser  radiation  is  on  transitions  above  v(U-3). 


A schematic  of  the  laser  configuration  is  shown  in  Figure  40.  Pre- 
mixed gases  continuously  flow  through  a gain  ceil  made  from  Kel-F  with 
dimensions  of  length  6l  cm  and  diameter  12  cm.  The  laser  cavity  was  formed 
with  a fiat  total  reflector  and  a 10  meter  radius  92  percent  reflector 
separated  by  96  cm.  The  second  surface  of  the  partially  transmitting  mirror 
was  AR  coated  for  2.7  microns.  The  347l5  pulse  was  injected  into  the  laser 
cavity  by  reflection  from  a sapphire  fiat  coated  to  reflect  S polarised  34~i 
radiation  and  transmit  P polarized  2.7  micron  radiation. 


Two  sets  of  laser  tests  were  run.  For  the  first  set,  the  optical  com- 

i * o 

ponents  shown  in  Figure  **0  were  used.  Sapphire  windows  held  at  91  to  the 

laser  axis  were  attached  to  the  ends  of  the  gain  cell.  Initially,  sapphire 

windows  AP.  coated  with  MgO  were  used,  but  the  coating  quickly  degraded  and 

optical  losses  became  unacceptable.  Sapphire  windows  held  at  Brewster's 

angle  were  also  tried  in  an  attempt  to  cut  optical  losses.  However,  this 

arrangement  rejected  about  50  percent  of  the  doubled  ruby  radiation  even 

though  the  polarization  of  this  beam  was  oriented  for  optimum  passage  by  the 

window.  A germanium  flat  at  Brewster's  angle  was  pl'aced  between  the  cell  and 

0 

the  totally  reflecting  mirror  to  reflect  3^71A  radiation  from  tne  optical  tra 
and  prevent  feedback  to  the  RDA  crystal  with  possible  crystal  damage.  Radia- 
tion from  a cw  probe  laser  was  used  to  align  the  optical  components  and 
measure  individual  losses.  Values  found  were:  347-5  ref.ector  - 5 percent, 

gas  cell  - 24  percent,  germanium  flat  - lo  percent,  transmitting  mirror  - 
S percent.  The  optical  components  were  further  aligned  for  peak  laser 

output  by  flowing  a SF  ,-H  mixture  through  the  cell  and  by  repetitively 

o 2 
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FIGURE  40  LASER  EXPERIMENT 


pulsing  the  mixture  with  an  electrical  discharge  to  achieve  rapid  laser 
pulsing.  Two  pins  were  inserted  in  the  Kel-F  tube  as  electrodes. 

Measurements  were  made  of  relative  intensities  of  the  total  pulse  of 
HF  laser  radiation  and  of  each  spectral  component  using  gas  mixtures  con- 
taining (a)  0.2  torr  H , 1.0  torr  F?,  98.8  torr  Ar;  lb)  0.5  torr  K_,  1.0 
torr  F 5 , 98.5  torr  Ar;~(c)  1.0  torr  H , 0.5  torr  ?o,  98.5  torr  Ar;  (d)  0.5 
torr  K“,  0.5  torr  F„,  99  torr  Ar. 

o 

An  InSb  detector  monitored  the  total  pulse  and  a AuGe  (77  K)  detector 
was  used  at  the  output  of  a Jarrell  Ash  1/2  meter  spectrometer  for  measuring 
spectral  components.  The  detector  signals  were  displayed  and  photographed 
as  oscilloscope  traces.  Examples  of  such  traces  are  shown  in  Figure  41. 

Tata  for  the  first  set  of  laser  tests  are  listed  in  Table  8- A.  The  most 
energetic  lasing  at  2.5  microns  was  obtained  with  a gas  mixture  of  0.5  torr 
K, , 1.0  torr  F and  98.5  torr  Ar  which  is  the  same  composition  that  exhibited 
tne  largest  gam  (see  Table  7).  Lasing  was  not  observed  on  transitions 
higher  than  v(2-l)  in  any  tests  of  this  set.  It  is  not  surprising  that  highe 
transitions  were  not  observed  considering  that  for  a gain  length  of  6l  cm 
the  gain  for  individual  lines  varied  from  about  60  percent  per  pass  downward 
and  optical  losses  in  the  cavity  were  53  percent. 

For  the  second  set  of  laser  tests  the  germanium  flat  was  removed  and  the 

o 

90  sapphire  windows  were  replaced  with  3aF^  windows.  Optical  losses  with 
this  configuration  were:  3^7lX  reflector  -""5  percent,  gas  cell  - 13  percent, 

transmitting  mirror  - 8 percent  for  a total  of  26  percent.  A team  splitter 
diverted  10  percent  of  the  input  3^*7lS  radiation  to  an  energy  meter  to  monito 
the  dissociation.  With  removal  of  the  germanium  flat  approximately  twice 
the  former  amount  of  3^71^  radiation  traverses  the  gain  cell.  The  .aser  puls 
energy  at  2.5  microns  was  measured  and  values  were  found  in  the  range 
of  0.001  to  0.002  J as  the  doubled  ruby  pulse  varied  from  0.2  to  0.25  J. 

Laser  lines  that  were  observed  with  the  second  configuration  are  listed  in 
Table  6-3  along  with  relative  values  for  the  strength  of  each  line.  In  con- 
trast to  the  first  tests,  laser  oscillation  was  observed  on  transitions  from 
v(l-0)  to  v(5-L)  with  a major  part  of  the  energy  accounted  for  by  a few 
v(l-0)  lines.  It  appears  from  these  results  that  additions*  lines,  including 
those  of  v(6-5  ) , would  also  be  observed  if  more  energy  in  the  doubled  ruby 
pu.se  were  available. 
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FIGURE  41  PULSED  LASER  OUTPUT 


TABLE  8.  PUISED  LASER  TEST  DATA 


A.  Tests  with  Germanium  Mirror 


Code 

Pulse 

Energy 

Line 

Relative 

Observed  Signal  Strength 

Laser 

Run  #1 

Mix  .5 

torr  Uq,  1.0  torr  Fg,  98*5  torr  Ar 

AA 

.071 

1P4 

* 1.9 

BB 

•> 

1P3 

* 6.4 

CC 

.075 

1P2 

DD 

.070 

LP5 

EE 

.077 

1P6 

pm 

.068 

2P3 

* 5.0 

FF2 

.075 

2F3 

* 6.5 

GG 

.075 

2P4 

* 3.0 

HH 

.075 

2P5 

II 

.068 

2P6 

JJ 

.080 

2P7 

KK 

.075 

2P2 

Laser 

Run  #2 

Mix  l.C 

1 torr  Hg,  .5  torr  Fg,  98.5  torr  Ar 

AA1 

.07 

1P3 

3B1 

.079 

1P4 

CC1 

.076 

2P3 

* 1.5 

DD1 

.080 

2P4 

EE1 

.076 

1P3 

* .9 

Laser 

Run  # 3 

Mix  .5 

torr  Ho,  .5  torr  Fg,  99  torr  Ar 

AA2 

.083 

1P3 

* .5 

332 

.077 

2P3 

* .45 

CCS 

.079 

1P4 

DD2 

• 

O 

ON 

2P4 

TABLE  8.  PUIS  ED  LASER  TEST  DATA  (COUT'D) 


Pulse  Relative 


Code 

Enersry 

Line  Observed 

Sisnal  Strength 

Laser 

Run  #L 

Mix  .2  torr  Hq,  1.0  torr 

Fg,  98.8  torr  A r 

AA3 

• 075 

1P3  * 

2.0 

3B3 

• 077 

1P2 

CC3 

.086 

lPL  * 

1.75 

DD3 

.075 

1P5 

EE  3 

• 07L 

2P3  * 

6.5 

FF3 

.075 

2PL  * 

1.25 

GG3 

.077 

2P5 

HH3 

.075 

1P2 

AAL 

.072 

3-2  PL 

BBL 

.079 

3-2  P3 

CCL 

.082 

3-2  P5 

DDL 

.078 

3-2  P8 

EEL 

.07L 

L-3  PL 

F^L 

.07L 

L-3  P6 

Laser 

Run  #5 

Mix  .5  torr  HI,,  1 torr  F2 

, 98.5  torr  Ar 

LrJ. 

.07 

1PL  * 

• 50 

LA2 

.C65 

1P5  3-2 

LA3 

.068 

PL( 3-2 ) 

LAL 

.062 

P3(3-2) 

LA  5 

.068 

P3(L-3) 

LA 6 

.065 

PL(L-3) 

LA7 

.065 

P5(L-3) 

LA6 

.068 

2P3 

83 


TA3LS  8.  PUIS  ED  LASER  TEST  DATA  (COUT’D) 


Code 

Pulse 

Energy 

Line 

Observed 

Relative 
Signal  Strength 

LA9 

.069 

P3(2-l) 

LAIO 

.06? 

P3(2-l) 

* 

• 75 

LA11 

.067 

P3(2-l) 

LA12 

• 07 

P4(2-1) 

* 

• 65 

LA13 

.068 

P5(2-l) 

LA14 

.065 

P5(2-l) 

* 

• 52 

LA15 

O 

CT\ 

03 

?6(2-l) 

LAlb 

.065 

P7(2-l) 

LA17 

.063 

P3( l-o) 

' 

LA18 

.06 

P4(  1-0 ) 

LA19 

.06 

P4(2-l) 

* 

• 35 

LA20 

.06 

P4(  1-0 ) 

* 

.82 

LA21 

.06 

P4(l-0) 

* 

1.5 

LA22 

.06 

P5(l-0) 

3.  Test: 

3 without 

Germanium 

Mirror 

Gas  Mixture  0.5  toj 

rr  Hg,  1.0 

torr  Fo, 

98.5 

torr  Ar 

Pulse 

Code 

Energy 
( Joules ) 

Line 

Signal  Strength 

v(  1-0) 

A1 

0.19 

?3 

400 

A2 

0.20 

P4 

500 

A3 

0.18 

P5 

A4 

0.15 

P6 

20 

34 


:« * * ' r>  *• 


d 


TABLE  8.  PULSED  LASER  TEST  DATA  (COUT'D) 


Code 

Raise 
Energy 
( Joules ) 

Line 

Signal  5t 

A5 

0.18 

F7 

40 

A6 

0.20 

P8 

2 

A7 

0.21 

P9 

1 

v(2-l) 

31 

0.20 

P4 

1.0 

32 

0.24 

P5 

0.7 

B3 

0.20 

P6 

2.0 

34 

0.20 

P7 

1.5 

v(3-2) 

Cl 

0.19 

P5 

0.5 

C2 

0.18 

P6 

0.3 

C3 

0.22 

P7 

— 

v(4-3) 

D1 

0.22 

P2 

0.05 

D2 

0.27 

P3 

0.1 

D3 

0.24 

P4 

0.05 

D4 

0.18 

P5 

— 

v(5-4) 

El 

0.29 

P2 

0.02 

E£ 

0.20 

P3 

— 

E3 

0.22 

P4 

— 

S4 

0.20 

P5 

... 
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TABLE  8.  PULSED  LASER  TEST  DATA  (COIJT'D) 


Code 

Pulse 
Energy 
( Joules ) 

Line 

Signal  Strength 

v(6-5) 

FI 

0.20 

P2 

— 

F2 

0.18 

P3 

— 

F3 

0.20 

P4 

— 

F4 

0.22 

P5 

— 

?5 

0.19 

P6 

— 

SECTION  IV.  DATA  ANALYSIS 


As  data  from  the  flow  tube  and  pulsed-fluorescenee  experiments  became 
available,  an  ongoing  analysis  effort  was  performed.  The  purpose  was  to 
obtain  the  best  self-consistent  set  of  kinetic  rate  coefficients  that  would 
fit  the  results  from  both  experiments.  Such  a rate  package  would  tne.n  have 
application  in  predicting  and  optimizing  the  performance  of  chain-reaction 
(and  also  cold-reaction)  HE  chemical  lasers.  The  results  of  this  analysis 
are  given  below. 


Chemical  Reactions 

At  the  start  of  this  contractual  effort  an  extensive  parameter  study  was 
performed  to  identify  the  major  problem  areas  related  to  predicting  perfor- 
mance of  HE  chain  reaction  lasers.  It  was  concluded  that  among  the  possible 
reaction  mechanisms,  the  key  uncertainties  were  in  the  forward  rates  of  the 
pumping  reactions  (F  + H2  -»  HE*  + H,  H + Fp  — HE^  + F).  Also  of  possible 
importance  were  wall  interactions  (at  least  in  the  envisioned  experiments) 
and  the  back  cold  reaction  (HF(v  a 3)  + H -•  H2  + F).  Other  reaction  mechan- 
isms - recombination,  chain  branching,  hot  atom  effects  - were  found  to  be  of 
secondary  importance.  Furthermore , such  effects  could  not  be  easily  iden- 
tified in  the  available  experiments. 

The  findings  from  the  experiments  have  not  contradicted  the  preliminary 
conclusions.  While  it  is  possible  that  hot  atom  effects  might  influence  the 
pulsed-fluorescenee  experiments,  it  appears  that  the  basic  hot  and  cold  reac- 
tion mechanism  is  predominant.  Some  weak  evidence  of  influence  from  one  back 
cold  reaction  is  present,  however,  and  (for  these  experiments  at  least1)  it 
seems  the  influence  of  the  test  cell  walls  cannot  be  completely  ignored. 

Flow  Tube 

Analysis  of  the  flow  tube  data  was  first  aimed  at  the  evaluation  of  the 
hot  reaction  rate  coefficient,  which  it  was  thought  might  be  slower  t..ar 
previously  reported  (see  Appendix  A for  a description  of  the  initial  rate 
package).  Four  test  conditions  were  considered  with  varying  proportions  of 
Ep  ar.d  Hg  (see  Table  $).  The  corresponding  HE  production  is  shown  in 
Figure  h2. 

in  the  Hp  rich  cases  the  rate  determining  s^ep  is  the  hot  reaction.  Thus 
analytic  predictions  of  the  production  of  HF  for  such  cases  will  be  ncs~  sensi- 
tive to  the  value  used  for  the  hot -reaction  rate  coefficient.  Tf  one  assumes 
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TABLE  9.  INITIAL  CONDITIONS 
Flow  Tube  Simulation 


Velocity  = 10000  cm/sec 
Pressure  = 20  Torr 

Species  Concentration, 

3 . 

xF2 

XHg 

XAr 

0.0001 

0.00435 

0.00055 

0.995 

0.0001 

0.00214-5 

0.00C775 

0.99675 

0.0001 

0.00095 

0.00155 

0.9974 

0 

b 

0 

25 

0. 0002 5 

0.00310 

0.99655 

0.0001 

0.00172 

0.00C88 

0.99731 

0.00026 

O.OOC59 

0.00088 

0.99327 

FLOW  TUBE  EXPERIMENT 


FIGURE  42  VARIATION  OF  TOTAL  HF 


-hat  '.vail  interactions  and  the  back  cold  reaction  are  unimportant , a 
preliminary  value  for  the  hot  reaction  rate  coefficient  can  be  obtained 
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Appendix  3' 


J Hn^/dt 


for  the  Hg  rich  cases  (3  and  4)  and  using  the  formula 


(at 


i » 


n h 


dt 


or^HF 


2(nv)  jCcnHo  - — 


i’esuj.t  o: 


ne 


tea 


cate  approximation. 
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approach,  the  value  for  the  hot  reaction  rate  coefficient  is  found  to  be 
approximately  one  half  the  commonly  accepted  value  of  Albright  (Hef.  20) 
et.  al.  fi.e.,  1.  : x 1C“ — : cc/molecule-sec  . 


(1) 


.he  quasi-one-dimens ional  computer  code  CKFMLZ  was  then  used  to  test  the 
validity  of  such  a rate.  To  account  for  two  dimensional  boundary  layer  effects 
a correction  factor  is  employed  in  the  calculation.  This  is  described  in 
detail  in  Appendix  C.  Initial  conditions  for  the  calculations  are  given  in 
Table  9.  A comparison  of  the  model  calculations  (Fit  F)  with  the  experimental 
results  are  shown  in  Figures  4-3  through  46.  It  can  be  seen  that  the  comparison 
is  reasonable  for  the  Hp  rich  cases  (Conditions  3 and  4) , but  for  the  Fp  rich  cases 
(Conditions  1 and  2),  the  predicted  increase  of  HF  with  distance  fails  below 
the  experimental  curves.  Use  of  the  nominal  value  for  the  hot  reaction  rate 
coefficient  improves  the  agreement  for  Conditions  I and  2 at  the  expense  of 
the  good  agreement  for  Conditions  3 ar.d  4. 


liote  that  in  all  cases  the  experimental  curves  start  off  mere  steeply 
than  the  model.  In  drawing  the  curves  through  the  experimental  points,  it  was 
assumed  that  at  the  start  cf  the  flow  tube  (i.e.,  X = 0)  there  is  no  HF.  The 
flow  in  the  vicinity  of  X = 0,  however,  is  quite  complicated;  recirculation 
zones  might  exist.  In  such  a case  it  would  be  more  appropriate  to  model  she 
experiments  assuming  a finite  initial  value  for  the  HF  concentration. 


If  hue  constraint  that  n-qp  = 0 at  X =0  were  relaxed  it  is  quite  reason- 
able to  extrapolate  the  experimental  carves  back  to  nonzero  initial  values. 
This  is  illustrated  by  the  dashed  lines  on  Figure  42.  A comparison  of  these 
modified  curves  with  the  model  is  shown  in  Figures  47  through  50.  The  initial 
concentrations  for  the  calculations  are  taken  equal  to  the  extrapolated  values 
from  the  experimental  curves. 


'.law  the  slopes 
curves  compare  cons 
'.he  poorest  fit  is 


and  magnit udes  of  -...e  exoerinen tal  and  theoretical  ’i"  t 
ideraoly  better  than  when  zero  initial  values  are  assumed, 
for  run  condition  4.  This  might  indicate  that  a value  for 
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Albright,  R.  3.,  et  al.  ” Mas s -3 pec t rone trie 
Constants  for  H-Atom  Reactions  with  Clp  and 
'•*32,  1969. 


Determination  of  Rate 
Fo . " .’hem.  Fhyo.,  50, 
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HF  VARIATION  IN  FLOW  TUBE 


FIGURE  43  COMPARISON  BETWEEN  MODEL  AND  EXPERIMENT 


HF  VARIATION  FLOW  TUBE 

RUN  CONDITION  2 


FIGURE  44  COMPARISON  BETWEEN  MODEL  AND  EXPERIMENT 
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FIGURE  46  COMPARISON  BETWEEN  MODEL  AND  EXPERIMENT 


HF  VARIATION  IN  FLOW  TUBE 


FIGURE  47  COMPARISON  BETWEEN  MODEL  AND  EXPERIMENT 


HF  VARIATION  IN  FLOW  TUBE 


FIGURE  48  COMPARISON  BETWEEN  MODEL  AND  EXPERIMENT 
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FIGURE  50  COMPARISON  BETWEEN  MODEL  AND  EXPERIMENT 


Gn  the  other  hand,  a better  match  might  also  be  achieved  if  the  true  initial 
HF  concentration  were  smaller  than  assumed.  While  it  appears  plausible  from 
the  above  result  that  some  HF  is  present  initially,  in  the  absence  of  any 
measurements  at  X = 0 it  cannot  be  conclusively  proven. 

An  alternative  explanation  to  the  experimental  results  might  be  that  the 
steep  initial  slopes  do  indeed  represent  faster  rates  for  the  hot  and/or  cold 
reactions , and  the  shallower  slopes  later  on  result  from  wall  interactions 
(in  particular  the  loss  of  H or  F atcms).  Support  for  such  an  approach  is 
provided  by  the  pulsed-fluorescence  experiments,  which  indicated  faster  rate 
coefficients  for  both  the  hot  and  cold  reactions.  (This  will  be  described 
in  more  detail  below. ) 

The  rates  that  apparently  fit  the  pulsed  experiments  best  are: 

kc  = 2.0  x nominal 
]a  a 1.3  x nominal 

These  rates  were  used  in  a flow  tube  calculation  in  which  it  was  further 
assumed  that  all  H atoms  colliding  with  the  wall  were  lost  (to  later  recom- 
bine to  form  Ho).  This  H-atom  loss  was  accounted  for  approximately,  based 
upon  estimates  arising  from  an  analysis  described  in  Appendix  D.  The  loss 
of  H atoms  was  expressed  as 


:k:  wans 


= _ rw’ 


,-9‘ 


(2) 


with  3 = 156.  This  value  was  obtained  from  a fit  of  the  curve  in  Figure  Dl. 
over  the  time  scale  of  the  experiment. 


A comparison  of  this  set  of  rates  (Fit  P2)  is  also  shown  in  Figures  43  through  46. 
As  can  be  seen,  even  with  the  inclusion  of  potential  H-atom  less  such  a set 
of  rates  overpredicts  the  HF  population  for  Conditions  2 through  4 in  all  but  the 
initial  portions  of  the  flow  tube. 


The  final  comparison  made  is  for  a fit  intended  to  give  good  general 
agreement  in  both  experiments.  The  rates  employed  for  the  hot  and  cold  reac- 
tions are 


k = 1.5  x nominal 
c 

k,  = 0.9  x nominal 
n 


oa 
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Again  the  loss  of  H-atoms  to  the  walls  was  assumed.*-  The  results  of  the 
comparison  are  shown  in  Figures  43  through  46  as  fit  Cl.  We  see  that  such  a fit 
matches  the  data  reasonably  well  although  the  slope  of  the  curves  is  a bit 
too  shallow  at  first  and  too  steep  later  on.  Use  of  such  rates  in  a comparison 
assuming  finite  initial  HF  concentrations  is  shown  in  Figures  47  through  50.  The 
fit  is  a bit  poorer  and  tends  in  all  cases  to  overestimate  the  HF  populations. 

Pulsed-Fluorescence 


These  experiments  were  run  for  six  different  sets  of  conditions,  which 
are  outlined  in  Table  10.  The  time  history  of  the  production  of  HF  is 
available  from  oscilloscope  traces  of  the  absorption  of  a(v=l)  -»  (v=0)  cw 
laser  beam.***  Attempts  to  fit  the  data  using  the  nominal  (Homann)  value 
(Ref.  21)  for  the  cold  reaction  rate  coefficient  and  half  the  nominal  value 
(Ref.  20,  Albright)  for  the  hot  reaction  rate  coefficient  as  had  been  done 
for  the  flow  tube  met  with  no  success.  The  predicted  temporal  development 
of  HF  fell  below  the  measured  values  in  virtually  all  cases.  The  general 
indication  was  that  faster  hot  and  cold  reaction  rate  coefficients  would  be 
required  to  fit  the  majority  of  the  data. 

In  an  attempt  to  facilitate  the  analysis  of  the  data  from  each  test 
series,  composite  plots  were  made  up  employing  a time  normalized  to  the 
pulse  energy  (i.e.,  t = t x E^).  The  rationale  behind  this  is  that  for  a 
quasi-steady-state  situation  in  the  absence  of  any  loss  mechanisms  all  the 
cases  would  fall  along  the  sam<=  '"urve . 

Consider  the  equation 


dn 


HF  = ^r-Fo  "C  kh  nF2  nH2 


dt 


(3) 


^c  nH2  + khnF2 


and  also 


F0  P 


(4) 


* This  had  an  effect  ranging  from  0 to  25  percent 

**  For  the  conditions  of  these  tests  the  contribution  from  upper  levels  to 
the  total  HF  population  rapidly  becomes  negligible. 


21.  Homann,  H.,  et  al. , "Fine  Methode  Zur  Snzenvang  von  Fluorstonier 
in  Ir.erten  Atmos phare,”  Per.  Bunse.  F.  Phys.  Them.,  '4,  555,  1970. 
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TABLE  10.  INITIAL  CONDITIONS 


Pulsed  Fluorescence  Simulation 
Pressure  = 100  Torr 
Species  Concentrations: 


Species 

xf2 

xHo 

xAr 

A 

.0094 

..0102 

.9804 

B 

.0102 

.0222 

.9676 

C 

.0099 

.0055 

.9846 

D 

.0049 

.0055 

.9896 

2 

.0049 

.0102 

.9849 

F 

.0101 

.0024 

.9875 

101 


J 


Thus  we  my  write 


^"HF  ^ *^h  n?2  ^Hp  (c'j 

dt  nH2  + kh  “F2 

Here  the  right  hand,  side  is  not  explicitly  dependent  on  pulse  energy.  The 
presence  of  heat  transfer  effects,  which  will  be  different  for  different  pulse 
energies,  and  of  wall  interaction  effects  will  cause  some  deviation  in  the 
curves  from  case  to  case.  However,  it  is  not  unreasonable  to  expect  all  the 
results  from  a given  series  to  start  off  with  essentially  the  same  slopes,  with 
the  higher  energy  cases  being  perhaps  slightly  steeper. 

The  resulting  composite  plots  are  shown  in  Figures  51  through  56.  The  results  for 
series  A,  B,  and  C show  relatively  little  scatter  near  t = 0,  although  it  is 
difficult  to  find  any  correlation  between  initial  slope  and  the  pulse  energy. 
Series  D also  gives  relatively  consistent  results  with  the  possible  exception 
of  DH5  and.  DH3.  On  the  other  hand,  the  series  E and  F results  show  consider- 
able scatter.  As  a consequence,  the  results  from  series  A,  B,  and  C have 
been  weighted  most  heavily  and  those  of  E and  F least  heavily  in  attempting  to 
find  a good  overall  fit  to  the  data. 

A complication  that  renders  data  interpretation  difficult  for  the  pulsed 
fluorescence  experiments  is  the  apparent  adsorption  of  HF  on  the  test  cell 
walls.  It  is  possible  to  observe  the  fall  off  of  HF  number  density  in  the 
latter  stages  of  many  tests.  Also,  the  maximum  number  density  produced  in 
any  given  test  tends  to  fall  below  the  expected  maximum  based  on  the  initial 
concentrations  of  Hp  and  Fp . Furthermore,  there  is  a correlation  between  the 
pulse  energy  of  a given  test  and  the  maximum  HF  produced  (Figures  57  through  62 ) , It 
can  be  seen  that  higher  pulse  energies  (corresponding  to  high  initial  F-atom 
concentrations  and  faster  reactions)  result  in  high  HF  concentrations.  Such 
a behavior  would  occur  if  there  were  competition  between  the  walls  and  the 
chemical  reactions.  Finally,  an  analysis  to  estimate  the  worst  possible  loss 
of  HF  to  the  walls  (Appendix  D)  shows  that  the  amounts  that  might  be  lost 
to  the  walls  are  not  unreasonable.  As  a result  of  these  considerations  a 
simple  model  including  the  loss  of  HF  to  the  wall  was  used  to  first  try  to  fit 
some  of  the  experimental  data.  It  was  assumed  that  the  adsorption  of  HF  on 
the  walls  was  proportional  to  the  room  available  on  the  walls,  (l  - r--HF  'W  ) } 

where  nvy  by  is  the  equivalent  number  density  of  HF  already  on  the  walls  and 
nHF I FAX  the  maximum  allowed.  The  adsorption  is  also  taken  to  be  propor- 
tional to  the  HF  available  in  the  ;as  phase  np-?|  . f?.ef.  22).  Thus  we  have 


22.  Moore,  W.  J.,  Physical  Chemistry,  Prentice-Hall,  Englewood  Cliffs 
HJ,  1955. 
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FIGURE  52  COMPOSITE  HF  TEMPORAL  BEHAVIOR  NORMALIZED  TIME 
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PULSE  CONDITION  D 
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FIGURE  54  COMPOSITE  HF  TEMPORAL  BEHAVIOR 


PULSE  CONDITION  E 
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FIGURE  55  COMPOSITE  HF  TEMPORAL  BEHAVIOR 
NORMALIZED  TIME 


PULSE  CONDITION  F 
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FIGURE  56  COMPOSITE  HF  TEMPORAL  BEHAVIOR 
NORMALIZED  TIME 


TEST  SERIES  B 
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FIGURE  58  VARIATION  OF  MAXIMUM  HF  PRODUCED  WITH  PULSE  ENERGY 
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FIGURE  60  VARIATION  OF  MAXIMUM  HF  PRODUCED  WITH  PULSE  ENERGY 

TEST  SERIES  DH 


FIGURE  61  VARIATION  OF  MAXIMUM  HF  PRODUCED  WITH  PULSE  ENERGY 

TEST  SERIES  EF 


TEST  SERIES  F 
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FIGURE  62  VARIATION  OF  MAXIMUM  HF  PRODUCED  WITH  PULSE  ENERGY 
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V|W  = - 

The  overall  production  of  HF  is  approximated  by  the  steady  state  formula 

2 ( np ) 0 nF  nH  1%  kc 


^HF I chem 


Vf2  + kcnH2 


(7) 


and  the  net  change  of  HF  in  the  gas  phase  is  then 


Vi  ~ V| 


HF 

chem  ~ 1 W 


(3) 


A simple  integration  of  the  above  equations  for  several  experimental  conditions 
was  performed  to  see  if  one  set  of  rate  parameters  could  be  ’used  to  reasonably- 
fit  a number  of  cases  (A6-S,  33-5,  C6-9,  DH1-2).  A reasonable  fit  was  obtained 
for  the  following  set  of  coefficient  values: 

k^  = 1.3  x nominal 
kc  = 2.0  x nominal 

k = 200  sec"3- 

w , 

. _ ,„l6  -3 

“hfImax  = 1,5  x 1°  cm 


The  results  of  this  exercise  are  given  in  the  contract  test  report  (Ref. 
23).  One  problem  with  this  fit  to  the  data  was  that  she  experimental  curves 
tend  to  bend  over  more  gradually  than  the  model  predicts  as  the  reaction 
nears  completion.  A possible  explanation  for  such  behavior  could  be  the  loss 
of  H atoms  to  the  wall.  An  examination  of  that  possibility  (Appendix  D)  has 
indicated  a potentially  significant  effect.  Hence  the  model  was  modified  to 
account  for  the  loss  of  H atoms  in  a manner  as  described  in  Appendix  B (see 
Eq.  317).  However,  the  form  of  kn  is  taken  to  be 

k^j  = Atm-1  (9) 

since  we  may  better  fit  the  analytic  curve  (Appendix  D,  Figure  32)  for  the  loss 
of  H atoms  to  the  form 


“H  = "Hq 


“A.t.n 


(10^ 


For  complete  loss  of  H atoms  the  coefficients  A and  m would  have  the  values 
17.325  and  0.603  respectively. 

i ■.  3oedeker,  L.  R. , et  al. , "Long-Range  Chain  Laser  Studies  Test  Report," 
for  Air  Force  Contract  F29601-75-C-0139,  UTRC  R76-952 198-4?,  East 
Hartford,  CT,  October  1976 • 
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The  modified  model  was  then  used  in  extensive  comparisons  with  experiment 
Using  the  composite  plots  (Figs.  51-56)  for  guidance,  representative  cases 
from  each  test  series  were  selected  for  analysis.  These  cases  are  listed  in 
Table  11  along  with  the  initial  conditions  used.  The  values  selected  for 
HFW  are  designed  to  res -alt  in  a final  HF  number  density  near  that  observed 
experimentally.  The  rationale  is  that  the  wall  may  be  saturated  with  KF  to 
different  degrees  from  test  to  test. 

Some  of  the  results  of  the  analysis  are  summarized  in  Figures  63  through  74.  The 
curves  labeled  Fif  ? is  for  a fit  that  gave  the  best  overall  fit  to  the 
different  test  series.  The  values  for  the  various  rate  coefficients  are: 

k^  = 1.3  x nominal 
kc  = 2.0  x nominal 
kw  = 100  sec_1  l6  , 

°HF 1 MAX  1-5  x 10  cm 

A = 8.0 
m = O.608 

The  value  of  8.0  for  A implies  an  H atom  adsorption  accommodation  coeffi- 
cient of  about  0.5-  It  is  possible  to  juggle  the  values  of  all  the  wall 
interaction  parameters  to  obtain  somewhat  better  fits  but  further  effort  was 
not  deemed  worthwhile. 

As  can  be  seen  from  the  curves,  the  fit  is  best  for  series  A,  B,  and  C. 
This  is  by  design  since  it  is  felt  that  those  are  the  most  reliable  data. 

The  relatively  fast  ra^e  for  the  cold  reaction  is  dictated  primarily  by  the 
initial  slope  from  series  C,  whereas  the  somewhat  fast  hot  reaction  rate  was 
indicated  by  series  A and  B for  which  'rhe  hot  reaction  is  the  rate  controlling 
step.  Note  that  even  with  the  inclusion  of  H atom  effects  it  was  not  possible 
to  match  the  observed  bending  over  of  the  series  C curves.  Much  faster  *han 
reasonable  wall  interaction  and  chemical  rates  would  be  necessary  to  provide 
a better  fit;  and  that  would  be  at  the  expense  of  good  fits  for  series  A and 
B. 

Further  no'-.e  ‘•hat  a single  se4-  of  rates  could  no4-  really  fit  all  the 
series  simultaneously.  The  predictions  for  series  D fall  consistently 
below  the  experimental  curves.  This  is  even  more  dramatically  true  for 
series  E.  On  the  other  hand,  those  same  rates  predict  curves  well  above 
the  measured  curves  for  series  F. 

The  influence  of  reverse  chemical  reactions  might  be  invoked  as  an 
exnlana4- ion  of  this  behavior,  but  inclusion  of  such  effects  only  showed  a 
."4rortf  influence  for  unreasonably  fast  ra^es  (greater  than  'as  kine4ic)  and 
wcul  i imply  •'as4  er  forward  rates  as  well. 


TABLE  11.  INITIAL  CONDITIONS  - PULSED  FLUORESCENCE  EXPERIMENT 


Case 

F 

F2 

H2 

^g 

HF 

w 

A6 

3.51  x 1013 

2.87  x 1016 

3.12  x 10l6 

0.3  X IO10 

O.65  x 10 

A7 

2.33 

2.67 

2.92 

0.7 

1.0 

AS 

3.28 

2.77 

3.02 

0. 

• 5 

0.7 

AlO 

3.35 

2.85 

3.10 

0. 

• 35 

1.1 

A15 

3.16 

2.35 

3.10 

0. 

■ 35 

1.5 

B3 

2.32 

2.94 

6.80 

0. 

70 

1.0 

B4 

4.66 

2.99 

6.85 

0.60 

1.2 

B5 

1.88 

2.79 

6 .65 

1.00 

0.3 

c6 

3.38 

3.173 

1.77 

0. 

0 

1.25 

C7 

3.55 

I 

I 

1.0 

C8 

5.49 

1 

1 

0.9 

C9 

5.09 

♦ 

f 

0.9 

DKL 

2.66 

1.58 

1.77 

0.0 

DH2 

2.54 

1 

1 

0.19 

DH7 

1.76 

T 

1 

0.4 

EE  2 

2.05 

1.58 

3-02 

1.0 

SF5 

2.69 

1 

1 

1.32 

EF6 

2.35 

1 

r 

1.5 

F7 

4.35 

3.25 

0.76 

1.0 

FIO 

6.21 

i 

i 

f 

1.3 

Fll 

6.04 

T 

T 

f 

1.1 

11' 


PULSE  CASE  A6 


o 


(gi_0t*  e-W3)  JHu  'A1ISN3Q  d39IAIDN 


ne 


FIGURE  63  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
OVERALL  HF  DEVELOPMENT 
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FIGURE  64  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
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FIGURE  65  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
OVERALL  HF  DEVELOPMENT 
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FIGURE  67  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


OVERALL  HF  DEVELOPMENT 
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FIGURE  68  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


FIGURE  69  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


OVERALL  HF  DEVELOPMENT 
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FIGURE  7C  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


OVERALL HF  DEVELOPMENT 


FIGURE  71  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
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FIGURE  72  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
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FIGURE  73  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
OVERALL  HF  DEVELOPMENT 
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FIGURE  74  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
OVERALL  HF  DEVELOPMENT 


Also  shewn  in  Figures  63  through  74  are  curves  for  the  same  chemical  rates  as 
above,  bu'-  without  any  wall  effects  (Fit  Pi)  and  ‘‘he  combined  rates  meant  to 
fit  both  se'-s  of  experiments  (Fi4-  C),  also  with  no  wall  effects.  We  see 
tha-  '"he  initial  slopes  for  Fit  PI  still  give  a reasonable  fit  to  the  data, 
bu4-  '•ha4'  naturally  the  final  HF  populations  are  higher  than  observed  experi- 
mentally. It  is  important  to  note  that  the  absence  of  the  inclusion  wall 
effects  does  no*-  necessarily  imply  slower  chemical  rates. 

An  alternative  approach  can  also  be  used  in  interpreting  the  pulsed- 
fluorescence  data.  In  utilizing  the  probe  laser  there  is  some  problem 
with  stability,  both  in  maintaining  a constant  input  intensity  and  remaining 
on  line  center.  Both  the  intensify  and  frequency  can  drift  from  run  to  run, 
and  during  the  course  of  a run.  Despite  the  care  that  was  taken  to  minimize 
these  effects  it  is  possible  that  they  have  influenced  the  data  i.nterpretatio 

Consider  the  following  example.  Two  runs,  Cll  and  C12,  had  the  same 
pulse  energy  (0.134J)  and  supposedly  the  same  initial  reactant  concentrations 
The  curves  of  number  density  calculated  from  the  probe  laser  absorption 
traces,  however,  are  dramatically  different.  These  are  shown  in  Figure  75* 

One  might  suppose  that  somehow  the  gas  mixture  changed  and  we  are 
actually  observing  two  different  HF  temporal  histories.  However,  observation 
of  the  measured  temperature  histories  for  the  two  cases  (Fig.  76)  implies 
4'ha4‘  we  are  actually  observing  ‘•he  same  reaction.  Furthermore,  normalization 
of  ‘•he  KF  traces  to  unity  (Fig.  77)  yields  essentially  the  same  curve  in 
both  cases. 

While  the  above  example  is  one  of  the  more  dramatic  cases,  it  lends 
strong  support  to  the  hypothesis  that  much  of  "-he  variation  in  final  HF  concen- 
tration that  has  been  observed  is  attributable  to  drifting  in  the  cw  laser 
signal. 

Consequently,  i4:  was  decided  vo  reanalyze  the  da'-a  under  the  assump- 
tion ‘■hat  no  HF  was  lost  to  the  walls  and  that  all  cases  within  a given 
series  produced  the  same  ultimate  HF  concentration,  namely  the  maximum 
allowable.  Samples  of  the  results  of  the  analysis  wi4h  these  adjusted  HF 
populations  are  shown  in  Figures  78  through  89.  Two  fits  are  shown,  PI  and  C, 
which  have  been  previously  described.  Notice  •■hat  now  fit  ?1  compares 
well  with  series  F as  well  as  A,  B,  and  C.  The  curves  for  fi4*  C tend  to  be 
consistently  below  the  da*-a,  though  perhaps  tolerably  so. 

The  ma^ch  with  series  D and  E remains  moor.  The  experimental  curves 
lie  well  above  ‘he  oredicted  curves.  The  one  thing  tha4-  *-hese  series  have  in 
common,  to  ‘he  exclusion  of  all  others,  is  a lower  7,  concentration  (l  4 err 
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FIGURE  75  COMPARISON  OF  HF  HISTORIES 
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FIGURE  76  COMPARISON  OF  TEMPERATURE  HISTORIES 
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FIGURE  78  COMPARISON  OF  MODEL  WITH  EXPERIMENT  ADJUSTED  HF  POPULATIONS 
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FIGURE  79  COMPARISON  OF  MODEL  WITH  EXPERIMENT  ADJUSTED  HF  POPULATIONS 


ADJUSTED  HF  POPULATIONS 


FIGURE  80  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
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ADJUSTED  HF  POPULATIONS 
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FIGURE  82  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


ADJUSTED  HF  POPULATIONS 
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FIGURE  83  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


ADJUSTED  HF  POPULATIONS 


FIGURE  84  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


ADJUSTED  HF  POPULATIONS 


FIGURE  86  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


ADJUSTED  HF  POPULATIONS 
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FIGURE  87  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


ADJUSTED  HF  POPULATIONS 


O 


lUU 


FIGURE  88  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


ADJUSTED  HF  POPULATIONS 


FIGURE  89  COMPARISON  OF  MODEL  WITH  EXPERIMENT 


vs.  1 -err' . If  fhe  laJ  a are  to  be  believed  this  implies  a nonlinear  depen- 
dence of  the  overall  HF  production  rate  on  F0  concentration,  in  fact  less  tha 
linear.  A completely  different  reaction  mechanism  than  the  usual  chain  would 
be  necessary.  For  example,  Fo  might  be  playing  the  role  of  an  inhibitor  for 
seme  branching  or  ho1-  a^om  effect  *"hat  is  less  important  as  its  concentration 
decreases . 

“here  is  no4:  enough  evidence  in  these  experiments  alone  to  strongly 
supper4-  such  possibilities,  although  it  is  a question  that  should  be 
addressed  in  the  future.  For  now  we  are  inclined  uo  lack  confidence  in 
the  series  D and  E data. 

The  rates  indicated  for  the  pulsed  fluorescence  tests  are  higher  than 
can  be  reasonably  used  to  fit  fhe  flow  tube  data.  Thus  it  is  instructive  t0 
try  to  find  a rational  explanation  for  that  finding  that  might  reconcile  the 
iif  • 

One  might,  infer  that  perhaps  the  initial  F-atom  concentration  is  somehow 
in  error.  If  it  were  underestimated  the  resultant  analysis  would  yield  rates 
that  were  too  fast.  One  possibility  is  that  the  absorption  coefficient  used 
to  evaluate  the  number  of  fluorine  molecules  dissociated  is  in  error.  There 
is  no  reason  to  believe,  however,  that  the  literature  value  (Ref.  11)  would 
be  significantly  in  error.  To  invoke  a chain  branching  mechanism  would 
require  rate  coefficients  many  orders  of  magnitude  larger  than  suggested  by 
Sullivan  et  al.  (Ref.  24)  to  provide  even  a perceptible  difference  in  the 
model  predictions.  The  remaining  reasonable  explanation  is  the  possibility 
of  hot  F atom  effects. 

When  one  photon  of  3471  « radiation  dissociates  a fluorine  molecule, 
p4  percent  of  the  photon  energy  is  left  over  to  go  into  translation  of  fhe 
fluorine  atoms.  These  highly  energetic  atoms  could  be  very  efficient  in 
dissociating  more  Fo . The  large  amount  of  diluent,  on  the  other  hand,  tends 
to  moderate  these  atoms  before  they  collide  with  other  fluorine  molecules. 
Some  calculations  have  beer,  made  to  evaluate  such  an  effect  approximately,  bu 
in  a manner  likely  demonstrate  the  maximum  possible  influence  on  the  results. 
If  is  assumed  that  upon  dissociation  all  the  excess  energy  is  transferred  co 
one  of  the  fluorine  atoms  (i.e.,  half  of  the  initially  dissociated  T atoms 
are  "hot"  with  a translational  temperature  of  15C00°K; . In  a collision  of  a 
fl  j 3 r.e  moXacuXa  7'  with  a . ot  atom  hq  as  sums  a 5 rcsnl  * a - 


24.  Sullivan,  J.  H.,  et  al.,  "Mechanisms  and  Types  of  Explosive  Behavior 
In  Hydrogen-Fluoride  Systems,"  J.  Them.  Rhys.,  n2,  1714,  19~5- 


of  dissociation.  The  maximum  energy  that  can  be  transferred  to  the  internal 
nodes  of  the  ?2  is(mF2/(mF  + mF2))Sp^  = 2/3  (44.4  kcal/mole)  = 29*6  kcal/mole. 
This  is  some  3 kcal/mole  below  the  dissociation  energy  of  Fg«  Thus  our 
assumption  of  a 50  percent  probability  per  collision  is  not  unreasonable 
(and  likely  too  great).  In  competition  with  the  dissociation  collisions, 
the  hot  atoms  are  also  thermalized  by  collisions  with  other  molecules,  in 
particular  the  diluent.  The  loss  rate  of  F*  is  taken  to  be  equal  to  the 
collision  rate  of  the  F*  divided  by  the  average  number  of  collisions  required 
to  thernalise  the  15P00°K  F atom.  For  He,  seven  collisions  are  required  to 
reduce  a hot  atom's  energy  by  l/3,  and  for  Ar  three  collisions  are  necessary. 
Fluorine  atoms  with  energies  that  low  would  be  expected  to  have  quite  low 
probabilities  of  dissociating  Fq . To  be  somewhat  more  conservative  we  have 
assumed  that  1C  collisions  are  necessary  for  He  to  thermalise  hot  F atoms, 
and  five  collisions  are  necessary  for  Ar.  A calculation  was  performed  for  the 
case  0.2  torr  H2,  1.0  torr  F25  50  torr  He,  with  a pulse  energy  of  0.3  J.  The 
result  was  an  increase  of  11  percent  in  the  net  number  of  atoms,  yielding  an 
equivalent  increase  in  the  net  reaction. rate.*  Substitution  of  Ar  for  the  He 
reduces  this  effect  to  virtually  negligible,  which  is  the  primary  reason  the 
final  tests  were  nerfonr.ed  using  Ar.  Furthermore,  the  diluent  nressure  in  the 
final  fests  is  nearly  100  torr,  thus  reducing  this  ho1:  atom  effect  still 
further.  It  is  expected  that  the  influence  of  highly  energetic  F atoms  on 
the  results  of  ^he  final  series  of  tes^s  is  at  wors1  less  than  a few  percent. 
Furthermore,  any  hot  F atom  effect  of  this  sort  would  be  contrary  to  the 
trend  imolied  by  the  D and  E series . 

Hence,  in  the  absence  of  a good  new  model  tha1  would  reduce  ‘-.he  surpris- 
ingly fast  rate  coefficients  deduced,  no  one  se1-  of  rates  can  be  found  to 
give  a good  fit  to  both  the  flow  tube  and  pulsed  fluorescence  experiments . 

Thus  we  must  assign  a range  of  possible  values: 

0.5  x nominal  < k^  < 1.3  x nominal 
nominal  < kc  < 2.0  x nominal 

If  we  pick  values  in  the  middle  of  these  ranges  we  have 

kh  =-  .9  x nominal 

k =“  1.5  x nominal 
c 

wn icih  * •-.«?  v^3. lu. : _■  usQ'i  in  *■. h®  "2 ctt.d norr.  1 c ° f'it_ 


* Results  were  similar  for  a Dulse  energy  of  O.i.J. 
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Vibrational  Deactivation 


Preliminary  calculations  performed  under  Task  I indicated  that  the  most 
important  deactivation  mechanisms  in  need  of  study  were  HP  self  deactivation 
( VT  and  W)  and  H-atom  deactivation.  Support  for  this  emphasis  was  provided 
by  work  at  .Aerospace  Corporation  indicating  that  VT  deactivation  scales 
rapidly  with  vibrational  quantum  number  (Kwok,  Ref.  25),  and  that  H atoms 
are  a strong  deactivator  of  HF  for  v ^ 3 (Bctt,  Ref.  26). 

Flow  Tube 

Use  of  the  original  rate  package  (Appendix  A)  led  to  predicted 
vibrational  population  distributions  far  too  great  for  the  higher  lying 
levels.  Thus  it  was  necessary  to  revise  the  vibrational  rate  package.  The 
computer  code  CHFM1Z  was  modified  to  allow  for  arbitrary  scaling  of  all  V-V 
and  V-T  transitions  with  vibrations'.  quantum  number.  Provision  was  also 
made  for  multiquanta  V-T  transitions  for  all  colliding  species,  and  for 
alternative  single  quantum  V-V  models  in  lieu  of  that  normally  in  use. 

Modif ications  to  the  rate  package  were  guided  primarily  by  the  findings 
| of  Kwok  (Ref.  25)  and  Bott  (Ref.  26  and  27).  Kwok  has  found  that  the 

rate  of  disappearance  of  excited  HF  scales  as  v^*3.  Bott  has  found  that  the 
deactivation  of  HF  by  H atoms  is  slower  than  our  nominal  values  for  v = 1 
and  2,  but  considerably  faster  for  v = 3* 

o q 

One  rate  package  tried  employed  the  -r'-  scaling  on  all  single  quantum 
V-T  transitions.  In  addition  the  net  deactivation  rates  by  H-atoms  for 
v = 1,  2,  and  3 were  revised  to  match  Bott's  values.  The  temperat-ure 
dependence  and  relative  probabilities  to  the  various  lower  levels  sciven  by 


25.  Kwok,  M. , "Studies  of  the  Vibrational  Relaxation  of  Upper  Levels  of 
HF  (v  s 6)  in  a Large  Diameter,  Medium  Pressure  Flow  Tube,"  presented 
at  the  TriService  Chemical  Laser  Symposium,  Kirtland  AF3,  M,  February 
1976. 

2b.  3ott , J.  F.  and  R.  F.  Heidner  III,  "A  Study  of  the  Vibrational  Deactivation 
of  HF  (v  = 3)  by  H Atoms,"  presented  at  the  TriService  Chemical  Laser 
Sympos ium,  Kirtland  AFB,  HM,  February  1976. 

u".  Heidner,  R.  F.  Ill  and  J.  F.  Eott,  "Vibrational  Deactivation  of  HF  (v  = l) 
ana  DF  (r  = l)  by  H and  D Atoms,"  J.  Chem.  Phys.,  63 , 1610,  197 5 • 
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Wilkins  (Refs.  28  and  29),  however,  were  maintained.  This  required  multi- 
plying the  nominal  v = 1 rate  coefficient  by  0.0^9,  the  v = 2 rate  coefficient 
by  0.072,  and  the  v = 3 rate  coefficients  by  3*62.  On  the  assumption  that 
all  deactivation  above  v = 3 is  rapid,  all  H-atom  rate  coefficients  for 
v > 3 are  also  multiplied  by  3*62.  Similarly,  all  F atom  multiquanta  rates 
were  multiplied  by  3*3  to  bring  the  1-0  transition  rate  into  agreement  with 
the  value  suggested  by  Cohen  (Ref.  30). 

Application  of  these  changes  to  the  rate  package  yielded  surprisingly 
good  agreement  on  the  vibrational  distributions  for  the  four  basic  cases 
(l-h)  of  the  flow  tube  experiment.  Somewhat  better  agreement  was  achieved 
by  modifying  the  V-T  formalism  slightly.  Consider  a case  in  which  transitions 
from  a level  v are  equally  probable  to  all  lower  levels  (v-1,  v-2,  ...,  0). 
Suppose  each  individual  transition  scales  linearly  with  v.  Then  the  net  rate 
of  disappearance  from  each  level  would  be  proportional  to  v2.  The  net  loss 
rate  of  quanta,  however,  would  be  considerably  greater  than  for  a single 
quantum  process. 

Such  a multiquanta  model  was  applied  to  all  collision  partners  except  H 
and  F atoms,  which  already  included  multiquanta  transitions  with  probabilities 
apportioned  according  to  the  theoretical  work  of  Wilkins.  The  net  rates  for 
the  H atom  interactions  were  again  revised  to  reflect  Bott's  measurements  as 
described  above. 

To  achieve  the  best  fit  with  the  model  as  formulated  in  the  computer 
code  it  was  found  that  the  V-V  model  normally  employed  could  be  used  if  the 
transitions  were  scaled  with  v"^“  rather  than  v.  With  such  scaling,  exo- 
thermic transitions  of  the  type  (v,v')  -»  (v-1,  v'+l)  are  roughly  independent 
of  v and  diminish  slowly  with  v'.  This  is  a consequence  of  changes  in  energy 
defect  for  the  various  transitions.  In  summary  the  rate  package  selected  as 
the  best  fit  to  the  flow  tube  data  includes  the  following  changes:  H atom 


28.  Wilkins,  R.  L. , "Monte  Carlo  Calculations  of  Reaction  Rates  and  Energy 
Distributions  Among  Reaction  Products  H + HF(v)  -vl^v')  + F and  H + 
HF(  v ) -*  HF(  v ' ) + H,"  J.  Chen.  Phys . , 58,  3C38,  1973- 

29.  Wilkins,  R.  L. , "Monte  Carlo  Calculations  of  Reaction  Rates  and  Energy 
Distributions  Among  Reaction  Products.  IV.  F + KF(v)  -»  HF(v')  + F 
and  F + DF(v)  -»DF(v')  + F,"  J,  Chem.  Phys.,  698,  1973. 

Cohen,  N.  and  J.  F.  Bott,  "A  Review  of  Rate  Coefficients  in  the  Ho-Fq 
Chemical  Laser  System,"  SAMSO  Report,  TR-76-82,  Aerospace  Coi’poration, 
1976. 
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rates  scaled  according  to  the  measurements  of  Boi*;  nul^iquanta  7-~  transi- 
••ions  individually  scaled  linearly  with  v:  and  V-V  transitions  scaled  with  d . 
The  comparison  with  experiment  for  run  conditions  1 through  U are  shown  in  Figures  90 
through  1C$.  At  first  the  calculations  were  performed  allowing  the  model  to  predict 
■he  HF  development  as  well  as  the  vibrational  distributions.  It  was  decided, 
however,  tha*-  since  the  actual  streamwise  HF  variation  is  known  from  the 
experimental  measurements,  and  since  HF  self  deactivation  is  so  important, 
error  in  fitting  the  vibrational  distribution  could  be  minimised  by  specifying 
the  measured  HF  development.  The  resultant  predicted  vibrational  distribu- 
tions are  also  shown  in  Figures  90  through  109.  Indeed  the  fit  to  the  data  is 
improved.  There  is,  however,  a tendency  in  cases  3 ana  h for  the  model  to 
over  predict  deactivation  in  the  latter  stages  of  the  flow  (Ports  U and  5). 

One  shortcoming  of  the  model  is  its  apparent  inability  +0  predict  the 
vibrational  distribution  at  early  stages  of  the  flow,  in  particular  at  mhe 
first  pert  location  (x  = U cm).  In  fhe  earlier  discussion  of  chemical  reac- 
tion rates  the  possibility  of  the  presence  of  seme  HF  initially  was  presented. 
If  we  were  to  consider  such  a possibility  new,  the  agreement  at  x = k improves 
considerably.  This  is  illustrated  in  Figure  110  for  run  condition  3 with  an 
initial  mole  fraction  of  0.0001S  of  HF  (see  Figs.  k2  and  1+9 )•  This  figure 
shoul  d be  compared  with  Figure  ICO.  The  match  for  vibrational  levels  0,  1, 
and  2 is  vastly  improved.  Similar  results  hold  true  for  cases  1,  2,  and  h. 

At  later  vct^s  *he  vibrational  distributions  are  comparable  with  these 
previously  predicted. 

Thus  we  see  that  if  recirculation  in  the  mixing  region  of  the  flow  t-^be 
dees  result  in  increased  concentrations  of  HF  initially  the  set  of  vibrational 
rates  selected  gives  good  agreement  with  experiment  at  all  measurement 
locations . 

An  interesting  consequence  of  the  changes  to  the  vibrational  rate  package 

is  fcha-  ‘he  rapid  V-T  rates  employed  render  the  comparisons  of  the  model  wi+h 

experiment  relatively  insensitive  to  the  scaling  of  the  single-quantum  V-V 

ra*  ss  in  the  model.  Previously  i‘  had  been  thought  that  the  lo ng- rang  • V-V 

me  iel  employed  in.  CHETfLZ  wcui  1 have  to  be  drastically  changed  to  avcii  -h~ 

severe  anharncr.ic  pumping  and  ~c~al  inversions  on  high  vibrational  levels 

that  had  been  predic*ed  in  the  Task  I studies.  To  avoid  this  possibility  the 

transition  probabilities  were  scaled  differently  with  v,  for  example,  as 
— T fO 

v t' “ "ather  "han  the  typical  linear  scaling.  This  reduced  upper  level 
* , - • 1 . ~ : ? . . . iq  , • . r • - : jurrs  r 
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FIGURE  91  COMPARISON  WITH  MODEL 


NUMBER  DENSITY.  ny  (cm- 


RUN  CONDITION  2 
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FIGURE  95  COMPARISON  WITH  MODEL 
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FIGURE  96  COMPARISON  WITH  MODEL 
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PORT  I - X = 4 cm 


FIGURE  100  COMPARISON  WITH  MODEL 
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RUN  CONDITION  3 
PORT  2 — X * 8 cm 


FIGURE  101  COMPARISON  WITH  MODEL 
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FIGURE  102  COMPARISON  WITH  MODEL 
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(Refs.  31  and  32  ).  In  the  course  of  parameter  variations  to  optimize 

» the  fit  of  the  model  to  the  data,  linear  scaling  with  v of  the  V-V  rates  was 

, also  tried.  A typical  comparison  of  the  model  with  linear  and  v"-'-/  2 scaling 

— is  shewn  in  Figures  m and  312.  The  differences  are  seen  to  be  relatively 

small,  despite  the  large  differences  in  the  upper  level  probabilities.  The 

faster  V-V  package  generally  tends  to  do  a better  job  of  matching  the 

experimental  distributions  for  the  lower  levels  at  short  times,  but  for 

longer  times  it  -ends  to  overestimate  the  ground  state  population  and  those 

above  v = 5>  while  underestimating  the  populations  for  v = 1 to  4 (Figs.  Ill  and 

31?).  '•.'hile  this  gives  generally  poorer  agreement  than,  say,  use  of  a 
-1  7o 

v scaling  for  these  flow  tube  tests,  m power  extraction  calculations  the 
faster  V-V  scaling  tends  to  give  more  pessimistic  (and  perhaps  more  realistic) 
predictions.  These  will  be  described  later  on.  This  is  in  part  because 
of  the  better  agreement  on  the  lower  levels  seen  for  short  times,  comparable 
to  those  in  lasing  systems.  The  poorer  agreement  at  later  times  and  on 
higher  levels,  which  do  not  lase  strongly,  typically  has  only  a secondary 
influence  on  predictions  of  power  extraction. 

To  further  test  the  model  and  isolate  HF  self  deactivation  tests  a com- 
parison was  also  made  for  cases  in  which  HF  was  added  to  the  flow.  In  par- 
ticular 50  cc/sec  and  100  cc/sec  were  added  to  run  conditions  1 through  3.  This  is 
I equivalent  to  initial. mole  fractions  of  0.0002  and  0.0004  of  HF  in  the  flow, 

comparable  with  the  amount  of  HF  produced  by  the  chemical  reactions . A 
comparison  of  the  model  with  experiment  is  shown  in  Figures  113  through  126.  Since 
no  v = 0 data  are  available  to  provide  experimental  values  for  the  streamwise 
development  of  HF,  calculated  values  are  used.  Even  so  the  agreement 
remains  generally  good.  The  tendency  is  for  the  model  predictions  to  fall 
along,  or  somewhat  below  the  experimental  measurements.  This  indicates 
that  HF  self  deactivation  is  perhaps  not  quite  as  strong  as  included  in  the 
model . 


To  better  isolate  the  effects  of  H atoms  two  additional  sets  of  runs 
were  performed  (Conditions  5 and  6).  These  conditions  are  also  listed  in 
Table  9.  The  object  was  to  find  two  cases  with  approximately  the  same 
rate  of  HF  production,  but  having  one  with  a significantly  higher  H-atcm 
concentration  (Condition  6)  than  the  other  (Condition  5).  Thus  differences 
in  ‘he  vibrational  populations  would  be  due  primarily  to  the  differences  in 
the  H-atom  concentration. 


of 


In  order  to  assess 
populations  for  give 


the  influence  of  H atoms,  comparisons  of  the  ratios 
n levels  between  the  *-wc  cases  were  made  at  each  eor‘ . 


31.  Kwok,  M.  private  commui.ication,  August  1975. 

32.  Shin,  H,  K. , "Vibration-Vibration  Energy  Transfer  in  HF  Dimers," 
J.  Chem.  Phys.,  6£,  2091,  1975- 
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FIGURE  117  COMPARISON  WITH  MODEL 
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In  the  ideal  situation  the  streairwise  variation  of  HP  would  be  the  same 
for  each  case.  Then  any  departure  of  the  ratios  from  unity  would  be  pre- 
dominantly the  result  of  H atom  effects.  If  such  effects  were  accurately 
modeled  the  calculations  would  compare  well  even  if  HF  self  deactivation  were 
poorly  modeled. 

Figure  127  shows  the  measured  development  of  HF  for  the  *:wc  cases. 
Unfortunately  differences  in  the  initial  reaction  rates  result  initially  in 
somewhat  larger  HF  concentration  for  Condition  6,  whereas  at  latter  stages 
of  the  flew  Condition  5 has  high  concentrations.  Thus  differences  in  the 
HF  self  relaxation  will  tend  to  mask  the  H atom  effects.  This  problem  is 
further  complicated  by  the  possibility  that  HF  is  initially  present,  as 
illustrated  by  the  dashed  lines  in  the  figure. 

The  experimen*al  population  ratios  were  evaluated  and  compared  with 
three  sets  of  calculations:  (l)  H atom  deactivation  scaled  t o Bott's 

measurements,  (2)  H atom  deactivation  set  at  the  nominal  values,  (3)  K atom 
deactivation  set  to  zero.  The  third  set  of  calculations  was  performed  to 
assess  the  influence  of  HF  self  relaxation  alone.  The  results  of  the  com- 
parison are  shown  in  Tables  12  through  16.  An  earlier  comparison 
using  calculated  net  HF  development,  had  given  fairly  good  agreement 
(Ref.  23).  Now,  however,  with  the  HF  variation  specified  from  measurement, 
the  comparison  is  less  favorable.  The  predicted  influence  of  HF  self  deacti- 
vation  alone  gives  ratios  comparable  to  those  measured  experimentally. 

Hence,  little  can  be  said  about  the  influence  of  the  H atoms. 

A clearer  illustration  of  the  problem  is  available  in  Figures  123  through  1:1. 
These  shew  vibrational  distributions  at  x = 8 and  x = 31  cm,  respectively, 
for  the  two  cases.  As  can  be  seen,  the  computation  assuming  no  H atom 
deactivation  under-predicts  the  upper  level  populations  for  both  cases  at 
the  latter  stages  of  the  flew.  I1  is  apparent  that  either  shor*  comings 
still  exist  in  the  self-deactivation  model  or  errors  in  the  measured  HF 
development  preclude  a valid  assessment  of  H atom  influences. 

“^lls  ed-Fluorescer.ce 

The  data  available  from  the  pulsed-flucrescence  experiments  that  can  be 
used  to  study  vibrational  relaxation  are  primarily  qualitative  in  nature,  in 
the  form  of  fluorescence  traces  from  the  vibrational  levels  v = 1 through  6.  Some 
quantitative  data,  namely  gain  measurements,  are  also  available. 

Due  to  the  short  scales  and  weak  signals  the  fluorescence  traces 

could  not  he  calibrated  to  provide  absolute  number  densities.  Thus  the 
traces  were  normalized  to  unity  at  the  maximum.  One  case  from  each  series 
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TABLE  12.  x = U cm  - PORT  1 


Population  Ratios:  Cond  5/Cond  6 

V 

Exptl. 

Bott  H-Atom 
Deactivation 

Nominal  H- 
Atom  Deactivation 

No  H-Atom 
Deactivation 

B 

1.3 

1.5 

1.3 

1.0 

9 

2.0 

2.4 

2.1 

1.5 

5 

2.9 

2.3 

2.2 

2.0 

6 

2.7 

2.2 

2.0 

1.9 
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TABLE  14.  x = 12.1  cm  - PORT  3 


Population  Ratios:  Cond  5/Cond  6 


Exptl. 

Bott  H-Atom 
Deactivation 

Nominal  H- 
Atom  Deactivation 

No  H-Atom 
Deactivation 

2.0 

3.1 

2.9 

1.6 

3.3 

3.6 

3.1 

2.1 

2.7 

3.3 

2.6 

2.2 

2.1 

3.2 

-3- 

C\J 

2.1 

192 
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table  15.  X = 21.4  ?rn  - PORT  4 


Population  Ratios:  Cond  5/Cond  6 

V 

Exptl. 

Bott  H-Atom 
Deactivation 

Nominal  H- 
Atom  Deactivation 

No  H-Atom 
Deactivation 

■ 

2.1 

4.3 

3-4 

2.2 

s 

— 

4.4 

3.4 

2.5 

5 

— 

4.2 

3.0 

2.5 

6 

— 

4.2 

3.0 

2.5 

L93 


TABLE  1 6.  x = 31  cm  - PORT  5 


Population  Ratios:  Cond  5/Cond  6 


Exptl. 


1.9 


3ott  H-Atom 
Deactivation 


4.1 

4.0 

3.7 


Nominal  H- 
Atom  Deactivation 


3.2 

3.0 

2.7 


No  H-Atom 
Deactivat ion 


2.2 

2.4 

2.2 


2.5 
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that  appeared  representative  of  that  series  was  selected  for  comparison 
with  the  computer  model.  The  cases  selected  and  the  input  conditions  used  in 
the  computer  calculation  are  given  in  Table  17. 

The  calculations  were  first  performed  using  '-he  ra'-e  package  giving  the 
best  fit  in  the  flow  tube  experiments:  (l)  H atom  rates  scaled  according 

'-o  Bot4:,  (2)  Multiquanta  V-T  transitions  individually  scaled  with  v,  and  (3) 
V-V  transitions  scaled  with  v-J-'-.  The  comparison  of  these  cal  cola4:  ions 
with  experiment  is  shown  in  Figures  132  through  lU8.  Generally  the  slopes  in  the 
diminishing  part  of  the  curves  agree  fairly  well  between  model  and  experi- 
ment, but  the  model  predictions  for  the  position  of  the  neak  lag  far  behind 
the  experimentally  observed  peaks.  This  more  rapid  turn  around  is  likely  due 
to  the  same  fast  initial  adjustment  of  the  populations  seen  in  the  flow  tube 
experiments  that  was  also  poorly  simulated  by  the  model. 

Additional  calculations  were  made  using  the  faster,  linear-v  scaling 
on  the  V-V  rates.  A comparison  of  these  (FAST  V-V)  with  the  slower  (SLOW  V-V) 
rates  and  experiments  are  shown  in  Figures  143  through  148.  The  results  are  typical 
for  all  the  other  cases.  The  peak  for  v = 1 and  v s 6 tends  to  occur  later 
than  the  "slow  W"  predictions,  and  that  for  v = 2 through  5 tends  to  occur  earlier. 
Also  the  negative  slopes  tend  to  be  shallower. 

Since  the  above  comparisons  are  only  of  relative  quantities  it  is 
difficult  to  draw  any  firm  conclusions.  There  is  no  way  of  knowing  wha4-.  con- 
centration the  experimental  maximum  represents  and  how  it  compares  with  the 
model.  It  does,  however,  appear  that  population  redistribution  as  trea'-.ed 
in  the  model  is  no4-  fast  enough  to  produce  the  behavior  observed  in  the 
experiments . 

The  measurements  of  gain  made  experimentally  provide  a more  quantita- 
tive means  of  gaining  perspective  on  the  adequacy  of  the  relaxation  model. 

An  experimental  study  of  the  maximum  gain  as  a function  of  transition  was 
made  for  conditions  corresponding  to  series  F (0.2  torr  Hp,  1.0  torr  ?0). 

A comparison  with  the  model  is  shown  in  Table  18.  Two  calculations  ,are 
shown  — slow  W,  which  is  the  "best  fit"  rate  package  including  v--1"'  - 
scaling  of  the  V-V  rates;  fast  W which  is  fhe  same  se4  of  rates  except 
linear  v scaling  is  employed  with  the  V-V  rates. 

We  see  tha4-  on  the  v = 1 -»  0 transitions  the  predicted  maximum  gains 
for  both  the  fast  and  slow  W calculations  are  lower  than  measured  experi- 
mentally, with  4he  slow  V-V  calculation  being  in  somewhat  better  agreement. 

Cn  4he  other  hand,  for  v = 2 -»  1 transitions  the  predicted  maxima  are  higher 
t.han  4he  measured  values,  with  the  fast  W calculations  being  in  closer 
agreement . 
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TABLE  17.  INITIAL  CONDITIONS  - FJIBED  FLUORESCENCE 
VIBRATIONAL  STUDIES 


Initial  Pressure  = 100  torr 
Initial  Temperature  = 300°K 

Mole  Fractions 


Case 

F 

f2 

«2 

H "■*' 

Ar 

A8 

1.02x10-5 

0.008477 

0.008677 

0.001863 

0.981003 

33 

7.20x10-° 

0.009130 

0.021180 

0.002174 

0.967509 

C3 

7.70 x 10'° 

0.009707 

0.005336 

0 

0 . 984668 

DH7 

5.47x10-° 

0.004907 

0.005497 

0 

0.969591 

EF6 

7.3  x 10‘° 

0.004938 

0.010155 

0 

0.964900 

F7* 

1.35x10-5 

0.01009 

0.002345 

0 

0.987539 

*The  traces  labeled  G i 

n the  figures 

are  for  essen 

tially  the  same 

conditions 

as  F7. 
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TABLE  18.  COMPARISON  OF  PREDICTED  AND  MEASURED  GAIN 


Line  Variation 


Mixture : 

0.2  Hg,  : 

L.O  ?2,  98.8  Ar 

Line 

Experiment 

Model 

Pulse  Energy 
(J) 

Gain 

$/cm 

Pulse  Energy 

(J) 

Gain 
Slow  W 

Gain 
Fast  W 

1P2 

.14 

.98 

.13 

.66 

.45 

1P3 

.13 

.78 

• 13 

• 89 

.67 

1P4 

.12 

1.01 

.13 

.86 

.63 

1P6 

.12-.  14 

.55 

.13 

• 36 

.22 

1P7 

.11 

.26 

• 13 

.16 

.091 

2P3 

.13 

1.2 

• 13 

2.59 

1.3 

2P4 

.12 

.66 

• 13 

2.36 

1.15 

2P5 

.11 

.48 

• 13 

1.63 

.70 

2P6 

.12-. 14 

.49 

.13 

.88 

.34 
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RELATIVE  POPULATION  DEVELOPMENT 


FIGURE  13?  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
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RELATIVE  POPULATION  DEVELOPMENT 
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PULSE  CASE  C3 
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FIGURE  140  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
RELATIVE  POPULATION  DEVELOPMENT 


PULSE  CASE  DH7 
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FIGURE  141  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
RELATIVE  POPULATION  DEVELOPMENT 
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FIGURE  142  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
RELATIVE  POPULATION  DEVELOPMENT 
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FIGURE  143  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
RELATIVE  POPULATION  DEVELOPMENT 


PULSE  CASE  F7 


FIGURE  145  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
RELATIVE  POPULATION  DEVELOPMENT 
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RELATIVE  POPULATION  DEVELOPMENT 
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FIGURE  147  COMPARISON  OF  MODEL  WITH  EXPERIMENT 
RELATIVE  POPULATION  DEVELOPMENT 


It.  should  also  be  noted  that  in  the  experiment  the  gain  was  observed 
to  have  a duration  of  from  20  to  40  usee.  This  is  considerably  less  than  pre- 
dicted by  the  model.  The  slow  W cases  had  durations  ranging  from  100  to 
300  s,sec  and  the  fast  W cases  from  90  to  200  ^sec.  These  facts  indicate 
that  the  model  is  inadequate  in  describing  the  rapid  redistribution  of  the 
primitive  distribution  on  short  time  scales.  The  strongly  pumped  v = 2 popu- 
lation remains  too  high  *:oo  long,  and  the  v = 1 population  remains  too  low 
(relatively) . 


The  strong  influence  of  initial  HF  contamination  must  be  nc'-ed.  In  some 
of  the  experiments  there  was  an  apparent  prereaction  resulting  in  small  con- 
centrations of  HF  at  the  tine  of  the  pulse  initiation.  To  assess  the  possible 
influence  of  such  HF,  fwo  comparable  cases  were  run  for  run  condition  C3 
(0.5  torr  R-,,  1.0  tcrr  Fo).  In  one  case  no  HF  was  present  at  initiation,  and 
in  the  other  0.03  tcrr  was  present.  This  represents  3 lo  of  the  maximum 
possible  HF.  The  results  are  summarized  in  Table  19>  We  see  that  the 
effect  is  quite  dramatic.  Far  smaller  concentrations  of  initial  HF  could 
strongly  influence  measured  gains  and  render  comparisons  very  difficult. 

A final  comparison  was  made  to  see  if  the  model  could  properly  predict 
trends  as  the  mixture  of  H-,  and  ?2  wa3  varied.  The  comparison  is  shown  in 
Table  20.  The  model  does  quite  well  in  this  respect,  correctly  predicting 

which  mixture  would  give  the  highest  gain.  The  only  large  discrepancy  is 
for  the  hydrogen  rich  case  (l.O  torr  rL,,  0.5  torr  FQ),  but  that  might  be  the 
result  of  some  HF  contaminant  in  the  experiment. 

.All  ir.  all  it  appears  that  the  relaxation  model  does  an  adequate  Job  of 
nredicting  vibrational  distributions  and  net  vibrational  energy  less  in  most 
cases.  However,  over  short  time  scales  the  model  does  not  duplicate  the 
rapid  redistribution  of  vibrational  oopulations  that  i=  acparerh  in  experi- 
ment ( flow  tube  and  pulsed  fluorescence).  This  indicates  that  a be'-ter  VY 
model,  including  multiple-quanta  transitions  is  needed  to  provide  the  rapid 
redistribution  while  avoiding  the  problems  of  anhamonic  pumping  evident  in 
single  quantum  models. 

Furthermore,  over  longer  time  scales  the  net  less  of  vibrational  energy 
is  occasionally  over  oredicted.  This  implies  that  a somewhat  slower  V-T  ra*  e 
model  might  appropriate  along  with  revisions  to  the  V1/  model. 


Performance  Calculations 

While  the  rate  nackages  employed  in  this  study  may  not  v,e  <-he  last  word 
with  regard  to  HF  chain  kinetics,  iv  is  of  interest  tc  see  what  influence  they 
may  have  on  predictions  of  laser  performance. 
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TABLE  19.  PREDICTED  GAIN 


Influence  of  Initial  HE 
Mixture:  O.5H2  1.0  Hg  98.5  A r 


Gain  (fo/cm) 


Line 

No  Initial  HF 

0.3  Torr  HE 

1P2 

1.10 

Absorption 

1P3 

1.25 

1P4 

1.07 

1P5 

0.71 

1P6 

0.37 

1P7 

0.15 

! 

2P3 

3.70 

0.0027 

2P4 

3.10 

0.0035 

2P5 

2.00 

0.0041 

2P6 

1.02 

0.0047 
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TABLE  20.  COMPARISON  OF  PREDICTED  AND  MEASURED  GAIN 


Mixture  Variation 
Line  1PU 


Experiment  Model 


Mixture 

Pulse  En. 

Gain 

Pulse  En. 

Gain 

(J) 

%/cm 

(J) 

%/cm 

0.2  Hg  1.0  F2 

.12 

1.01 

.13 

0.86 

0.5  % 0.5  F2 

.11 

0.51* 

.11 

O.U5 

0.5  Hg  1.0  F2 

.11 

1.40 

• 09 

1.10 

1.0  Hg  0.5  F2 

.08 

0.09 

.10 

.62 
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Two  se^s  of  rates  were  used  for  “-he  calcula4" ions , designated  as  C and 
These  are  rates  designed  to  give  the  best  overall  fits  bo*h  ‘he  flow  tube 
and  pulsed  fluorescence  experiments.  The  set  designated  as  C nay  be  s'unr.a- 
rized  as  follows: 

(1)  Hot  reaction  rate  = 0.9  x nominal 

(2)  Cold  reaction  rate  = 1.5  x nominal 

(3)  Multiquanta  V-T  deactivation  scaled  linearly  with  v 

(U)  Single  quantum  Sharma-model  W exchange  scaled  with  v”  ' 

(5)  H atom  deactivation  scaled  according  to  the  measurements  of  Bott. 
The  set  designated  C7  is  identical  except  the  V-V  scaling  is  linear  with  v. 

I4’  was  first  of  interest  to  assess  the  influence  of  the  new  rate  package 
on  predicted  performance  of  cold  reaction  lasers.  Sample  conditions  of 
calculations  for  the  long-range  cold-reaction  laser  (matrix  hole  nozzle  LRCl) 
are  given  in  Table  21.  The  resulting  calculations  are  summarized  in 
Table  22.  We  see  tha~  generally  little  change  in  *:he  predicted  laser  out- 
put results  from  use  of  the  new  rates.  The  largest  variation  is  between  the 
two  V-V  models  for  the  5 torr  case.  This  is  an  indication  tha*:  the  V-V  model 
selected  can  strongly  influence  power  extraction  predic'ions,  even  if  there 
is  little  apparent  difference  in  unloaded  population  distributions. 

A comparison  was  also  made  with  some  calculations  previously  performed 
under  a U.S.  Army  sponsored  chain-reaction  demonstration  contract.  (Ref.  33)* 
Two  cases  in  particular  were  considered:  (l)  The  initial  design  point,  and 

^ 2 ) a later,  mere  realistic  design  poinT.  The  input  conii-icr.s  for  these 
are  given  in  Tables  23  and  2k.  The  resulting  predictions  are  shown  in  Tables 
25  and  26.  These  calculations  are  admittedly  optimistic  since  aerodynamic 
and  optical  effects  likely  to  degrade  performance  are  not  adequately  included 
in  the  model.  Even  so  the  predicted  performance  is  in  the  range  which  render 
the  chain  reaction  at  test  marginally  competitive  with  the  cold  reaction  as 
a laser  pumping  mechanism. 


00 . 


Kepler,  C.  E.,  et  al.,  "Demonstration  of  HF  Chain  Reaction 
iJTF.C  Technical  Report  RK-CR-75-^,  UTRC,  East  Hartford,  CT, 


Laser, " 
July  1974. 
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TABLE  21.  INPUT  CONDITIONS  - LONG  RANGE  COLD  REACTION  LASER 


O'  = 0.90 

He  = 0.157 

p = 1.0,  5.0 

gQ  = 0.003  cm"1 

L mix  = °-77 

Constant  Pressure  Option 


Cavity  Entrance  Conditions 


X 


He 


*DF 


u (cm/sec) 


T ( °K) 


Primary 

0.15652 

0.1U311 

0.00626 

0.11271 

IIU3U5 

268 


Secondary 


0.58114.0 

36481+5 

339 
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TABLE  22.  PREDICTED  PERFORMANCE  - COLD  REACTION  LASER 


Case 

Rate  Package 

X -cm 
Li 

a-J/g  ( kJ/lb ) 

P * 1 

Nominal 

0.79 

723 

( 329) 

C 

0.77 

71U 

( 325  ) 

CF 

0.71 

688 

(313) 

P * 5 

Nominal 

0.66 

271 

(123) 

C 

1.18 

381 

(173  ) 

CF 

0.78 

275 

(125) 
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TABLE  23.  IUPTJT  CONDITIONS  - DESIGN  POINT 
CHAIN  REACTION  LASER 


Qf  = 0.05 
He  = 0.85 
p = 20  torr 


g0  = 0.001  cm' 


Lmi x = 5 cm 


Constant  Pressure  Option 


Cavity  Entrance  Conditions 


X- 


He 


72 


XH2 

u( cm/ sec ) 
T ( °K  ) 


Primary 

.63512 

.07098 

.OOJbJ 

212000. 

150 


226 


Secondary 

.21173 

.07^70 
212000 . 
75 


• ■ *• 


TABLE  2k.  INPUT  CONDITIONS  - REPRESENTATIVE 
CHAIN  REACTION  LASER 


or  = 0.10 
He  =0.525 


p = 20  torr 

go  = 0.001 


W = 10  cm 


Constant  Pressure  Option 


Cavity  Entrance  Conditions 


Primary 

Secondary 

xHe 

•69375 

.23125 

Xp2 

.03068 

— 

xF 

.00682 

— 

XH2 

— 

.03750 

u (cm/sec) 

226600 

195770 

T (°K) 

157 

75 

227 


TABLE  25. 

PREDICTED  PERFORMANCE 
CHAIN  REACTION  LASER 

- DESIGN  POINT 

Rate 

Package 

xL- 

cm 

c- 

J/g, (kj/lb) 

v_>_2 

a 

Nominal 

22 

1012  (46o) 

0.39 

C 

20 

572  (260) 

0.21 

CF 

20 

480  (218) 

0.28 

TABLE  2 6.  PREDICTED  PERFORMANCE  _ REPRESENTATIVE 
CHAIN  REACTION  LASER 
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To  see  whether  a more  competitive  set  of  conditions  cculi  be  found  for 
the  chain  reaction,  some  conditions  were  used  in  calculations  ever  a range 
of  higher  pressures  (+c  aid  in  recovery).  The  input  conditions  for  these 
are  shown  in  Table  27.  Included  are  some  fluorine  rich  cases  reflecting  the 
experimental  observation  tha*~  they  produce  the  best  gain.  The  results  of 
these  calculations  are  summarized  in  Tables  28and29.  We  see  that, 
on  paper  af  least,  higher  pressure  performance,  especially  fluorine  rich, 
provides  more  favorable  performance  from  the  chain,  indicating  a possible 
regime  in  which  such  lasers  might  have  seme  specialized  applications. 

Discussion 

Chemical  Reactions 

Some  difficulty  was  encountered  in  trying  to  find  a set  of  chemical 
reaction  rates  that  would  fit  both  the  flow-tube  and  the  pulsed  flourescence 
data.  The  flow-tube  results  implied  somewhat  lower  rates  than  did  the 
pulsed-flourescence  data.  While  some  systematic  error  or  unaccounted  for 
reaction  mechanism  might  explain  the  differences,  we  were  unable  to  find 
a satisfactory  resolution  to  the  discrepancies.  Consequently  a fairly  wide 
range  must  be  applied  to  the  values  for  the  rate  coefficients  we  have 
evaluated.  For  the  hot  reaction  we  find  over  the  range  300  < T <t  380°K 

k^  = 1.8  x 10' 10  exp  (-2400/RT)  cc/molecule-sec 

with  an  error  band  of  ± 40  percent.  This  is  90  percent  of  the  nominal  value 
of  Albright, et  al.  Because  of  the  limited  temperature  range,  no  conclusion 
could  be  drawn  about  the  temperature  dependence  of  the  rate,  hence  Albright's 
value  has  been  used.  There  is  a wide  enough  range  in  temperatures,  however, 
to  require  the  inclusion  of  temperature  dependence  in  order  to  achieve  a 
realistic  fit  to  the  data. 

For  the  cold  reaction  we  find  that 

^ = 4.0  x 10  exp  ( -1600/RT ) cc/molecule-sec 

with  an  error  band  of  ± 35  percent.  The  nominal  value  of  Horamann,et  al., 
lies  at  the  lower  end  of  this  range.  Again,  no  conclusions  could  be  drawn 
with  regard  to  the  temperature  dependence. 

In  the  pulse  experiments,  very  fast  back  reaction  rate  coefficients 
that  would  be  most  important  in  the  higher  concentration  cases  (A,  B,  and  C) 
and  those  with  stronger  upper  level  populations  (C  and  F)  can  be  invoked  to 


* 
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TABLE  27.  INPUT  CONDITIONS  - CHAIN  REACTION 
PRESSURE  VARIATION  STUDIES 


a = O.05 

p =20,  60  torr 

go  = 0.001 

Lmix  = 10,  30  cm 

Constant  Pressure  Option 
CAVITY  ENTRANCE  CONDITIONS 


Case  - ?2.5i  He,  F2/H2  = 1 


Primary 

Secondary 

xHe 

.69375 

.23125 

XF2 

.03393 

— 

xF 

.00357 

— 

OJ 

& 

— 

0.03750 

u 

226600 

195770 

T 

157 

75 

Case  - 95 % He,  F2/H2  = 1 

xHe 

.7125 

.2375 

xF2 

.02262 

— 

xF 

.00238 

— 

xH2 

— 

.002500 

u 

226600 

195770 

T 

157 

75 

Case  - 88.75 i He, 

F2/H2  = 2 (a  = 0.025) 

xHe 

.55625 

* .23125 

Xp2 

.071^3 

— 

Xp 

.00357 

— 

%2 

— 

0.0375 

u 

226600 

195770 

T 

157 

75 

TABLE  28. 


PREDICTED  PERFORMANCE  - 
EFFECT  OF  PRESSURE  VARIATION 


Rate 

P 

XL 

a 

av>3 

P 

Q 

Te 

Case 

Package 

(torr ) 

(cm) 

(j/g) 

C7 

(torr) 

fHF 

(°K) 

M 

92.5$  He 

C 

20 

56 

229 

0.20 

74 

0.88 

808 

1.48 

60 

18 

286 

0.21 

220 

0.89 

872 

1.43 

He 

20 

13 

14 

0 

304 

0.05 

165 

3.15 

60 

48 

128 

0.32 

250 

0.93 

654 

1.58 

92.5$  He 

Cp 

20 

51 

136 

0.27 

95 

0.62 

614 

1.70 

60 

17 

202 

0.29 

292 

0.61 

693 

I.60 

He 

20 

20 

13 

0 

306 

0.05 

165 

3-15 

60 

39 

77 

0.46 

289 

0.75 

561 

1.60 

TABLE  29.  PREDICTED  PERFORMANCE  - 
EFFECT  OF  FLUORINE 


Rate 

D 

xL 

c 

CTv  >3 

po 

oe 

Case 

Package 

(torr) 

(cm) 

(J/g) 

c 

(torr) 

fHF 

(°K) 

M 

F2/K2=l 

Cp 

20 

51 

136 

0.27 

95 

0.62 

6l4 

1.70 

60 

17 

202 

0.29 

292 

0.61 

693 

1.60 

F2/H2=2 

20 

33 

341 

0.36 

96 

0.85 

783 

1.70 

60 

13 

321 

0.38 

337 

0.69 

832 

1.65 

231 
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explain  the  anomalously  fast  reaction  rates  seen  for  the  lower  concentration 
cases  (D  and  E).  But  that  would  require  correspondingly  faster  forward  rates. 
Furthermore  the  back  reaction  rates  now  are  already  taken  to  be  nearly  gas 
kinetic  for  the  higher  vibrational  levels.  While  there  is  not  enough  firm 
evidence  here  to  justify  support  of  such  a possibility,  it  is  perhaps  some- 
thing that  bears  further  study. 

Supplemental  chemical  reactions  were  not  found  to  be  important.  Hot 
atom  effects  could  be  important  in  specialized  situations  used  to  create 
F atoms  and  might  come  into  play  in  pulsed  devices.  Such  a situation  could 
best  be  handled  in  a model  by  appropriately  modifying  the  initial  F atom 
concentration  to  reflect  additional  dissociation  from  F*-F2  collisions. 

Any  hot  atom  effects  arising  from  the  chain  reaction  are  likely  to  be  masked 
in  the  overall  hot  and  cold  reaction  rates.  In  the  highly  dilute  cases 
tested  here,  such  effects  would  be  minimal.  Less  dilute  situations  might 
result  in  slightly  faster  effective  rates  (Ref.  3*0  • 

Wall’ interactions , which  are  highly  specific  to  individual  devices, 
can  in  some  cases  be  important.  At  the  very  least  in  these  experiments 
H-atom  loss  appears  to  have  played  a role  in  the  flow  tube  and  possibly 
also  in  the  pulse  cases  (to  help  explain  the  lessening  rate  of  reaction 
at  later  times).  In  any  case  care  should  be  taken  in  evaluating  any  device 
to  assess  the  influence  of  the  walls  if  possible. 

The  currently  recommended  chemical  reaction  rates  for  H2-F2  chain 
reaction  lasers  are  listed  in  Table  3 0. 

Vibrational  Deactivation 


Since  there  are  so  many  potential  transitions  involving  HF  vibrational 
states  and  since  there  is  a fair  amount  of  error  in  any  experimental  measure- 
ments, no  set  of  rates  that  might  be  proposed  is  necessarily  unique.  Widely 
divergent  approaches  as  to  which  transitions  are  important  can  lead  to  equally 
good  fits  of  the  available  data. 

In  this  program  we  have  taken  approach  that  assumes  V-V  interactions 
to  be  relatively  unimportant  for  upper  vibrational  levels,  while  V-T  deac- 
tivation scales  up  rapidly  with  increasing  vibrational  level.  For  such  a 
model  we  find  relatively  good  agreement  with  the  experimental  data.  Some 
of  the  salient  features  of  the  model  are  summarized  below. 


k.  Sullivan,  J.  H.  and  R.  C.  Feter,  "The  kinetics  of  an  ^-Fo  Explosion  as 
Initiated  by  a Pulse  of  Fluorine  Atoms,"  Report  LA-5‘+09,  Los  Alamos 
Scientific  Laboratory,  Los  Alamos , X!,  1 '7*-. 


TABLE  30.  RECOMMENDED  CHEMICAL  RATES 


l’vkf 

1.  H + ?2  - HF(v)  + F AH  = -98  kcal/mole 

kj-  = 1.8 x 10"^ exp  (-2400/RT)  cc/molecule  sec 

k^,  =2.96  x 10  exp  _-(l0040C  - EV)/RT]  cc/molecule  sec 

(fv,  v = 0,  11)  = .012,  .027,  .033,  .039,  .134,  .265,  .298, 
.134,  .058,  0,  0,  0 


fv*f 

2.  F + H2  -»  HF(v)  + H AH  = -31.6  kcal/mole 

-1C 

k+,  =4.0  x 10  exp  (-I6OO/RT)  cc/molecule  sec 
-10 

k^  =4.4  x 10  exp[-(33200-Ey) /RT]  cc/mclecule  sec 

( fv,  v = 0,  11)  = 0,  .1639,  -5574,  .2787,  0,  0,  0,  0,  0,  C,  0,  0 

Optional 

3.  H + F + M -*  HF(O)  + M AH  = -129.6  kcal/mole 

kj.  = 6.6  x ?0  "°  T cc^ /molecule^ sec 

-6-1.0 

= 7*73  x 10  T exp  (-I296OO/RT)  cc/molecule  sec 


5. 


F + F + M - F,  + M AH  = -37-7  kcal/mole 

, -33  -0.5  2 . 2 

k,  = 4.99  x 10  T cc  /molecule  sec 

k.  = -j.3-  ::  10  ^ exp  ( -35100/KT)  :c/molecule  sec 


+ H+  M->h2+M  AH  =-104.2  kcal/mole 

= 2.75  X 10  T cc_/molecule  sec 
^ = e.96  x 10  ° T 1'°"5  exp  (-104200/RT)  cc/molecule  sec 
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TABLE  30.  RECOMMENDED  CHEMICAL  RATES  (COIJT'D) 


c.  HF( v > 3)  + F-  - HF(O)  + 2F  AH  = 37.7  kcal/mole 

“17 

= 3*5  x 10  cc/ molecule  sec 

kb  s ° 

7.  Ho(v  k l)  + ^(°)  + H + F AH  = 5.92  kcal/mole 

-19 

kf  = 8.0  x 10  cc/molecule  sec 

kb  s o 


8.  F2  + F*  -*  3F  AH  = -7  kcal/mole 

kf  = .7  x 10  ' cc/molecule  sec 
kb  = 0 

9.  F*  + Ar  -»  F + Ar  AH=  -44.7  kcal/mole 

kp  = 3 x 10  ^ cc/molecule  sec 
kb  = 0 


1C.  F*  + He  -»  F + He 


AH  = -44.7  kcal/mole 


kp.  = 1 x 10 

X 


-10 


cc/molecule 


sec 


kb  = 0 

Note:  reactions  3-10  have  been  used  as  an  upper  limit  on  hot  F atom  effects 

in  low  diluent  cases  when  F atoms  are  generated  by  a doubled  ruby  laser 
pulse.  It  has  been  assumed  that  in  each  dissociation  one  atom  gets  all  the 
translational  energy  (i.e.,  half  the  atoms  are  hot  - T =*■  15000°K).  It  is 
assumed  that  F*  + ?2  collisions  have  a 50%  probability  of  dissociating  the 
?2m  '^ie  r*tes  for  thermalization  of  F*  assume  it  requires  5 collisions  with 
Ar  or  10  with  He  to  reduce  the  F atom  translational  energy  to  a value  (~1500°K) 
at  which  enhanced  dissociation  would  be  minimal.  In  any  situation  in  which 
the  initiation  of  a pulsed  laser  might  produce  hot  F atoms,  a scheme  of  this 
sort  Cor  more  sophisticated  miwht  be  required  to  adequately  model  -.he  results. 


V-T  deactivation  of  HF  is  treated  as  a multiquanta  process.  All  tran- 
sitions from  a given  level  to  lower  levels  are  taken  to  have  equal  probability. 
Each  is  scaled  linearly  with  vibrational  quantum  number.  Upper  level  rate 
coefficients  can  thus  be  evaluated  by  scaling  the  best  literature  values 
for  the  1-0  transition.  An  alternative  approach  that  also  gives  good  general 
agreement  with  these  experiments  is  to  treat  the  V-T  deactivation  as  single 

p O 

quantum  and  scale  upper  level  rate  coefficients  by  v J.  Such  an  approach 
would  be  easier  to  implement  in  existing  models. 

With  such  fast  V-T  deactivation,  the  influence  of  single -quantum  V-V 
exchange  is  not  as  strong  as  originally  thought.  Here  we  have  modified  our 
existing  Sharma  formulation  by  scaling  with  v-1/2  rather  than  v.  An  equally 
good  alternative  would  be  to  use  a constant  probability  for  all  exothermic 
transitions.  The  probabilities  for  endothermic  transitions  would  be  found 
by  applying  detail  balance.  The  appropriate  rate  coefficients  would  be 
obtained  from  measurements  for  lower  level  transitions  (Refs.  35  and  36). 

Direct  assessment  of  H-atom  deactivation  proved  difficult  as  the  tests 
were  too  strongly  influenced  by  differences  in  HF  self  deactivation.  In 
the  absence  of  convincing  contrary  evidence,  H-atom  rates  scaled  to  Bott's 
measurements  are  being  employed. 

A summary  of  rates  that  gave  the  best  overall  fit  to  the  vibrational 
population  data  within  the  constraints  of  the  model  used  (in  particular  the 
limitation  to  single -quantum  V-V  exchange)  is  given  in  Table  31. 

While  the  above  described  model  gave  fairly  good  general  agreement 
with  experiment,  certain  shortcomings  were  evident.  First  the  model  does 
not  appear  to  adequately  describe  the  rapid  readjustment  of  primitive 
vibrational  distributions  to  nearly  Boltzmann  distributions.  Also,  a "knee" 
which  is  not  evident  in  experiment,  appears  in  the  predicted  vibrational 
distributions  for  v=5  and  v=6  that  tends  to  persist  regardless  of  hew  the 
various  rates  were  modified.  Finally  there  was  the  tendency  in  some  flow- 
tube  cases  (3,4,5  and  6)  to  underpredict  the  populations  of  intermediate 
levels  (v  = 2-4)  during  the  latter  stages  of  the  flow. 


35.  Airey,  J.  R.  and  I.  W.  M.  Smith,  "Quenching  of  Infrared  Chemiluminescence 
Rates  of  Energy  Transfer  from  HF  (v  £ 5)  to  CO2  and  HF,  and  from  DF  (v  £ S'1 
to  COg  and  HF,"  J . Them , Phyc . , 57,  1669,  1972. 

36.  Osgood,  R.  M. , et  al.,  "Measurements  of  V-V  Exchange  Fates  for  Excited 
Vibrational  Levels  (2  £ v £ 4)  in  Hydrogen  Fluoride  las,"  J.  :hem.  Fhyo . , 
60,  1464,  1974. 
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TABLE  31.  RECOMMENDED  CHANGES  TO  ORIGINAL 

VIBRATIONAL  RELAXATION  RATE  PACKAGE 


I.  V-T  Multiple-Quanta  Transitions;  v -»  v-n,  n arbitrary 

A.  For  H and  F Atoms 

k , = aT  b exp(-c/RT)  sec-^  atm-"'' 
v -*  v’ 

(Note  values  for  a in  following  tables  are  multiplied  by  10”^) 


H atoms 

F atoms 

V' 

a 

cr 

c 

a 

b 

e 

B 

0 

.06 

1.275 

1172 

29.87 

1.766 

3099 

v=2 

0 

.045 

1.084 

896 

0 

1 

1 

1 

.0046 

1.201 

722 

27.26 

1.708 

3077 

■ 

0 

1.88 

1.063 

761 

0 

HI 

1 

1 

1.36 

1.078 

675 

27.89 

3718 

H 

2 

1 80 

1.09 

889 

43.2 

IHil 

3645 

0 

7.16 

1.324 

806 

0 

1 

1 

v=4 

1 

.321 

• 91 

635 

30 

1.75 

37C0 

2 

.606 

.966 

667 

66 

1.7 

3600 

3 

.786 

1.001 

519 

100 

1.7 

3600 

0 

4.63 

1.234 

636 

0 

1 

1 

1 

I.92 

1.149 

610 

30 

1.75 

3700 

v=5 

2 

.967 

1.043 

711 

33 

1.7 

3700 

*3 

1-53 

1.058 

668 

100 

1.7 

370C 

4 

1.06 

1.047 

527 

165 

1.7 

37CC 

3.69 

1.189 

620 

0 

1 

1 

1.92 

• 837 

568 

30 

1.75 

3700 

1.97 

1.155 

751 

34 

1.658 

3810 

.821 

1.01 

671 

62 .4 

1.699 

3589 

3-20 

1.211 

596 

160 

1.741 

3638 

5 

4 .41 

1.222 

626 

226 

1.712 

3o44 

TABLE  31.  RECOMMENDED  CHANGES  TO  ORIGINAL 

VIBRATIONAL  RELAXATION  RATE  PACKAGE  (CONT'D) 


^v,v'-n  “ vic10 

k^  as  given  in  Appendix  Cl 

2 # ^ 

(alternative  - single  quantum  with  k , = v~‘Jk,/v) 

v,v-l  10 
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TABLE  31.  RECOMMENDED  CHANGES  TO  ORIGINAL 

VIBRATIONAL  RELAXATION  RATE  PACKAGE  (CONT'D) 


II.  V-V  Transitions 


_/ — L_\  7 3+1  \ 2 rr- 

',r-l  \l-er/  \l-(s+l  )*)  N 

where  K =>  0.8  x 10-11+ 


S 5 + 1 

(Alternative  - k ’ = 4.4K  for  all  exothermic 

v r».  J_ 

reactions ) 


A 


*This  value  Is  selected  such  that  the  transition  (0,2)  -*  (l,l)  has  the 
value  recommended  by  Cohen  (Ref.  30)'. 
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Such  problems  illustrate  the  need  to  further  improve  the  model.  One 
such  approach  would  be  to  employ  arbitrary  scaling  of  the  V-T  rates  (e.g., 
a bit  slower  for  v = 2 through  4 and  faster  for  v = 5 and  6),  but  that  would  be 
tedious  with  little  chance  of  general  application.  Arbitrary  changes  in 
the  single -quantum  V-V  rates  tend  to  introduce  new  problems  while  alle- 
viating old  ones.  For  example,  faster  V-V  rates,  which  would  help  provide 
rapid  early  adjustment  of  primitive  distributions,  tend  to  depress  inter- 
mediate level  populations  and  enhance  upper  level  populations  too  much. 

It  appears,  in  retrospect,  that  a new  mechanism -- namely  multiquanta 
V-V  exchange  — will  be  required  to  improve  theoretical  predictions  of 
vibrational  distributions.  An  approach  assuming  multiquanta  VV  exchange 
in  conjunction  with  less  strongly  scaled  VT  interactions  has  been  proposed 
recently  by  Clendening, et  al.,  (Ref.  37).  This  approach,  based  upon 
surprisal  analysis,  provides  a rational  means  of  evaluating  all  applicable 
transition  rates.  A model  such  as  this  should  be  applied  to  these  experi- 
mental results  in  the  near  future. 

Chain  Reaction  Laser  Performance 


Performance  calculations  with  the  new  rate  package  do  not  indicate  an 
optimistic  future  for  the  HF  chain  reaction  laser,  in  basic  agreement  with 
the  findings  of  Clendening  (Ref.  38)*  However,  it  appears  that  a high- 
pressure  fluorine-rich  operating  condition  can  provide  reasonable  performance 
that  might  be  useful  in  some  applications.  Such  a possibility  should  be 
examined  in  more  detail . 


37.  Clendening,  C.,  et  al.,  "information  Theory  Analysis  of  reactivation 
Rates  in  Chemical  Lasers,"  AFWL-TR-76-lhU,  Kirtland  AF3,  NCC,  1976. 

33.  Clendening,  C.,  "Current  Status  of  HF/DF  Chain  Chemical  Laser," 
AFWL-TR-76-194,  Kirtland  AF3,  NIC,  1976. 
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SECTION  V.  ESR  TITRATION  EXPERIMENT 


Summary 

It  is  shown  that  titration  of  F,  Fg  with  HC1  is  described  principally 
by  the  reaction,  F + HC1  ->  HF  + Cl.  Accordingly,  F concentration  may  be 
determined  by  titration  of  F with  a known  quantity  of  HC1  and  a measurement 
of  HF  concentration.  The  required  HF  concentration  may  be  determined  by 
laser  optical  absorption  methods  (Refs.  39  > l) • The  F concentration  is 
relevant  to  understanding  the  governing  reaction  sequences  in  HF/DF  chain 
reaction  lasers  (Ref.  4o).  The  conclusion  concerning  the  above  reaction  is 
arrived  at  by  observing  the  effect  of  HC1  titration  on  the  intensity  of  an 
F ESR  (Electron  Spin  Resonance)  signal,  the  intensity  being  proportional  to 
concentration.  It  is  found  that  the  intensity  decreases  linearly  with 
increasing  HC1  flow  and  that  1.07  - 0.12  HC1  molecules  are  required  to 
remove  1 F atom.  Since  the  flow  contains  F^  as  well  as  F,  these  latter 
results  demonstrate  the  lesser  importance  of  competing  reactions,  like  Cl  + 
Fg  -*  C1F  + F,  to  the  titration  method  of  F concentration  determination. 

With  regard  to  the  latter  reaction,  present  results  verify  that  the  rate 
constant  associated  with  it  is  small,  which  is  a conclusion  reached  by 
others  earlier  using  a method  different  than  the  present  one  (Ref.  hi). 

Introduction 

In  order  to  verify  that  HC1  titration  of  F is  described  by  F + HC1  -* 
HF  + Cl,  without  significant  interference  from  competing  reactions,  it  is 
required  that  the  relative  amounts  of  F and  HC1  are  known  at  the  titration 
region.  Since  the  species  move  in  a flow  tube,  the  flow  rates  in  (Torr  cc 
sec  x)  are  compared,  the  latter  unit  being  proportional  to  the  number  of 
molecules  (or  atoms)  passing  a given  cross  section  of  the  tube  per  unit 
interval  of  time. 


39.  Bonczyk,  P.  A. , ’’Determination  of  Einstein  A and  Linewidth  Dependence 
on  Pressure  for  DF  Spectra.  I.  Absorption  Measurements  for  Pure 

DF,"  Phys . Rev. , 11A,  1522,  1975- 

hO.  Bronfin,  3.  R.,  Fifteenth  Symposium  (international)  on  Combustion,  The 
Combustion  Institute,  Pittsburgh,  935-950,  197^+. 

hi.  Oanguli , P.  3.  and  M.  Kaufman,  "The  Rate  of  Homogeneous  Recombination 
of  Fluorine  Atoms,"  Chem.  Phys.  I/?tt.,  2S,  221,  197h. 


The  F concentration  may  he  determined  from  intensities  of  F ESR  spectra 
provided  that  calibration  of  the  ESR  spectrometer  is  done  first  with,  for 
example,  0o.  The  use  of  C>2  for  calibration  is  known  from  earlier  work  by  Ultee 
(Ref.  42)  wherein  a linear  dependence  of  ESR  signal  intensity  on  02  pressure 
is  demonstrated. 


The  relationship  between  species  partial  pressure  and  ESR  signal  inten- 
sity is  contained  in  (Ref.  4-3) 


P 


i 


2(kT)" 

hv@ 


(ID 


In  Eq.  (11),  is  species  partial  pressure,  k is  Eoltzmann's  constant,  T 
is  the  gas  temperature,  h is  Planck's  constant,  v is  the  ESR  transition 
In  Equation  (11),  is  species  partial  pressure,  k is  Boltzmann's  constant, 

T is  the  gas  temperature,  h is  Planck's  constant,  v is  the  ESR  transition 
frequency,  B is  one  Eohr  magneton,  is  a dimensionless  quantity  appropriate 
to  species  i,  f is  a spectrometer  factor  independent  of  i,  is  the 
imaginary  part  of  the  complex  molecular  susceptibility,  and  H is  the  variable 
external  magnetic  field  strength.  The  quantity  [’“y/'dH  is  in  principle 
calculable,  being  independent  of  f.  However,  Aj°  = f^-Z'dH)  is  used  herein, 
which  is  the  area  bounded  by  y."  vs  H,  and  is  referred  to  as  the  ESR  signal 
intensity.  Accordingly,  equation  (11)  may  be  rewritten  in  the  form 
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Appropriate  expressions  in  the  case  of  F and  0.  are 


2(kT)‘ 

hv@ 


QpAp 


(13) 


and 


ssjst*  A 

hvg  \>2  0_ 


d1 


Dividing  equation  (13)  by  equation  (lb)  and  solving  for  Pp,  it  follows 


be  . 


^7 

P-.  = P — — - 

°2  Q02  *02 


(15) 


’Jltee,  C.  J. , "Determination  of  Absolute  Concentrations  by  Gas  Phase 
Meetroa  Spin  Resonance  Spectroscopy/'  J . Appi.  Phys.  37,  1477,  1966. 
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Westenberg,  A. A.,  Progress  in  Reaction  Kinetics  edited  bv  K.  R.  Jennings 
and  j.  Cundall,  Pergammon  Press,  New  York,  Vol.  7/1,  p.  55,  Eq,  20,  lb 


The  quantities  Pq0  ana  A q0  are  known  from  the  spectrometer  calibration.  The 
factors  Qp  and  are  conveniently  tabulated  in  Reference  43.  The  Ap  is  deter- 

mined from  the  measured  F spectra.  Having  Pp  and  the  F volume  flow  rate 
in  cc/sec,  the  torr  cc  sec"-1-  of  F at  the  mixing  region  is  known.  This  is 
then  compared  with  the  torr  cc  sec*-  of  HC1  required  to  extinguish  the  F 
s ignal. 


Apparatus 

A schematic  of  the  apparatus  is  given  in  Figure  149.  The  (Fg  + Ar) 
mixture  is  taken  from  a cylinder  in  which  the  concentration  of  Fp  is  10 
percent.  The  volume  flow  rate  of  this  mixture  in  cc/sec  is  measured  at  1 
atm  pressure  with  a calibrated  rotameter.  The  pressure  at  the  HC1-F  mixing 
region  is  considerably  less  than  1 atm  owing  to  flow  of  the  mixture.  The 
pressure  is  determined  by  interpolating  between  values  read  from  two  Bourdon 
type  gauges  placed  upstream  and  downstream  from  the  mixing  region.  The  volume 
flow  rate  at  the  mixing  region  falls  within  the  range  112  to  334  cc/sec. 

The  flow  tube  has  a cross-sectional  area  of  0.5  'cm^;  hence,  the  molecular 
flow  speed  is  224  to  668  cm/sec. 

The  FP  dissociation  occurs  in  a microwave  discharge  with  the  nominal 
frequency  and  power  being  2450  MHz  and  100  W,  respectively.  The  dissociation 
could  be  made  to  vary  from  15  to  50  percent  while  retaining  an  F signal, 
tut  data  for  less  than  15  percent  dissociation  are  not  obtainable  owing  to  increased 
difficulty  with  the  stability  of  the  discharge. 

The  HC1  apparatus  consists  of  valves  to  control  the  HCl  flow  and  a 
calibrated  mass  flowmeter.  The  Matheson  8116-0250  meter  gives  HCl  flowrate 
directly.  The  calibration  is  checked  before  putting  the  instrument  into  use. 

The  ESR  apparatus  is  a Varian  V-4502  system.  Listing  of  its  components 
is  omitted.  Additions  to  the  system  include  electronics  for  measuring  the 
microwave  frequency  and  power,  and  a NT®  gaussmeter  for  calibration  of  the 
magnetic  field  sweep. 


Results 


(a)  Op  calibration  spectra 


An  Op  line  is  given  inFigurel50.  The  signal  corresponds  to  an  Op  pressure 
of  6.33  torr  as  measured  with  a McCleod  type  manometer.  The  signal  shape  is 
the  derivative  of  a Lorentzian  absorption  signal.  The  derivative  function 
is  due  to  100  KHz  modulation  of  the  variable  magnetic  field  and,  hence, 
the  absorption  line-shape,  along  with  phase-sensitive  detection  of  the  result - 
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FIGURE  149  ESR  F-ATOM  TITRATION  APPARATUS 
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ing  modulated  ESR  signal.  The  calibration  of  the  abscissj.  in  Figure  150  is  done 
with  the  NMR  gaussmeter.  From  the  magnetic  field  value  at  line  center 
*(  = 9U39  G)  and  the  9257  MHz  klystron  frequency,  the  transition  is  identified 
as  K = 3,  J = 4,  Mj  — Mj'  = 3 — 4 in  the  Oo  (3Tg~  ) spectrum.  This  line  is 
identical  with  line  18  in  Table  III  of  Reference  42.  It  occurs  at  a somewhat 
lower  value  of  magnetic  field  since  the  klystron  frequency  in  this  case  is 
lower  than  that  in  Reference  42. 


For  the  observed  spectra,  the  small  modulation  amplitude  used  relative 
to  the  nominal  spectral  linewidth  permits  the  absorption  lineshape  to  be  taken 
as  Lorentzian.  (For  variations  of  this,  see  Ref.  42.)  In  the  Lorentzian 
case,  the  signal  intensity  is  given  by  Reference  42. 
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f - (^) 
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In  Equation  (l6),  W is  signal  width  in  Gauss,  L is  signal  amplitude,  S is  an 
electronic-in-origin  signal  amplification  factor,  and  V is  modulation 
amplitude  in  rms  volts.  Equation  (16)  above  is  equivalent  to  Equation  (18)  of 
Reference  42.  The  W and  L are  identical  with  &H  and  S_p,  respectively  as  given 
in  Figure  1 of  Reference  42 . Since  spectra  are  taken  for  different  values  of  S 
and  V,  we  normalize  all  area  determinations  with  respect  to  S and  V. 


Intensity  vs  02  pressure  is  given  in  Figure  151.  This  is  the  calibration 
curve  for  the  spectrometer.  In  order  for  the  calibration  to  be  applicable 
for  a longer  period  of  use  of  the  spectrometer,  the  Oo  intensities  are 
referred  to  a signal  similar  to  that  in  Figure  150  but. due  instead  to  a ruby 
sample  permanently  fixed  in  the  spectrometer  cavity.  Ey  using  this  signal, 
changes  in  calibration  owing  to  daily  temperature  variations  are  compensated 
for  by  quick  reference  to  the  ruby  signal  rather  than  resorting  to  the  more 
lengthy  procedure  of  O2  recalibration.  Determination  of  a typical  ordinate 
value  of  normalized  intensity  in  Figure  151  which  corresponds  to  a pressure  of 
6.33  torr  is  given  in  detail  in  Example  1 . The  slope  of  the  straight  line 
inFigure  151  is  5.108  x 10“3  torr_'L  . Accordingly,  the  O2  intensity  which 
corresponds  to  1 torr  of  Op  is 


= (5.108  x 10":)A 


ruby 


(17) 


where  the  resulting  A qq  is  used  in  Equation  (15)  in  order  to  determine  ?y.  In 
Equation  (17),  Aruby  is  determined  by  a relationship  among  W,  L,  S,  and  V appro- 
priate to  ruby  and  given  as  in  the  case  of  0?  by  Equation  (16).  In  this  regard, 
see  also  Example  1 . 


Cb)  F-spectra 

The  F-line  used  in  this  work  was  (Mj,  My)  -*  (Mt  , Mr  ) =»  (4,  -)  -*(32,’) 
for  the  E(-F^/2)  state.  This  is  line  "A"  in  Table ' 13,  Reference  13.  .-,s  above, 
this  identification  is  verified  from  the  value  of  magnetic  fieli  at  line 


'-'l 


4b 


EXAMPLE  1 


(Clarification  of  the  Ordinate  Values  of  Figure  15l) 
Values  of  W,  L,  S,  and  V for  a typical  0 line  at  6.33  torr  w^re: 
W = 11.3  (div)  x 1.724  (KHz/div)  x 0.60437  (g/KHz) 

= 11.77  G (div  = 1 unit  of  recorder  paper  length) 

L = 62  div 
S = 262 


V = 38.53  rms  volts 

Accordingly,  from  Equation  (l4)  and  for  this  02  line,  (w2l/S7)  = [(11.77)2  x 62]/ 
r2b2  x 38.53]  = O.S508  G^.  At  6.33  torr,  9 such  .Lines  were  measured  and 
there  resulted  a final,  average  value  of  intensity  given  by,  (T^L)^ 
r average]  = 0-9346  G^.  - 

For  the  ruby  signal,  the  values  of  W,  L,  S,  and  V were:  i 

W = 16.6  (div)  x 4.39  (KHz/div)  x 0.23487  (G/'KHz) 


= 17.12  G 


L = 80 
S = 20.6 

V = 38.38  r.m.s.  volts 


l 


Accordingly,  the  ruby  intensity  is  . = [(17.12)2  x 
29*66  G^.  The  ordinate  of  Figure  151  at  0.^3  torr  is  then 
0.9346/29.66  = 31.5  x 10"3. 


3o]/r20.6  x 38.38]= 
simply  Ao2/Aruby  = 
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center  and  the  measured  klystron  frequency.  A representative  F-line  is  given 
in  Figure  152.  The  area  computation  applied  to  it  is  identical  in  procedure  to 
that  given  for  Op  and  ruby  above.  The  signal  is  appropriate  to  an  F pressure 
of  0.233  torr.  The  details  of  this  computation  are  outlined  inExample2. 

(c)  ( F+HC1) -spectra 

The  effect  of  HC1  titration  on  the  intensity  of  an  F signal  is  shown  in 
Figure  153.  The  HCL  is  introduced  into  the  flow  tube  through  four  0.5  dia  holes 
in  the  wall  of  the  tube  and  at  a point  12.5  cm  upstream  from  the  microwave 
cavity.  The  data  of  Figure  153  correspond  to  an  F pressure  of  0.233  torr  at  a 
flowrate  of  33^  cc/sec,  which  yields  0.233  x 33*+  = 73  [torr  cc  sec*^-]  for  the 
F "flow"  at  the  mixing  region.  From  the  77  [torr  cc  sec'^J  HC1  "flow" 
required  to  extinguish  fully  the  F signal,  it  follows  that  77.-/7S  = 0.99  HCl 
molecules  are  required  to  remove  one  F atom  from  the  flow  stream.  The  linearity 
of  F signal  intensity  with  KC1  flow  and  the  above  ratio  of  nearly  unity 
validate  the  reaction  KC1  + F -•  HF  + Cl  and,  moreover,  the  linearity  excludes 
the  significant  presence  of  interfering  reactions  like  Cl  + F0  - C1F  + F. 

(Ref.  4l)  Table  32  summarizes  titration  results  for  varying  conditions  of 
F and  HCl  flow. 


Conclusion 

Table  32  results  and  the  linearity  of  F signal  intensity  with  HCl  "flow" 
demonstrate  that  the  governing  reaction  is  HCl  + F - HF  + Cl  and,  therefore, 
that  HF  determination  by  absorption  of  light  measurement  gives  the  F 
concentration. 

k 

The  rate  constant  k associated  with  HCl  + F - HF  + Cl  was  determined 
earlier  from  the  decay  of  HF  chemical  laser  signals  and  was  found  to  be  k = 
1.5  x 1013  [cm^  mole--*-  sec--]  (Ref.  44).  The  rate  k'  associated  with  Cl  J- 

Fp  * C1F  + F was  determined  to  be  k'  < 1 x 10 9 [cm3-  mole"^  sec-'1']  (Ref.  4l) 
Hence,  these  latter  results  are  consistent  with  conclusions  reached  herein. 


44.  Kompa,  K.  L.  and  J.  Wanner,  "Study  of  Some  Fluorine  Atom  Reactions  Using 
a Chemical  Laser  Method,"  Chem.  Rhys.  Lett.,  12,  560,  1972. 
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EXAMPLE  2 


(Pp  Determination  for  a Typical  Case) 

For  Figure  152  and  4 other  similar  lines,  it  follows  that  i^J~)  = 
0.8706  G2.  Since  the  ruby  intensity  was  = 29. 08  G2,  it  follows 

that  the  02  intensity  which  corresponds  to  1 torr  of  O2  is 

2 

(^)q2  = (5.108  X 10-)  X 29.06  = 0.1465  G2. 

Next,  it  follows  from  Equation  (13)  that 


P 


F 


2.14  0.8706 

58.6  X 0.1485 


0.214  Torr 


where  values  of  and  0 are  taken  from  Reference  43. 

°2 

The  ?p  above  is  the  pressure  at  the  microwave  cavity,  which  is  12.5  am 
downstream  from  the  region  of  HC1  mixing.  What  is  required  is  pressure  at 
the  mixing  region.  In  order  to  correct  ?_,  to  a point  12.5  cm  upstream  from 
the  cavity,  measurements  of  F intensity  versus  cavity-mixing  separation  are 
done,  which  give  a 9 percent  decrease  in  F intensity  in  going  from  the  point 
of  mixing  to  the  cavity.  Applying  this  correction,  the  F pressure  at  the 
point  of  mixing  is  1.09  x 0.214  = 0.233  torr. 
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FIGURE  152  ATOMIC  FLUORINE  ESR  SIGNAL 
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ATOM  SIGNAL  INTENSITY  (ARBITRARY  UNITS) 


FIGURE  153  HCI  TITRATION  OF  F-ATOM  ESR  SIGNAL 
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APPENDIX  A:  NOMINAL  RATE  PACKAGE 


Chemical  Reactions 

The  possible  reactions  included  in  the  model  and  their  associated  rate 
coefficients  are  listed  in  Table  Al.  For  the  cold  reaction  (l)  the  reaction 
rate  used  is  that  of  Homann  (Ref.  21).  et.  al.  The  reverse  reaction  rate,  and 
that  for  all  reactions  (unless  otnerwise  stated)  is  evaluated  by  using  equili- 
brium constants  derived  from  the  JANAF  Tables  (Ref.  Al).  The  division  of  products 
between  thermal  atoms,  H and  hot  atoms  H*  is  arbitrary.  The  vibrational  pro- 
duct distribution  is  based  upon  the  work  of  Jonathan,  et  al.  (Ref.  A2). 

The  "hot"  reaction  rate  (2)  used  is  from  the  work  of  Albright,  et  al. 

(Ref.  20),  and  the  vibrational  distribution  is  from  Jonathan,  et  al.  (Ref. 

A3).  For  reaction  (3)>  the  "hot"  reaction  with  hot  H atoms,  Albright's 
expression  is  used,  but  with  the  assumption  of  no  activation  energy  for  the 
forward  reaction. 


Coupled  with  the  hot  H-atom  reaction  is  the  deactivation  of  the 
reaction  (4).  Here  the  assumption  is  made  that  the  loss  rate  of  hot 
to  thermal  H atoms  occurs  at  10  percent  of  the  collision  rate  of  the 


hot  atoms , 
H atoms 
atoms  with 


other  gas  molecules . 


The  chain  branching  reactions,  (5)  and  (6),  have  rates  based  on  the  work 
of  Sullivan,  et  al.  (Ref.  11). 


Al. 

Stull,  D.  R. 

2nd  Ed. , U.S. 

A2 . 

Jonathan,  N. , 

Determined  by 

A3. 

J onathan , N . , 

Determined  by 

and  Prophet,  H. , Project  Directors  JANAF  Thermochemical  Tables, 
Dept,  of  Commerce,  National  Bureau  of  Standards,  June  1971. 

et  al.,  "Initial  Vibrational  Energy  Level  Distributions 
Infrared  Chemiluminescence,"  Mol.  Fhys.,  22,  561,  1971. 

et  al.,  "Initial  Vibrational  Energy  Level  Distributions 
Infrared  Chemiluminescence,"  Mol.  Phys.,  24,  1143,  1972. 
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TA3LE  Al.  CHEMICAL  REACTIONS  CONSIDERED  AND  NOMINAL  PATES 


f/f 

1.  F + H2  «•  HF(v)  + H,  H*  AH  = -31.6  kcal/mole 
kf  = 2.63  x 10-^  exp  (-I600/RT)  cc/tnolecule  sec 
kb  = 2.95  x 10"10  exp  (-33200/RT)  cc/molecule  sec 

(fv,  v = 0,11)  = 0,  .1639,  .5574,  .2787,  0,  0,  0,  0,  0,  0,  0,  0 

[Note  - fractional  division  between  H and  H*  (hot  atoms)  is  arbitrary] 
fvkf 

2.  H + F2  - HF(v)  + F AH  = -98  kcal/mole 

kf  = 2 x IQ'10  exp  (-2400/RT)  cc/moleouie  sec 

^b  = 3.29  x 10  exp  (-100h00/RT)  cc/molecule  sec 

(fv,  v = 0,11)  = .012,  .027,  .033,  .039,  .13^,.  .265,  .298,  .134,  .058,  0,  0,  0 
* kffv 

3.  H + F2  - HF(v)  + F AH  + -129.6  kcal/mole 

kf  = 2 x 10-^0  cc/molecule  sec 

kb  = 3.29  x 10"^-  exp  (-I296OO/RT)  cc/molecule  sec 
fy  arbitrary,  usually  taken  as  for  rctn  1. 

4.  H*  + M - H + M AH  = -31.6  kcal/mole 

kf  = 2.44  x 10"~2T^/2  cc/molecule  sec 

kfa  » 0 
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TABLE  Al.  CHEMICAL  REACTIONS  CONSIDERED  AND  NOMINAL  RATES  (CONT'D) 


5.  HF (v>3 ) + F2  - HF  (0)  + 2F  AH  = 37.7  kcal/mole 

kf  = 3.5  x lO"1?  cc/molecule  sec 

kb  - 0 

6.  H2 (vy.)  = F2  - HF(O)  + H + F AH  = 5.82  kcal/mole 
kf  = 8.0  x 10“^-9  cc/molecule  sec 

kb  = 0 

7.  F + F + M-F2+M  AH  = -37.7  kcal/mole 

kf  = 4.99  x io-33t-°-5  (cc/molecule  sec)2 

kb  = 8.32  x 10"11  exp  (-35100/RT)  cc/molecule  sec 

8.  H+H+M-Hg+M  AH  = -104.2  kcal/mole 

kf  = 2.75  x 10"3®T~^  (cc/molecule  sec)^ 

kjj  = 6.96  x 10‘6T-1-025  eXp  (-104200/RT)  cc/molecule  sec 

9.  H + F + M - HF(0)  + M AH  = -129.6  kcal/mole 

kf  = 6.6  x 10“3°t"'3“0  (cc/molecule  sec)2 

kb  = 7.78  x 10“^T~^  exp  (-129600/RT)  cc/molecule  sec 


26l 


The  F2  dissociation  reaction  rate  (7)  is  based  on  the  measurements  of 
Shui  (Ref.  A4),  et  al.  That  for  Hg  (8)  is  based  on  the  work  of  Jacobs,  et  al. 
(Ref.  A5).  Finally,  the  HF  dissociation  (9)  rate  coefficient  is  based  on  the 
work  of  Blaur,  et  al.  (Ref.  A6)  with  a revised  activation  energy. 


V-T  Relaxation 

The  relaxation  times  for  the  1-0  V-T  transition  are  listed  in  Table  A2. 
In  all  cases  a Landau-Teller  fit  has  been  used.  For  HF-HF  deactivation  it  has 
been  necessary  to  fit  the  data  in  four  temperature  ranges  due  to  the  anomalous 
temperature  dependence  exhibited.  The  fit  is  based  upon  the  composite  data 
from  many  sources  as  summarized  in  Cohen  (Refs.  28  and  29). 

The  Hg  V-T  rate  is  based  on  the  upper  limit  deduced  by  Bctt  (Ref.  A7). 
Finally,  the  Fg,  Ar,  and  He  relaxation  times  are  conservative  fits  to  the 
composite  data  summarized  in  Cohen  (Refs.  28  and  29). 


The  special  cases  of  relaxation  by  H and  F atoms  have  been  treated  as 
multiple-quanta  transition  processes.  The  relaxation  times,  based  on 
Wilkins  (Ref.  27)  calculations  are  summarized  in  Table  A3. 


V-V  Relaxation 

The  V-V  formulation  for  HF-HF  self  relaxation  is  based  primarily  on  the 
long-range  dipole-dipole  interaction  model  of  Sharma  (Ref.  A8)  as  adaoted  by 


A4.  Shu,  V.  H. , et  al.,  Thirteenth  Symposium  (International)  on  Combustion, 
Combustion  Institute,  Pittsburgh,  1971. 

A5.  Jacobs,  T.  A.,  et  al.,  "Kinetics  of  Hydrogen  Halides  in  Shock  Waves, 

II.  A New  Measurement  of  the  Hydrogen  Dissociation  Rate,"  J.  Chem.  Phys., 
47,  54,  1967. 

A6.  Blauer,  J. , "The  Kinetics  of  Dissociation  of  Hydrogen  Fluoride  behind 
Incident  Shock  Waves,"  J.  Phys,  Chem.,  72,  79,  1968. 

A7.  Bott,  J.  F. , "Vibrational  Relaxation  of  HF  (V  = 1)  and  DF  (V  = 1)  by  Hg 
and  Dg,"  J.  Chem.  Phys.,  6l,  2530,  1974. 

AO.  Sharma,  R.  D.  and  H.  Schlossberg,  "On  the  Assumptions  Underlying  the 

Hear  Resonant  Theory  of  Inergy  Transfer,”  Chem.  ih~,s.  Lett.,  CC,  5,  1973- 
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TABLE  A2.  V-T  RELAXATION  TIMES 


1.  HF-HF  - k , = v k 

v,v-l  10 

Anomalous  temp  dependence  treated  by  breaking 
up  data  into  1+  separate  regions  and  fitting  each 
to  a Landau- Teller  form: 

-l/R 

PT  = a exp  (bT  ) =■  l/k 


Temp  Range 

a (sec-atm) 

b 

T a 1371  °K 

.876  x 10“9 

65.6 

1351  2 T 2 751 

.223  x 10'7 

29.5 

751  2 T 2 578 

.87  x lO-3 

-66 

578  2 T 

.0256 

-9^.2 

HF-M  - k = vk 

v,v-l 

10 

PT  = a exp 

(bT-1^3)  = l/k 

M 

a(sec-atm) 

b 

H2 

3.5  x 10"° 

a 

0 

* O 

5.929  X 10“ 

79.*+ 

At 

-8 

5.929  x 10 

79. ^ 

-8 

tie 

5.929  x 10 

79.  ^ 
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TAB  Li  A3.  V-T  MULTIPLE-QUANTA  TRANSITIONS:  v v-n,  n ARBITRARY 


For  H and  F Atoms 

_ v«  = aT  b exp(-c/RT)  sec  ^ atm  ^ 

(Note  values  for  a in  following  table  are  multiplied  by  10“12) 
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Rockwood  (Ref.  A9).  Both  short  and  long  range  terms  are  included  such  that 
the  net  rate  coefficient  k£;|It  = L§;§tl  + S^’^.  The  short  range  term  is 
given  by 


where 


~s,s+l 
°r ,r-l 


(7^-)  ( 

1-er 


l-t.Il).>  SS  8xp  <15rs/21c)  F««"W 


(Al) 


e 


CUT 


z 


z 

d 


u 

k 

M 


Z 


Molecular  Anhsirmonicity 


1/2  (3  - exp (-2/3  [cut])) 
2TT2  . .1/2 


exp 


A- 


■rs 


(— Z—Y 

v2rrkT ' 


(-2/3  [our]) 


u,kT 

&nZV?(v10)Z&2 

d2(8rrkT  )i/2 
y» 

gas-kinetic  collision  diameter  = 2.7  x 10“'-'’  cm 

reduced  mass  of  collision  pair 

Boltzmann's  constant 

reduced  mass  of  KF  oscillator 

frequency  of  fundamental  transition 

range  parameter  for  repulsive  interaction. 


It  turns  out  that  this  term  makes  only  a small  contribution. 


The  long  range  term  is  given  by 


. s ,s+l 
r,r-l 


s+1 


1-er'  l-(s+l)e 


<rb><i 


lw  ^ <*! 


*rs 


;2kT)  exp  (-(^S 


rs- 


A)2) 


4.A2 ) 


A9. 


Rockwood,  S.  D.,  et  al.,  "Time -Dependent  Calculations  of  Carbon  Monoxide 
Laser  Kinetics,"  IE55  J,  Quant.  Sleet.,  £>  120,  1973. 
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.'here 


and 


ii.i 

z T 

2 1 /2 
* = BTX' 


The  parameters  that  can  be  adjusted  in  the  model  are  the  range  parameter 
l used  in  the  short  range  term,  A,  and  a2.  The  values  we  have  used  are 

l = 2 x 10"9  cm  (A3) 

A = 27°K  (AU) 

62  = 10-25  (A5) 

These  parameters  were  selected  such  that  the  rate  of  the  transition 

HF(0)  + HF(2)  - 2 HF(1)  (A6) 

is  pegged  to  the  value  reported  by  Airy  and  Fried  (Ref.  33).  Thus 


at  300  n. 


= K = 1.82  x 10  ^ cc/molecule-sec 


(A7) 


A similar  treatment  is  used  for  HF-H2  V-V  interactions.  Here  the  short 
range  term  is  taken  as  zero, 

L10  ■ 5 x 10"13  (A6> 


6"  = 1.6  x 10"-°, 


w 


in  which  the  direction  of  the  transition  is  endothermic.  The  exothermic  rates 
are  found  from  detail  balance.  The  values  used  match  the  rate  for  the  exo- 
thermic reaction 


H2(l)  + HF(0)  - Hg(0)  + HF( l) 


{ A10 ) 


?67 


reported  by  Bott  and  Cohen  (Ref.  A10)  (i.e.,  1.33  x 10 


cc/molecule-sec ) . 


-12 


V-T  Relaxation 

The  rate  coefficients  employed  are  derived  from  the  recommendations  cf 
Cohen  (Refs.  28  and  29).  They  are  summarized  in  Table  A4. 


A10.  Bott,  J.  ?.  and  N.  Cohen,  "Temperature  Dependence  of  V-V  and  V-R,  T 

Energy  Transfer  Measurements  in  Mixtures  Containing  HE,"  J.  Chem.  Phys., 
^539,  1973- 
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.APPENDIX  B 


STEADY- STATE  SOLUTION 


In  complex  systems  such  as  HF  chain  reaction  lasers  a potentially  useful 
but  often  neglected  analytic  tool  is  the  development  of  approximate  closed 
form  solutions.  For  example,  in  typical  chain  reaction  conditions  computer 
calculations  indicate,  in  the  absence  of  significant  atom  loss  processes,  a 
quasi-steady  state  develops.  In  this  steady-state  condition  the  number  of 
H and  F atoms  remain  essentially  constant.  It  happens  that  this  situation 
is  amenable  to  a closed  form  solution. 

Consider  a situation  in  which  the  relative  rate  of  change  of  H2  and  F? 
is  slow  compared  to  the  time  scale  over  which  steady  state  is  achieved. 

Assume  also  that  the  temperature  varies  slowly  and  that  only  the  forward  hot 

kin  y 

and  cold  reactions  are  important  (i.e.,  H + ?0  HF  + F and  F + H2  "S'HF+H). 
Equations  for  ^he  rate  of  change  of  H and  F may  be  written: 


[H]  = kc[H2MF]  - (31) 

[F]  = -[H]  (32) 


Now  k,,[H2]  and  kv^Fo]  may  be  taken  as  constants,  thus  allowing  the  simple 
solution  of  '■wo  linear  differential  equations  with  constant  coefficients. 

The  initial  conditions  are  [F]  = [F]q,  [H]  = 0,  [F]  = k [H2][F]  , [H]  = -[F] . 
The  solution  becomes 


[F]  _ Kh  * KC  exo  (-[Kc  + Kh]t) 

[Flo  = Kc  + ^ 


fH~l  _ Kc  - Kh  exp  (-The  + Kh]t)  (~4) 

[F]  = K + K. 

L Jo  c n 


where  K q = kc[H2]  .and  kh=[F2]. 
s‘a‘e  solu'- ions^are  achieved: 


As  (K  + K )t  becomes  larger, 
- h 


he 


■ eady 


TFl  _ khrF?  i 

[F]q  ’Kh[*2l  + kcLH2l 


fHl  _ kofij^] 

[F]0  = kh[F2]  + keCHg] 


(55^ 
(B  6) 
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The  rate  of  increase  of  HF  is  simply 


[HF]  = 2 


kcfH2lkh[F2l 

W - W 


(B7) 


The  above  solution  could  have  equally  well  been  achieved  by  equating  31 
to  zero  wi*-h  the  constraint  [F]  = [H]  + [F]. 


Comparison  of  this  solution  with  computer  calculations  demonstrates 
excellent  agreement  (less  than  1 percent  error)  when  all  back  reactions 
are  set  equal  to  zero,  and  local  values  are  used  for  k^,  k^,  [H^l  and  [Fp]. 
Inclusion  of  the  back  cold  reaction  in  computer  calculations  can  lead  to 
significant  error,  depending  upon  the  rates  used  for  the  upper  levels  of  HF. 
Consequently  it  has  been  necessary  to  include  the  back  cold  reaction  in  the 
formulation.  This  is  easily  done.  The  rate  equations  become 

[H]  = kc[H2]*[F]  - (^[Fg]  + Z [HF.]k^)[H]  (b8) 

[F]  = '-[H]  (B9) 

Equating (B8)  t0  zero  with  [H]  + [F]  = [F]  leads  to 


[F]  _V_ 

[F]0  “ D 

Wo‘  D 

[HF]  = 2 rp-j 


(310) 


(311) 

( 312 ) 


where 


D = Kc  + Kh  + ZCHF^k^ 

(313) 

and 

K = K + r[KF. ]k 
h h 1 ci 

(31^) 

K and  K are  as  defined  »bove. 
n ^ 
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One  further  refinement  is  possible.  A simplified  diffusive  wall  loss 
model  may  be  treated.  Suppose  H atoms  are  lost  to  the  walls  at  fhe  ra*-e 


[SJw  = - *D  LHJ  (315) 

If  the  rate  of  change  of  F^  and  is  slow  we  again  may  treat  the  equations 
as  linear.  If  we  further  assume  that  the  loss  of  atoms  is  slow  relative  to 
the  other  reactions  (i.e.,  k^/D  = 6 «l)  the  equations  may  be  solved  simply 
*:o  give  the  quasi-steady  solution 


22  - a 


-K  61 
e c 


(316) 


(317) 


[HF] 


[F]o  :<hKc 


D 


2 


(B13) 


Comparison  of  *his  model  with  computer  calculations  is  excellent. 

Tables  31  and  B2  show  comparisons  of  the  model  with  a corresponding  com- 
puter calculation  for  cases  with  and  without  the  loss  of  H atoms  to  t-he  walls 
respectively.  Local  values  are  used  for  [F?]  and  [HQ]. 

Thus  we  see  tha^  this  simple  model  can  serve  well  as  a means  of  gaining 
quick  InsighT  into  experimen*al  results  and  provide  guidance  for  efficient 
use  of  the  more  sophisticated  computer  models. 


TA3LE  Bl.  STEADY  STATE  SOIXJTION  - COMPARISON 
WITH  COMPUTER  CALCULATION 


Case:  Back  react  ion  included 


x = 12.5  cm 


T = 312 


[F]0  = 1.53  x 10 

[H2]  = 1.1*0  x 1015 

[F2]  = 5.39  x 1015 


Calculation 

— 

Approximation 

Error 

(%) 

[K] 

6.23  x 1012 

6.l6 

X 

1012 

1.1 

[F] 

9.09  x 1012 

9.15 

X 

1012 

0.6 

[HF] 

17 

2.77  x 10  1 

2.Jh 

X 

101T 

1.1 

07U 


TABLE  B2.  STEADY  STATE  SOLUTION  - COMPARISON 
WITH  COMPUTER  CALCULATION 


Case:  Diffusion  included 


™Diff  ' -15°  M 


x = 25  cm 
T = 334°K 


[F]  = 1.434  x 1013 
[H2]  = 1.178  x 1015 

[F2]  = 4.911  x 1015 


Computer Approximation  Error 


[H] 

4.655  x 1012 

4.483  x 1012 

3.8$ 

[F] 

8.023  x 1012 

8.007  x 1012 

0.2 % 

[H+F] 

1.268  x lO1^ 

1.249  x 1013 

i.H 

[HF] 

2.415  x 101'  sec"1 

2.330  x lO1" 

iM 
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APPENDIX  C 


ACCOUNTING  FOR  TWO-DIMENSIONAL  VELOCITY  VARIATIONS  IN  PIPE  FLOWS 


The  use  of  a quasi-one-dimens ional  model  to  analyze  reacting  flows  within 
a circular  pipe  can  lead  to  errors  in  interpretation  of  the  data.  First,  the 
presence  of  boundary  layers  can  strongly  influence  the  pressure  and  mainstream 
velocity.  A simple  minded  attempt  to  employ  the  true  pressure  distribution 
and  mainstream  velocity  distribution  in  a one-dimensional  calculator,  however, 
will  still  lead  to  erroneous  conclusions  about  the  rates  one  is  attempting  to 
evaluate.  The  flow  residing  in  the  slower  moving  boundary  layer  tends  to 
increase  the  effective  reaction  time  for  the  gas  found  at  a particular  axial 
location.  Thus  a simple  .analysis  based  on  the  mainstream  velocity  would  tend 
to  overestimate  the  rate  coefficient  for  the  reaction  of  interest. 

It  is  possible,  however,  to  account  for  the  effective  r te  of  reaction  in 
a two-dimensional  flow  with  a one  dimensional  model.  Under  certain  assumptions, 
a multiplying  factor  for  the  rate  coefficients  may  be  defined,  based  upon  the 
2-D  velocity  profile,  that  is  incorporated  into  the  rate  equations. 

Consider  the  following  situation:  (l)  The  rate  equation  can  be  expressed 

as  linear  function  of  the  species  of  interest.  (2)  The  velocity  profile  can  be 
expressed  as  u = u^  ( 1 - (r/R)n),  where  is  the  center  line  velocity  and 

R the  radius  of  the  tube.  (3)  Diffusion  is  treated  as  for  a binary  system. 

The  first  assumption  is  applicable  to  most  relaxation  rate  equations  and 
is  also  approximately  applicable  to  the  production  of  HF  in  the  absence  of  back 
reactions . Then  we  have 


HF  = [Fj[Hoj  kc  + [H][F2j  kjj 


(Cl) 


where  kn  and  kc  are  the  hot  and  cold  reaction  rate  coefficients.  We  can  re-express 
^H2j  and  as(i.H2]0  - [HFj/2)  and  ([F2j0  - [HFj/2).  This  yields  the  rate 

equation 


1HF] 


= (Te  TE  1 If  TrIFE  1 If.  1 (-•-  ~ An-i  rRVl 

<■  -*•  J L «*0  ^2  * L^j  L^2  Jq  j “ p L^-  J 


(C2) 


For  quasi-steady  state  situations  the  terms  in  braces  will  be  slowly  varying 
and  may  be  treated  as  constants  so  that 


[HFj  = Kx  - Kg  [HFj 


(C3) 


The  second  assumption  reasonably  represents  the  profiles  of  a pipe  flow 
in  various  stages  of  development  from  an  initially  uniform  velocity.  Arbitrary 
profiles  can  also  be  treated  as  will  be  shown  later. 

Since  our  flow  is  primarily  diluent,  the  third  assumption  is  certainly 
valid. 


We  may  now  write  the  rate  equation  for  an  arbitrary  species  as 

/,  3s  D 19  /»  3s,  t -r  \ 

( 1 - r ) tt  = — ? - 7 (r  7 + (K-s) 

9x  k?.^  r 9r  dr 


whe  re 


r = r/R 
x = kx/u£ 
u = u£  (1  - ?n) 


'rate  proc 


= k(  K-s ) 


We  impose  the  initial  condition 


and  the  boundary  condition 


s = sQ  at  x = 0 


, D 9s  „ 

k s = - - — : at  r = 1 
w r 9r 


(cU) 

(C5) 

(c6) 

(C7) 

(08) 

(C9) 

(CIO) 


Let  us  now  consider  two  limiting  cases:  (l)  infinitely  fast  diffusion 

(2)  infinitely  slow  diffusion 

Infinitely  Fast  Diffusion 

For  this  case  we  define 


and 


kR2 


e = 


0 


Then  we  assume  the  rate  equation  has  a solution  of  the  form 

s = cpQ  + ecp1  + 0(  e2 ) 

Substitution  of  this  into  Equation  (cU)  yields  the  solution 

3 = *o(*)  + + i (If1  - K + " 777 


(n+2)^  9x 


?n+2  + o(e2){ 


(Cll) 


(C12) 


(C12) 


278 


Application  of  the  boundary  condition  with  = 0 gives 

- n+2  ~ 

c?q(x)  = K + C e n X 

We  may  now  apply  the  initial  condition  to  find  that  as  e -»  0 

. . . n+2  k i 

s = K + (sn-K)  exp  [-  " xj 


(ClU) 


(C15) 


This  solution  is  of  the  same  form  as  the  solution  for  a one -dimensional  flow 
at  velocity  u^  . The  only  difference  is  that  the  exponential  argument  contains 
the  additional  factor  (n+2)/2.  Thus  we  consider  the  solution  of  our  problem 
to  be  the  solution  of  the  one -dimensional  problem  referenced  to  the  centerline 
velocity,  but  with  an  effective  rate  coefficient  defined  by 


n+2  „ , . 

kgff  = k = ft»(n)k 

Cil  n 


(Cl6) 


where  3®(n)  is  defined  as  the  effective  rate  multiplier  for  infinite  diffusion. 
In  the  limit  of  fully  developed  flow  (n  = 2)  this  reduces  to 


keff  = 2k 

a solution  previously  obtained  numerically  (Ref.  Cl). 


(C17) 


Strictly  speaking  this  solution  applies  to  a case  in  which  the  velocity 
profile  is  invariant  with  axial  distance.  However,  we  may  apply  the  solution 
to  a developing  flow  in  which  the  fit  parameter  n varies  by  using  the  value 
for  Soo(n)  appropriate  to  the  local  velocity  profile.  In  the  instantaneous 
diffusion  limit,  this  approach  is  rigorously  valid  since  the  initial  condition 
for  every  differential  interval  is  a uniform  concentration  profile. 

Infinitely  Slow  Diffusion 

We  now  turn  our  attention  to  the  case  of  no  diffusion  (i.e.,  D -+o).  In 
this  limit  we  must  solve  the  simplified  equation 


( 1 - r ) —w  = K-s 
Ox 


(C18) 


Since  there  is  no  communication  among  streamlines  we  find  an  independent  solution 
along  each: 


s(?)  = K + (sQ-K)  e 


(C19) 


Cl.  Poirier,  R.  V.  and  R.  W.  Carr:  "The  Use  of  Tubular  Flow  Reactors  for 

Kinetic  Studies  over  Extended  Pressure  Ranges,”  J . Phys . Chem . 75 , 1593, 

1971. 
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An  effective  rate  coefficient  may  be  defined  by  integrating  s(r)  and  equating 
it  to  the  one -dimensional  form  of  the  solution.  Thus  we  have 

keff  = - ^ (j  2r  exp  (— ~)  dr]  (C20) 

0 

and 

ec(n,x)  = keff/k  (C21) 

A plot  of  the  variation  of  0O  with  n for  various  values  of  x is  shown  in 
Figure  Cl.  We  see  that  as  the  extent  of  reaction  increases  the  rate  multiplier 
for  any  given  velocity  profile  decreases  towards  unity.  This  is  a consequence 
of  the  lack  of  communication  between  streamlines.  As  we  move  further  and  further 
downstream  the  reaction  on  the  outer  streamlines  approaches  completion  and  the 
effective  rate  of  reaction  becomes  more  a reflection  of  the  chemistry  near  the 
center  of  the  flow. 

As  for  the  infinite  diffusion  case,  the  solution  above  applies  for  a 
particular  velocity  profile,  invariant  with  axial  distance.  It  is  possible  to 
apply  it  to  developing  flows  by  using  locally  appropriate  values  for  30(n,x). 

This  approach  will  not  be  rigorously  correct  in  this  case  since  the  initial 
condition  for  any  interval  will  not  be  a uniform  concentration  profile.  Since  our 
interest,  however,  is  in  bracketing  the  true  solution,  such  an  approach  will 
suffice.  The  tendency  will  be  to  underestimate  30  and  thus  overestimate  the 
true  rate  coefficient  k in  an  analysis  of  experimental  data.  The  instantaneous 
diffusion  solution,  on  the  other  hand,  would  underestimate  k.  Hence  the  two 
solutions  should  effectively  bracket  the  results. 

Test  of  Model 

To  test  this  approach  to  simulating  reacting  pipe  flows  with  a one- 
dimensional model,  a comparison  was  made  with  a two-dimensional  computer  code, 
MARBL,  for  a relatively  simple  case  involving  the  chain  reaction  without  any 
vibrational  level  detail.  Two  cases  were  considered  (Table  Cl).  Each  had 
the  same  concentration  of  reactants,  but  in  one  case  had  an  additional  15  torr 
of  argon  to  bring  the  initial  pressure  up  to  20  torr.  The  calculated  velocity 
profiles  from  the  2-D  MARBL  program  were  used  to  develop  3 vs  x curves  for  use 
in  the  1-D  CHEMIZ  code. 

The  results  of  the  calculations  are  shown  in  Figures  C2  through  CL. 

The  first  two  show  the  HF  variation  for  the  5 and  20  torr  cases.  The  area 
between  the  infinite  and  zero  diffusion  solutions  is  shaded.  Also  shown  are 
the  solutions  for  two  more  simple  minded  one-dimensional  approaches -a  constant 
area  solution  and  specification  of  the  centerline  velocity  distribution.  We 
see  that  the  effective  rate  coefficient  approach  provides  the  best  agreement 
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TABLE  cl.  INPUT  CONDITIONS  - EFFECTIVE  RATE  TEST  CASES 


T = 300°K 
= 0 

u = lOOOO  cm/sec 


Case 

*f2 

P=5 

0.0025 

0.009 

p=20 

0.000625 

0.00225 

XJ 

0.000025 


*Ar 

0.98825 


0.00000625 


0.997063 


o m o 

<n  , - 

aandinnw  aiva 
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5 TORR  CASE 
HF  PRODUCTION 


FIGURE  C2  ONE  DIMENSIONAL  SIMULATION  OF  REACTING  PIPE  FLOW 


203 


20  TORR  CASE 
HF  PRODUCTION 


FIGURE  C3  ONE  DIMENSIONAL  SIMULATION  OF  REACTING  PIPE  FLOW 


2& 


in  both  cases.  This  is  also  true  for  the  predicted  temperature  as  is  shown  in 
the  third  figure . 

Application  to  Experiment 

In  the  actual  experimental  situation  the  velocity  profiles  are  not  the 
smooth  shapes  predicted  by  the  2-D  MARBL  program  (Figs.  6 - " ).  These 

profiles  cannot  be  easily  approximated  by  a profile  of  the  form  u = Uc  (l  - rn). 
A series  summation  of  the  form 


U ■ u£  <1  - f”) 


(C22) 


can  reasonably  fit  all  the  measured  profiles,  however.  The  solution  developed 
can  be  easily  extended  to  such  a case,  giving  the  result 


from  which  an  effective 


1 

soo = 

n+2 

fit  parameter  n can  be  defined  as 

61  j. 


2 

neff  = 3»-l 


(C23) 


(C2U) 


This  value  may  then  be  used  to  determine  the  appropriate  value  for  30« 


Based  upon  the  experimentally  measured  velocity  profiles,  the  effective 
rate  multipliers  have  been  calculated.  The  variation  of  3*  with  x is  shown 
in  Figure  C5.  The  values  for  80  were  not  significantly  lower  and  it  was 
found  that  little  predicted  difference  was  found  between  the  two  possible 
solutions.  Thus,  for  the  majority  of  calculations  the  infinite  diffusion 
solution  was  used  with  the  understanding  that  derived  rate  coefficients  might 
be  up  to  5/J  low. 


References 


> 


Cl. 


t 


Poirier,  R.  V.  and  R.  W.  Carr:  "Hie  Use  of  Tubular  Flow  Reactors  for 

Kinetic  Studies  over  Extended  Pressure  Ranges,"  J.  Fhys.  Chem.  75,  1593, 

1971- 


288 


r-xr-M 


APPENDIX  D 


ESTIMATES  OF  WALL  ADSORPTION  EFFECTS 


The  experimental  data  available  from  the  flow-tube  and  pulsed-fluorescence 
experiments  has  indicated  that  wall  interactions  may  not  be  negligible.  In 
particular  the  potential  loss  of  atoms  (H  or  F)  or  HF  on  the  walls  can  influenc 
the  interpretation  of  the  data.  It  is  thus  instructive  to  try  to  estimate  the 
maximum  possible  loss  by  such  processes. 


Suppose  that  any  molecule  of  the  species  of  interest  that  strikes  the 
walls  of  the  test  cell  is  adsorbed  (or  in  effect  destroyed).  Suppose  also  that 
our  cylindrically  shaped  test  cells  are  very  long  relative  to  their  diameter, 
and  thus  can  be  idealized  as  an  infinite  cylinder.  Then  the  rate  of  change  of 
the  subject  species  will  be  described  by  the  equation 


5n( r,t) 


D (1  ln(rr-) 

' r dr 


with  the  boundary  condition 


(Dl) 


where 


n(R,t)  =0 

n = the  number  density  of  the  subject  species 


(D2) 


n = the 

chem 


rate  of  production  by  chemical  reactions 


^ ~ the  dit  fusion  coetficient  o*  the  subject  species  relative 
to  the  entire  gas 

Solution  of  this  equation  will,  in  principle,  give  the  maximum  adsorption 
that  could  occur  for  the  conditions  specified.  Let  us  consider  two  highly 
simplified  cases:  (l)  the  species  concentration  is  initially  constant,  with 

a zero  chemical  production  rate;  (l)  the  species  concentration  is  initially 
zero  and  the  chemical  production  rate  is  a constant. 


The  first  case  approximates  the  disappearance  of  the  H or  F atoms  which 
would  otherwise  maintain  essentially  constant  concentrations  over  a large 
fraction  of  the  reaction  time.  The  second  case,  if  we  restrict  attention  to 
the  characteristic  reaction  time  can  approximate  the  loss  of  the  HF  being 
produced  to  the  walls. 


Potential  Loss  of  Atoms  to  the  Walls 


We  now  obtain  the  solution  for  the  first  case.  All  that  is  required  is 
the  solution  of  the  homogeneous  part  of  equation  (Dl) . This  standard  solution 
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may  be  found  in  any  textbook  on  differential  equations  (Ref.  Dl)  and  is 
given  by 

,R  - X 2©t 

n (r,t)  = -2  1 Jn  No(r)  Jp(-<nr/R)  dr  j (xn£)  g (D3 

n=  1 [Ji(xn)]2  0 R 

where  NQ(r)  is  the  initial  distribution  of  the  species  considered,  JQ  is  the 
zeroeth  order  Bessel  function,  is  the  first  order  Bessel  function,  and  xr 
is  the  n-th  2ero  of  JQ. 

If  we  let  N0(r)  = rig  = const,  we  may  integrate  equation  (D3)  to  obtain 


n(r,t)  = 2n0  I 


Jo(xnr/R) 


Xj^Dt 


1 xnJl(*nr/R) 

This  equation  can  be  nondimens ionali zed  by  letting 

r = r/R 

« .i&t 


Then  we  have 


whe  re 


n = n/n0. 

£=21  CnJ0(xnr)  e X1 


n xnJl^xrJ 

With  this  nondime ns ional  form  of  the  equation  it  is  possible  to  develop 
•universal  curves  of  the  species  distribution  as  a function  of  time.  Then,  with 
knowledge  of  the  cylinder  radius  and  gas  diffusivity  D,  the  time  history  of  a 
specific  case  can  be  evaluated.  The  concentration  profiles  may  also  be  averaged 
over  the  cylinder  for  various  times  to  provide  a plot  of  the  average  concen- 
tration in  the  cell  as  a function  of  time.  This  has  indeed  been  done;  the 
results  are  shown  in  Figure  Dl.  Given  are  the  centerline  variation,  the 
variation  of  the  average  along  a diameter  (called  the  linear  average)  and 
the  average  over  a circular  cross  section  (the  cylindrical  average).  The 
linear  average  is  essentially  what  one  might  observe  using  some  radiative  probe 
technique,  whereas  the  cylindrical  average  represents  the  true  average  concen- 
tration for  this  idealized  case. 


Dl.  Martin,  W.  T.  and  E.  Reissner,  Elementary  Differential  Equation,  Addisor.- 
Wesley,  Reading,  MA,  1961. 
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Consider  now  the  cases  of  interest  to  us.  For  the  pulsed- fluorescence 
experiment,  the  H atoms  are  immersed^in  a bath  of  Ar  at  100  torr.  For  this 
the  diffussivity  is  roughly  D =“  8 em~/sec.  The  radius  is  0.6  cm;  hence,  we 
have 

£ = 122  t (d6) 

The  time  scale  of  the  experiments  is  from  10  to  40  milliseconds.  Thus  the 
normalized  times  range  from  1.2  to  5*0.  We  see  that  over  such  a time  scale  a 
considerable  fraction  of  the  H atoms  can  be  lost. 

For  the  flow  tube  experiment  the  argon  pressure  is  only  20  torr,  yielding 
a diffussivity  D =»  40  cm^/sec.  Since  the  test  cell  radius  is  2.5  cm  we  have 

£ = 37  t (D7) 

Now  the  time  scale  of  a fluid  element  in  the  flow  tube  experiments  over  the 
range  of  measurements  is  4 to  6 milliseconds , giving  normalized  times  from 
0.15  to  0.22.  Even  for  these  relatively  short  times  30$  or  more  of  the  atoms 
can  potentially  be  lest.  Hence,  atomic  loss  to  the  walls  may  be  important 
for  both  experiments . 

Potential  Loss  of  HF  to  the  Walls 


We  now  consider  the  case  in  which  there  is  a constant  production  rate  K 
over  some  time  interval  t.  The  equation  to  solve  is  then 

/ 1 jn(r,~)  5^n(r,  th  (D8) 

at  V r dr  dr2  r K 

The  solution  to  the  homogeneous  part  of  this  is  given  by  equation  (D3).  In 
addition  we  obtain  the  particular  solution  for  a steady  state  condition  (dn/dt  = 
0): 

np(  r)  = C0  + (?!  Xnr  + Cgr2  (D9) 

Since  n is  finite  at  r = 0,  we  set  C-,  =0.  The  boundary  condition  (n  = 0 on 
r = R)  and  the  differential  equation  yield  the  complete  solution 

■ ^ (ff'-r2)  (DIO) 


?92 


■ 1 T 


«■'  • 


We  may  superpose  the  two  solutions  and  satisfy  the  initial  condition 


oy  letting 


n(r,o)  = nQ(r) 


"o  = - IS 


(Dll) 

(D12) 


in  the  homogeneous  part.  Thus  we  have 

[4Dr 


nQ(r)  + r(r2-R2)  J0(i^L)dr 


'Xnr, 


CJl(*n» 


- xn" Dt 


^0(¥)e  r2 


(D13) 


In  prineiple'we  could  apply,  this  in  successive  segments,  each  having  a 
different  average  chemical  production  rate  K,  with  the  final  profile  of  one 
segment  used  as  the  initial  condition  for  the  next.  For  our  purposes,  however, 
it  will  only  be  necessary  to  consider  the  simple  case  for  which  n0  = 0,  using 
an  average  chemical  rate  over  the  duration  of  the  reaction. 


For  the  case  N0  = 0 equation  (D13 ) car.  be  easily  integrated  to  give 

, X L2  2,  0n2  “ J0(xnr/R)  ?y  ! 

n(r,t)  = H7;  ‘k(R  "r  ) - 8R  1 -^rrr—r  e ( 

4D  I 1 xn  ^i'xn'  ) 


(Dl^) 


We  may  define  a universal  curve  applicable  to  any  combination  of  diffusivity 
D and  reaction  rate  K by  letting 


r = r/R 


£ = (x1/R)^Dt 

.KR2 
n « n/~ 


(D15) 

(D16) 

(D17) 


to  give 


n - 4 * * 


ICnJo(xn?)e  J 


(D1S) 


vhe  re 


(D19) 
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Equation  (Dl8)  can  then  be  used  to  generate  curves  of  n against  ? for  various 
values  of  £,  and  these  in  turn  can  be  integrated  to  define  average  number 
densities  as  a function  of  £.  This  has  been  done;  the  results  are  presented 
in  Figure  D2.  Two  averages  are  presented:  the  linear  average 


,1 

(n)  = j n dr  (D20) 

o 

and  the  cylindrical  average 

(n)  =2  r n dr  (D21) 

o 

The  "average"  number  density  seen  by  a probe  laser  passing  along  a diameter  of 
the  cell  with  a finite  beam  width  will  lie  between  these  two  values. 

It  is  expected  that  for  the  flow  tube  experiments  HF  adsorption  on  the 
walls  will  have  little  influence  on  the  results.  These  experiments  are  per- 
formed over  long  time  periods  such  that  the  walls  would  be  saturated  with  HF. 
The  pulsed  fluorescence  experiments,  on  the  other  hand,  are  of  short  duration, 
with  emptying  and  refilling  of  the  test  cell  between  each  run. 

For  the  pulsed-fluorescence  experiments  the  diffusivity  of  the  HF  through 
the  argon  is  D 2 cmd/sec  and  the  radius  of  the  test  cell  is  0.6  cm.  Thus 
we  have 


£ = 3^+t  (D22) 

Consider  now  test  series  A.  The  maximum  possible  HF  is  6.Ch  x 10*D  cm_: ; the 
reaction  occurs  over  approximately  20  milliseconds  at  £ = 0.6h.  By  referring 
to  the  figure  we  see  that  the  average  number  density  is  between  55  and  655$  of 
the  maximum  possible.  This  gives  a number  density  range  of  3*3x1^°  to  3*3x10 
cm~3,  well  below  the  observeu  maxima.  Similar  we  car.  consider  the  other  test 
series.  The  results  are  summarized  in  Table  Dl.  As  can  be  seen,  the 
observed  number  densities  generally  lie  in  a range  for  which  HF  could  have  been 
lost  to  the  walls.  A more  realistic  treatment,  allowing  for  faster  initial 
chemical  rates  and  slower  rates  at  later  times  would  have  given  even  lower 
predicted  number  densities.  Hence  even  with  diffusion  coefficients  as  low 
as  2.0  cm  /sec  the  potential  for  significant  wall  loss  exists. 
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TABLE  Dl.  POSSIBLE  HF  WALL  LOSS  EFFECTS  HI  PUISED-FLUOEESCENCE  TESTS 


Series 

Maximum 
No.  Density 
(cm"3x10-1°) 

Reaction 

Time 

(ms) 

Average  Number  Density 

Observed 

Cylindrical 

(cm-'xl0-1° 

Linear 

( cm'^xlO"1^) 

No.  Densities 
(cm'-'xlO  ^ ) 

A 

6.04 

20 

3-3 

3-9 

5. 2-6.1 

3 

6.6 

20 

3.6 

4.3 

5-9-6. 5 

C 

3-5 

10 

2.4 

2.8 

2. 7-3-3 

D 

2.8 

20-40 

1.6-0. 9 

1.8-1. 2 

1. 5-2.0 

E 

3-2 

10 

2.1 

2.6 

2. 6-3. 2 

7 

1.5 

.. 

20-30 

0.8-0. 7 

1.0-0. 9 

.80-1.4 

* 

Pulse  energies  greater  than  0.1  J 
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LIST  OF  SYMBOLS 

P„  Atomic  fluorine  partial  pressure 

P Molecular  oxygen  partial  pressure 
°2 

k Boltzmann's  constant 

T Temperature  of  gas 

h Planck's  constant 

v ESR  transition  frequency 

8 Bohr  magne  ton 

H Magnetic  field  strength 

f Spectrometer  factor 

X ' Molecular  susceptibility 

Q Theoretical  ESR  intensity  factor 
Ap  Area  bounded  by  F ESR  signal 
Aq  Area  bounded  by  0C  ESP.  signal 
W Width  of  ESR  signal 

L Height  of  ESR  signal 

S ESR  signal  amplification  factor 

V Field  modulation  voltage 
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